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Abstract 

Intensity interferometry and in particular that due to Bose Einstein correla- 
tions (BEC) constitutes at present the only direct experimental method for the 
determination of sizes and lifetimes of sources in particle and nuclear physics. The 
measurement of these is essential for an understanding of the dynamics of strong 
interactions which are responsible for the existence and properties of atomic nuclei. 
Moreover a new state of matter, quark matter, in which the ultimate constituents 
of matter move freely, is within the reach of present accelerators or those under con- 
struction. The confirmation of the existence of this new state is intimately linked 
with the determination of its space-time properties. Furthermore BEC provides 
information about quantum coherence which lies at the basis of the phenomenon 
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of Bose-Einstein condensation seen in many chapters of physics. Coherence and 
the associated classical fields are essential ingredients in modern theories of particle 
physics including the standard model. Last but not least besides this "applica- 
tive" aspect of BEC, this effect has implications for the foundations of quantum 
mechanics including the understanding of the concept of "identical particles" . Re- 
cent theoretical developments in BEC are reviewed and their application in high 
energy particle and heavy-ion reactions is analyzed. The treated topics include: 

a) a comparison between the wave-function approach and the space-time approach 
based on classical currents, which predicts "surprising" particle anti-particle BEC 

b) the study of final state interactions c) the use of hydrodynamics d) the relation 
between correlations and multiplicity distributions. 
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1 Introduction 



The method of photon intensity interferometry was invented in the mid fifties by Hanbury- 
Brown and Twiss for the measurement of stellar dimensions and is sometimes called the 
HBT method. In 1959-1960 G.Goldhaber, S.Goldhaber,W.Lee and A.Pais discovered that 
identical charged pions produced in p — p annihilation are correlated (the GGLP effect). 
Both the HBT and the GGLP effects are based on Bose- Einstein correlations (BEG). 
Subsequently also Fermi-Dirac correlations for nucleons have been observed. Loosely 
speaking both these correlation effects can be viewed as a consequence of the symme- 
try (antisymmetry) properties of the wave function with respect to permutation of two 
identical particles with integer (half-integer) spin and arc thus intrinsic quantum phe- 
nomena. At a higher level, these symmetry properties of identical particles are expressed 
by the commutation relations of the creation and annihilation operators of particles in 
the second quantization (quantum field theory) . The quantum field approach is the more 
general approach as it contains the possibility to deal with creation and annihilation of 
particles and certain phenomena like the correlation between particles and antiparticles 
can be properly described only within this formalism. Furthermore, at high energies, 
because of the large number of particles produced, not all particles can be detected in 
a given reaction and therefore one measures usually only inclusive cross sections. For 
these reactions the wave function formalism is impractical. Related to this is the fact 
that the second quantization provides through the density matrix a transparent link be- 
tween correlations and multiplicity distributions. This last topic has been in the center 
of interest of multiparticle dynamics for the last 20 years (we refer among other things to 
Koba-Nielsen-Olesen (KNO) scaling and "intermittency" ) . Furthermore one of the most 
important properties of systems made of identical bosons and which is responsible for the 
phenomenon of lasing is quantum statistical coherence. This feature is also not accessible 
to a theoretical treatment except in field theory. 

The present review is restricted to Bose Einstein correlations which constitute by far 
the majority of correlations papers both of theoretical and experimental nature. This is 
due to the fact that BEG present important heuristic and methodological advantages as 
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compared to Fermi-Dirac correlations. Among the first we mention the fact that quantum 
coherence appears only in BEC. Among the second, one should recall that pions are the 
most abundantly produced secondaries in high energy reactions. 

In the last years there has been a considerable surge of interest in boson interferometry. 
This can be judged by the fact that at present there is no meeting on multiparticle 
production where numerous contributions to this subject are not presented. Moreover 
since 1990 |l| meetings dedicated exclusively to this topic are organized; this reflects the 
realization that BEC are an important subject in its own. This development is due in part 
to the fact that at present intensity interferometry constitutes the only direct experimental 
method for the determination of sizes and lifetimes of sources in particle and nuclear 
physics. Since soft strong interactions which are responsible for multiparticle production 
processes cannot be treated by perturbative QCD, phenomenological approaches have to 
be used in this domain and space-time concepts are essential elements in these approaches. 
That is why intensity interferometry has become an indispensable tool in the investigation 
of the dynamics of high energy reactions. 

However this alone could not yet explain the explosion of interest in BEC if it were 
not the search for quark-gluon plasma which has mobilized the attention of most of the 
nuclear physics and of an appreciable part of the particle physics community. For several 
reasons the space-time properties of "fireballs" produced in heavy ion reactions are an 
essential key in the process of undertanding whether quark matter has been formed. 

The main emphasis of the present review will be on theoretical developments which 
took place after 1989-1990. Experimental results will be mentioned only in so far as they 
illustrate the theoretical concepts. For a review of older references cf. e.g. the paper 
by Boal, Gelbke, and Jennings 0. In the nineties the single most important theoretical 
event was in our view the space-time generalization of the classical current formalism. For 
a detailed presentation of this generalization and its applications up to 1993 the reader is 
referred to ref.[^. For a review of experimental results in e"*" — e~ reactions cf. Finally 
a more pedagogical and more complete treatment of the theory of BEC can be found in 
i- 

There are two categories of papers not mentioned: (i) Those which the reviewer was 
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unaware off; he apologizes to the authors of these papers for this, (ii) Those which he 
considered as irrelevant or as repetitions of previous workQ. The large number of papers 
quoted despite these restrictions shows that an exhaustive listing of references on BEC is 
not trivial. 

2 The GGLP effect 

In the period 1954-1960 Hanbury-Brown and Twiss developped the method of optical 
intensity interferometry for the determination of (angular sizes) of stars (cf. e.g.[^). 
The particle physics equivalent of the Hanbury-Brown Twiss (HBT) effect in optics is 
the Goldhaber, Goldhaber, Lee and Pais (GGLP) effect 0, |§] which we shall describe 
schematically below. 

However when going over from optics to particle physics the following point has to 
be considered: In particle physics one does not measure distances r in order to deduce 
(differences of) momenta k and thus angular sizes, but one measures rather momenta in 
order to deduce distances. This explains why GGLP were not aware of the HBT method 



Up to a certain point there are two ways of approaching the relationship between 
optics and particle physics: the wave function approach and the field theoretical approach. 
Although the first one is only a particular case of the second one it is still useful because 
it allows sometimes a more intuitive understanding of certain concepts and in particular 
that of the distinction between boson and fermion correlations. We shall start therefore 
with the wave function approach. 

^As a matter of fact, "rediscoveries" are not rare in this field and it would be desirable that journals 
and referees should do a better job in trying to avoid such a situation, not to mention the prevention of 
publication of misleading statements and/or wrong results. 

^For a comparison of optical and particle physics intensity interferometry cf. e.g.|^]. 
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2.1 The wave function approach 

Let us consider for the beginning a source which consists of a number of discrete emission 
points i each of which is characterized by a probabihty amphtude 

Fi{r) = FiS{r-ri). (2.1) 

Let ipkir) be the wave function of an emitted particle. The total probability -P(k) 
to observe the emission of one particle with momentum k from this source is obtained 
by summing the contributions of all points i. If this summation is done incoherently the 
single particle probability reads 

P{k)=J:\FMr^)\' (2.2) 

i 

Instead of discrete emission points consider now a source the emission points of which are 
continuously distributed in space and assume for simplicity that the wave functions are 
plane waves ^k(i") ^ exp{ikr). The sum over Xi will now be replaced by an integral and 
we have 

P(k) ^ J \F{r)\^£r. (2.3) 
Similarly the probability to observe two particles with momenta ki, k2 is 

P(kik2) = I |Vi,2r/(ri)/(r2)dVidV2 (2.4) 

Here ■0i^2 = '0i,2(ki, k2; ri, r2) is the two-particle wave function and we have introduced 
the density distribution / = 

Suppose now that the two particles are identical. Then the two-particle wave function 
has to be symmetrized or antisymmetrized depending on whether we deal with identical 
bosons or fermions. Assuming plane waves we have 

V2 

with the plus sign for bosons and the minus sign for fermions. With this wave function 
one obtains 

P(ki,k2) = |/,(0)r±|/,(ki-k2)r (2.6) 
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The incoherent summation corresponds to random fluctuations of the amphtudes Fi. 

Following the GGLP experiment [Q consider the space points ri and r2 within a 
source, so that each point emits two identical particles (equally charged pions in the case 
of GGLP) with momenta ki and k2. These particles are detected in the registration 
points Xi and X2 so that in Xi only particles of momentum ki and in X2 only particles of 
momentum ^2 are registered (cf.Fig.l). Because of the identity of particles one cannot 
decide which particle pair originates in ri and which originates in r2. 

Assuming that the individual emission points of the source act incoherently GGLP 
derived eq.( |2.6D which for bosons leads to the second order correlation function 

C2(ki,k2) = P2(kl,k2)/Pi(ki)Pl(k2) = 1 + |/(q)P (2.7) 

where / is the Fourier transform of / and q the momentum difference ki — k2. 

Eq. ( |2.7| ) shows quite clearly how in particle physics momentum (correlation) measure- 
ments can yield information about the space-time structure of the source 

In the case of coherent summation one gets instead C2 = 1 (cf. e.g. 0). We retain 
herefrom (cf. also below) that coherence reduces the correlation and that a purely coherent 
source has a correlation function which does not depend on its geometry^. 

2.1.1 Newer correlation measurements in p — p annihilation 

The original GGLP experiment [0] measured the correlation function in terms of the 
opening angle ^ of a pion pair. 



The GGLP experiment has been repeated in the last years at LEAR (cf.e.g. [12] 
where also older references are quoted Q). In these newer experiments the three momenta 
of particles were measured, although one continued to use as a variable for the correlation 
function the invariant momentum difference Q, as suggested already in P]. 

•^Cf. also 1^,(1^ for an attempt to approach the issue of coherence and chaos by using wave packets. 
^It may be amusing to note that the use of the opening angle as a kinematical variable in BEC studies 

was rediscovered in 1994 in ref. where it was recommended as a tool for the investigation of final 

state interactions. 

^For a reanalysis of the results of an older annihilation experiment cf. [|l3|. 



10 



One of the remarkable observations made in all annihilation reactions is the fact that 
the intercept of the second order correlation function C2{k,k) appears consistently to 
exceed the canonical value of 2 reaching values up to 4. (This effect was possibly not 
seen in the original experiment [0] because of the averaging over the magnitudes of the 
momenta.) 

In this effect was attributed to resonances while subsequently |]r^,||T^ a (non- 



chaotic) Skyrmion type superposition of coherent states was proposed as an alternative 
explanation. 

Another possible explanation of this intriguing observation may be that in annihilation 
processes squeezed states (cf. section 2.2) are produced, while this is not the case in other 
processes. Indeed, it is known that squeezed states can lead to overbunching effects 
(cf. below). Furthermore, as shown in ref.([|TB|), squeezed states may be preferentially 



produced in rapid reactions. Or, according to some authors |T^, annihilation is a more 
rapid process than other reactions occuring at higher energies. Given the importance 
of squeezed states, further experimental and theoretical studies of this issue are highly 
desirable f\ 

2.1.2 Resonances, apparent coherence and other experimental problems; the 
A factor 

It is known that in multiparticle production processes, an appreciable fraction of pions 
originate from resonance decays. Resonances act in two opposite ways on the correlation 
function of pion pairs. On one hand the interference between pions originating from short 
lived resonances and "directly" produced pions leads to a narrow peak in the second 
order correlation function at small values of q [|T7|, \TE\. On the other hand long lived 



resonances give rise to pions which are beyond the range of detectors and this leads to 
an apparent decrease of the intercept of the second order correlation function C2{k,k). 
These modifications of the intercept are important among other things, because their 



^In ref. the effect of chaotic superpositions of squeezed states in BEC was also studied. It is 
shown that such a superposition leads always to overbunching, while pure squeezed states can lead also 
to antibunching. 
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understanding is essential in the search for coherence through the intercept criterion. As 
mentioned aheady one of the most immediate consequences of coherence is the decrease 
of the correlation function at small q. 

Because of the large number of different resonances produced in high energy reactions a 
quantitative estimate of their effects is possible only via numerical techniques. In sections 
4.9.1 and 5.1.3 we will present a more detailed discussion of the influence of resonances 
on BEC within the Wigner function formalism. For older references on this topic cf. also 

i, [0, m- 

Another experimental difficulty is that, because of limited statistics or certain techni- 
cal problems, sometimes not all degrees of freedom can be measured in a given expriment. 
This leads to an effective averaging over the non-measured degrees of freedom and hence 
also to an effective reduction of the correlation function. As a matter of fact BEC ex- 
periments have shown from the very beginning that the extrapolation of the correlation 
function to g = almost never led to the maximum value 6*2(0) = 2 permitted by eq.( p.7| ). 

To take into account empirically this effect experimentalists introduced into the cor- 
relation function a correction factor A. Thus eg. ( p.7|) e.g. was modified into 



C2(ki,k2) = l + A|/(q)p (2.8) 

Because one initially thought that formally this generalization offers also the possibility 
to describe partially coherent sources, the corresponding parameter A was postulated to 
be limited by (0,1) and called the "incoherence" factor: indeed A = leads to a totally 
coherent source and A = 1 to a totally chaotic one. Unfortunately this nomenclature 
is not quite correct, as will be explained below. It should also be mentioned that there 
exists a strong correlation between the empirical values of A and and those of the "radius" 
which enters / (cf. e.g. |2T| for a special study of this issue.) 



2.1.3 The limitations of the wave function formalism 

The wave function formalism presented in the previous subsection has severe shortcom- 
ings: 
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a) The correlation function ( p.7|) depends only on the difference of momenta q and not 
also on the sum ki + k2, in contradiction with experimental data. Although we will see in 
section 4.3 that this limitation disappears automatically in the quantum statistical (field 
theoretical) approach, it can be remedied also within the first quantization formalism 
by using the Wigner function approach. In this approach the source function is from 
the beginning a function both of coordinates and momenta and therefore the correlation 
function depends on ki and k2 separately. There is of course a price to pay for this 
procedure as it involves a semiclassical approximation []. 

b) The wave function formalism may be useful when exclusive reactions are considered 
as was the case e.g. in the GGLP paper. Indeed ipi^2 in eq. (p^) is just the wave function 
of the two-boson system i.e. the assumption is made that two and only two bosons are 
produced. At low energies, i.e. low average multiplicities, this condition can be satisfied. 
However at high energies the identification of all particles is very difficult and up to now 
has not been done. Therefore one measures practically always inclusive cross sections. 
This means that instead of (|2.5| ) one would have to use a wave function which describes 
the two bosons in the presence of all other produced particles. To obtain such a wave 
function one would have to solve the Schrodinger equation of the many body, strongly 
interacting system, which is not a very practical proposal. Related to this is the difficulty 
to treat higher order correlations within the wave function formalism. 

c) The correlation function C2 in the wave function formalism is independent of isospin 
and is thus the same for charged and neutral particles. We shall see in section 4.4 that in 
a more correct quantum field theoretical approach, this is not the case. This will affect 
among other things the bounds of the correlations and will lead to quantum statistical 
particle-antiparticle correlations which are not expected in the wave function formalism. 

d) Coherence cannot be treated adequately (cf. below). The correlation functions 
derived in this subsection refer in general to incoherent sources and attempts to intro- 
duce coherence within the wave function formalism are rather ad-hoc parametrizations. 

^ String models |^j2^ also use a Wigner type formalism. Here it is postulated that there exists a 
"formation" time r and therefore the particle production points are distributed around — — t'^. 
This implies among other things a correlation between particle production points and momenta. 
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However coherence is the most characteristic and important property of Bose-Einstein 
correlations among other things because it is the basis of the phenomenon of Bose Ein- 
stein condensates found in many chapters of physics, hke superconductors, superfluids, 
lasers, and the recently discovered atomic condensates [^. It would be very surprising if 
coherence would not be found also in particle physics given the fact that the wave lengths 
of the emitted particles are of the same order as that of the sources. Furthermore as 



pointed out in this connection in modern particle physics is based on spontaneously 
broken symmetries. The associated fields are coherent. That is why one of the main moti- 
vations of BEC research should be the measurement of the amount of coherence in strong 
interactions. For this purpose the formalism of BEC has to be generalized to include the 
presence of (partial) coherence and this again can be done correctly only within quantum 
statistics i.e. quantum field theory. 

We conclude this subsection with the observation that the wave function or wave 
packet approach may nevertheless be useful in BEC for the investigation of final state 
interactions (cf. below) or for the construction of event generators, where phases or 
quantum amplitudes are ignored anyway. Also for correlations between fermions where 
coherence is absent the wave function formalism may be an adequate substitute, although 
here too a field theoretical approach is possible. 



2.2 Quantum optical methods in BEC 

In high energy processes in which the pion multiplicity is large, we may in general expect 
the methods of quantum statistics (QS) f\ to be useful. For BEC in particular they turn 
out to be indispensable. These methods have been applied with great success particularly 
in quantum optics (QO), superfluidity, supercoductivity etc. What distinguishes optical 
phenomena from those in particle physics are conservation laws and final state interactions 
which are present in hadron physics. At high energies and high multiplicities the first are 
unimportant. Neglecting for the moment also the final state interactions, QS reduces 
then to QO and we may take over the formalism of QO to interpret the data on multipion 
^By QS we understand in the following the density matrix formalism within second quantization. 
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production at high energies, provided we consider identical pions. Given the general 
validity of QS (or QO), it is then clear that any model of multiparticle production must 
satisfy the laws of quantum statistics and this has far reaching consequences, independent 
of the particular dynamical mechanism which governs the production process. The main 
tools in the QO formalism are defined below |. 

Coherent states and squeezed coherent states Coherent states \a > are eigenstates 
of the (one-particle) annihilation operator a 

a\a >= a\a > . (2.9) 

Squeezed coherent states are eigenstates |/5 of the two-particle annihilation operator 

b = fia + ua'^ (2.10) 

with 

_ |zy|2 = 1 (2.11) 

SO that 

. (2.12) 

One of the remarkable properties of these states which explains also their name is that 
for them the uncertainty in one variable can be squeezed at the expense of the other so 
that 

{Aq)l < ^, {Apr, > |, (2.13) 

or vice versa. The importance of this remarkable property lies among other things in 

the possibility to reduce quantum fluctuations and this explains the great expectations 

associated with them in communication and measurement technology as well as their 

interest from a heuristic point of view. 

It has been found recently in [0 that squeezed states appear naturally when one deals 

with rapid phase transitions (explosions). Indeed consider the transition from a system 
review of the applications of quantum optical methods to multiparticle production up to 1988 can 
be found in p6| . 
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a to a system b and assume that it proceeds rapidly enough so that the relation between 
the creation and annihilation operators and the corresponding fields in the two "phases" 
remains unchanged. Mathematically this process will be described by postulationg at the 
moment of this transition the following relations between the generalized coordinate q 
and the generalized momentum p of the field: 




2 




(b^ - b) = tJ ^{a^ - a) (2.14) 



a\ a 



are the free field creation and annihilation operators in the "phase a" and b^,b the 



corresponding operators in the "phase 6" . Eq. (|2.14| ) holds for each mode p. Then we 



get immediately a connection between the a and b operators, 

a=b cosh r + b^ sinh r , 

a^=b sinh r + b^ cosh r (2-15) 

with 

r = r{p) = Uog{EjE,). (2.16) 



The transformation ( |2.15| ) is just the squeezing transformation ( |2.10| ) with 



/i = coshr; i/ = sinhr (2-17) 

which proves the statement made above. 

The observation of squeezed states in BEC may thus serve as a signal for such rapid 
transitions. Furthermore the existence of isospin induces in hadronic BEC certain effects 
which are specific for squeezed states. This topic will be discussed in section 4.4. 

IFtoiu. the point of view of BEC what distinguishes ordinary coherent states from 
squeezed states is the following: for coherent states the intercept C2{k,k) = 1 while for 
squeezed states it can take arbitrary values. In Fig.2 from ref.||2^ one can see such an 
example. 
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Expansions in terms of coherent states Coherent states form an (over) complete set 
so that an arbitrary state \f > can be expanded in a unique way in terms of these states. 

Of particular use is the expansion of the density matrix p in terms of coherent states. 
For a pure coherent state the density operator reads 

p = \a >< a\. (2-18) 

For an arbitrary density matrix case we have 



P 



j V{a)\a><a\d?a. (2.19) 



Here P is a weight function which usually, but not always, has the meaning of a proba- 
bility. The normalization condition for the density operator translates in terms of the V 
representation as 

trp = J V{a)d?a = 1. (2.20) 
Eq.( |2.19D is called also the Glauber-Sudarshan representation. 



The knowledge of V{a) is (almost) equivalent to the knowledge of the density matrix. 
However in most cases the exact form of V{a) is not accessible and one has to content 
oneself with certain approximations of it. Among these approximations the Gaussian form 
is priviledged because: 

(i) one can prove that V{a) is of Gaussian form for a certain physical situation which 
is frequently met in many body physics. 

(ii) its use introduces an enormous mathematical simplification. 
Proposition (i) is the subject of the central limit theorem which states that if 

1. the number of sources becomes large; 

2. they are stationary in the sense that their weight function V{a) depends only on 
the absolute value 

3. they act independently, 

then V{a) is Gaussian. These conditions are known to be fulfilled in most cases 
of optics and presumably also in high energy physics. Chaotic fields and in particular 
systems in thermal equilibrium are described by a Gaussian density matrix. 
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One of the reasons why the Gaussian form for V plays such an important part in 
correlation studies is the fact that for a Gaussian 'P(a) all higher order correlations can 
be expressed in terms of the first two correlation functions. 

On the other hand the coherent state representation is particularly important for 
correlation studies because in this representation all correlation functions can be expressed 
in terms of the creation and annihilation operators a and of the fields (particles). This 
follows from the Fourier expansion of an arbitrary field in second quantization 

vr(x) = ^Ke-*^^ + 4e''="]. (2.21) 

k 

This property will be used extensively in section 4.2 within the classical current for- 
malism. 

Correlation functions The first order correlation function reads 

G^^\x, x) = Tr[pn\x)n{x')] (2.22) 
Higher (n-th) order correlation functions are defined analogously by 

G'^"'\xi...Xn,Xn+l...X2n) = Tr[p7C^ (xi) . . .Tt'' {Xn)n{Xn+l) .7^{x2n)] (2.23) 

In quantum field theory because of the mathematical complexity of the problem exact 
solutions of the field equations are available only in special cases. One such case will 
be discussed later on. However for strong interactions Q even for this case one has to 
use phenomenological parametrizations of the correlation functions and determine the 
parameters (which have a definite physical meaning) by comparing with experiment. 

In optics for stationary chaotic fields two particular parametrizations are used: 

l.Lorentzian spectrum: 

G(i)(xi, X2) =< n,h > e-l"i-"^l/« (2.24) 
^"^Strong interactions are present not only in hadronic physics but also in quantum optics. 
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2. Gaussian spectrum: 

G^^\xi,X2) =< n,h > e-l"^-"^l'/«' (2.25) 

^ is the coherence length in x-space and < rich > is the mean number of particles 
associated with the chaotic fields. 



In it was proposed to use the analogy between time and rapidity in applying the 
methods of quantum optics to particle physics. Indeed in optics processes are usually 
stationary in time while in particle physics the corresponding stationary variable (in the 
rapidity plateau region) is rapidity 

Pure coherent or pure chaotic fields are just extreme cases. In general one expects 
partial coherence i.e. a superposition of coherent and chaotic fields 

^ '^coherent ~l~ '^chaotic- (2.26) 

This leads for the Lorentzian case e.g.to a second order correlation function of the 
form 

C2{x, x') = l + 2p{l - p)e-l^-^'l/« + p2g-2|x-x'|/^ (2.27) 

where p is the chaoticity, which varies between (for purely coherent sources) and 1 (for 
totally chaotic sources). Eq. (|2.8|) is seen to be a particular form of the above equation 
for X = p = 1 and it is clear herefrom that A does not describe (partial) coherence as its 
name would imply. The presence of coherence introduces a new term into C2. However it 
is remarkable that the number of free parameters in eq.(|2.27|) is the same as in eq.(^.8D. 



Formally it appears as if there would act two sources rather than one, but the "weights" 
and the space-time characteristics of these two sources are in a well defined relationship. 



This circumstance had been forgotten up to 1989 ||2^ both by experimentalists and 
theorists. The reason for this collective amnesia is the fact that during the eighties the 
wave function formalism was dominating the BEG literature, especially the experimental 
one. 

/^From the foregoing discussion it should be clear that there are various reasons besides 
coherence why the bunching effect in BEG is reduced. However it should also be clear that 
^^This property does not hold for other variables. 
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the empirical description of this state of affairs through the A factor is possible only for 
totally chaotic sources. Since in an experiment this is never known a priori, this implies 
that the fitting of data with a formula of this type is misleading and should be avoided, 
the more so that the correct formula (|2.27|) does not contain more free parameters than 
eq.(|2.8|). In the example presented above these free parameters are p and ^ for eq. ( |2.27|) 
and A and the effective radius (which enters in /) for eq.( |2.8| ) respectively. 

Using the rapidity-time analogy of ref.[^| for a partially chaotic field the second order 
correlation function in rapidity is then given by eq. ( p.27|) , with x being replaced by rapidity 

y- 

The fact that the last two terms in eq.( p.27D are in a well defined relationship and 
depend in a characteristic way on the two parameters p and ^ is a consequence of the 
superposition of the two fields (coherent and chaotic) and distinguishes a partially coherent 
source from a source which is a superposition of two independent chaotic intensities. 



Because of this the form ( |2.27| ) was proposed in p8| to be used as a signal for detection 
of coherence in BEC Q. As a matter of fact, an attempt in this direction was done in an 



experimental study by Kulka and Lorstad p9[. In this analysis BEC data from pp and pp 
reactions at ^/s = 53GeV were used to compare various forms of correlation functions. 
Among other things one considered formulae of QO type for rapidity 



C2 = 1 + 2p{l - p)e-\y'-y^\/^ + p2g-2|3/i-3/2l/? (2.28) 
(corresponding to a Lorentzian spectrum) and 

C2 = 1 + 2p{l - p)e-\y-^-y^\yi' + p2g-2|j/i-j/2lVC^ (2.29) 

(corresponding to a Gaussian spectrum) as well as arbitrary superpositions of two chaotic 
sources of exponential or Gaussian form respectively. 

C2 = 1 + Aie-l^i-^2l/« + X^e-^ly^-y^l/^ (2.30) 
^^Eq. ( ^.27|) is a special case of superposition of coherent and chaotic fields; it can be considered 
as corresponding to point-like coherent and chaotic sources and a momentum independent chaoticity; 
superpositions of more general finite-size sources are considered in section 4.4 (eq. [f.4C| ) and section 4.8. 
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(2.31) 



Here Ai and A2 represent arbitrary weights of the two chaotic sources. 

Because of the hmited statistics no conclusion could be drawn as to the preference 
of the QO form versus the two-source form. Similar inconclusive results were obtained 
when one replaced in the above equations yi — y2 by the invariant momentum difference 

2.2.1 Higher order correlations 

We have mentioned above that a characteristic property of the Gaussian form of density 
matrix (not to be confused with the Gaussian form of the correlator or the Gaussian form 
of the space-time distribution) is the fact that all higher order correlation functions are 
determined just by the first two correlation functions. Since all BEG studies in particle 
physics performed so far assume a Gaussian density matrix, the reader may wonder why 
it is necessary to measure higher order correlation functions. 
There are at least three reasons for this: 

(i) The conditions of the applicability of the above theorem and in particular the 
postulate that the number of sources is infinite and that they act independently can 
never be fulfilled exactly. 

(ii) In the absence of a theory which determines from first principles the first two 
correlation functions, models for these quantities are used, which are only approximations. 
The errors introduced by these phenomenological parametrizations manifest themselves 
differently in each order and thus violate the above theorem even if (i) would not apply. 

(iii) In experiments, because of limited statistics and sometimes also because of theo- 
retical biases not all physical observables are determined, but rather averages over certain 
variables are performed, which again introduce errors which propagate (and are amplified) 
from lower to higher correlations. 

Gonversely, by comparing correlation functions of different order one can test the appli- 
cability of the theorem quoted above and pin down more precisely the parameters which 
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determine the first two correlation functions (e.g. tlie cfiaoticity p and tfie correlation 
length ^ in eqs. ( [^.28D , ( |2.29| )), which is essentially the purpose of particle interferometry. 

The phenomenological application of these considerations will be discussed in the 
following as well as in section 6.3 for the particular case of the invariant Q variable, but 
the arguments (i),(ii),and (iii) have general validity. It would be a worthwhile research 
project to compare the deviations introduced in the relation between lower and higher 
correlation functions, due to (i) with those introduced by (ii) and (iii). 

The simplification brought by the variable Q can be enhanced by a further approxi- 



mation proposed by Biyajima et al. |^0|- With the notation Qij = ki — kj the analogue 
of eq. ( p. 29 ) can be written 

C2 = l + 2p{l - p) exp{-R^Ql^ + p^ exp{-2Ql2R^) (2.32) 
For the third order correlation function one obtains 

C,=l + 2p{l - p)[exp{-R^Ql,} + exp{-R^Ql,} + exp{-R'Ql,}] 
V[exp{-2i?2g2j ^ exp{-2i?2g2j + exp{-2R^Ql^}] 

+2p\l - p)[exp{-R\Ql, + Ql,)} + exp{-R'{Ql, + Ql,)} 

+ exp{-R'iQl, + Ql,)}] + 2p'[exp{-R'iQl, + Ql, + Ql,)}] (2.33) 
Ref. 10 proposed to use symmetrical configurations for all two-particle momentum differ- 



ences i.e. to consider Qij independent of (z, j). For C2 this assumption does not of course 
introduce any modifications. However for higher order the simplification is important. 
Thus e.g. for the third order correlation with Q12 = Qa = Q23 and with the definition 
Qthree = Qu + Qls + Qls eq.(g;33|) becomes 

^3=1 + Qp{l - p) exp |-li?2g2^^^^| 

+3p2(3 - 2p) exp |-^i?2g2^^^J ^ 2p3exp{-i?2Q2^_} (2.34) 



In ref.[0 similar expressions for C4 and C5, again for a Gaussian correlator, were 
given. 
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These relations for higher order BEC were subjected to an experimental test in ref. ^ 
using the UAl data for pp reactions at ^/s = 630 and dOOGeV. For reasons which will 
become clear immediately, we discuss here this topic in some detail. 

The procedure used in |^T|] for this test consisted in determining R and p separately 
for each order q of the correlation and comparing these values for different q. It was found 
that a Gaussian correlator did not fit the data. Next in one tried to replace the 



Gaussian correlator by an exponential (cf. eq.( |2.28| ). To do this one substituted simply 
in the expressions for the correlation functions of ref. the factor exp(— with 
exp{~RQ). Such a procedure was at hand given the fact that for C2 the QO formulae 
both for an exponential correlator and a Gaussian correlator were known and their 
comparison suggested just this substitution, as can be seen from eq.(p^.29[). 



In |31] one used then for the exponential correlator the relations 



Cr"'""^' = 1 + 2p(l - p) expi-RQu) + p" exp(-2gi2i?) 



(2.35) 



and 



^empirical^^ + 6p(l - p) CXp (^-^^RQthre^j 

+3p'(3 - 2p) exp [-'^RQthre)j + 2p' exp {-RQthree) (2.36) 



With these modified formulae one still could not find in M^ a unique set of values 



p and R for all orders of correlation functions. However now a clearer picture of the 
"disagreement" between the QO formalism and the data emerged. It seemed that while 
the parameter p was more or less independent of q, the radius R increased with the order 
g in a way which could be approximated by the relation 



R, = Rj\q{q~l). (2.37) 



^•^This reminds one of Polonius's bewilderment in Shakespeare's Hamlet: "Though this be madness, 
yet there is method in't" . 
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Indeed in it was shown that the findings of and in particular eq. (|2.37| ) not 
only did not contradict QS but on the contrary constituted a confirmation of it. While 



eq. (|2.35|) for the second order correlation function coincides with that derived in quantum 
optics for an exponential spectrum, this is not the case with the expressions for higher 
order correlations C^'^^^^^'^"'^ (eq.( |2.36] )). The formulae for C3, C4 and C5 corresponding in 
QS to an exponential correlator and derived in differ from the empirical ones used in 
|31| . As an example we quote 



C3 = l + 6p(l - p) exp (^--^RQthre^ 

+3p\3 - 2p) exp l^-^RQ^j^,^)^ + 2p=' eM-^RQthree) (2.38) 

As observed in |^ and as one easily can check by comparing eq. ( |2.36| ) with eq. ( |2.38| ), 
one can make coincide the empirical formulae for Cq used in with the correct ones, by 
replacing the parameter R with a scaled parameter Rg and the relation between R and i?^, 
is, nothing else but Rs = Rg, where Rg is given by eq. ( ^.371 ) . The fact that this happens 
for three different orders, i.e. for C3, C4, and C5 makes a coinicidence quite improbable. 

By empirically modifying the formulae of higher order BEC for the exponential case, 
paper |31] had explicitely violated QS and the "phenomenological" relationship ( 2.37|) 
between R and q just compensated this violation. 

The fact that this compensation and the final agreement between theory and experiment was not 
perfect is not surprising and is discussed in p3]. Besides the reasons (i), (ii), (iii) mentioned above, one 



has to take into account that the QO formahsm on which eqs.( 2.35 ),( p.33| ) are based assumes stationarity 
in Q, i.e. assumes that the corrrelator depends only on the difference of momenta fci — k2 and not also 
on their sum. As mentioned already this condition is in general not fulfilled in BEC. Furthermore, the 
parameter p, if it is related to chaoticity, is in general momentum dependent (cf. section 4.8). Also, the 
symmetry assumption, Q{i,j) independent of may be too strong. Besides these theoretical caveats, 
there are also experimental problems, related to the fact that the U Al experiment is not a dedicated 
BEC experiment and thus suffers from specific diseases, which are common to almost all particle physics 
BEC experiments performed so far. Among other things, there is no identification of particles (only 85% 
of the tracks recorded are pions), and the normalization of correlation functions is the "conventional" 
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one, i.e. not based on the single inclusive cross sections as the definition of correlation functions demands 



(cf.(2.7) and section 4.11), but rather uses an empirically determined "background" ensemble. 

The very fact that QS was effectively "reinvented" by an unbiased experimental team 
is in our view overwhelmingly convincing. It may enter the history of physics as a classical 
and quite unique example of the intimate relationship between theory and experiment. 
Ref.p2| concludes that experiment |^ supports QS in general and the standard form 
of the density matrix in particular. From an epistemological point of view, the mes- 
sage of this "conspiracy" is highly significant: QS is robust enough to resist attempts of 
falsification. 

In the mean time further theoretical and experimental developments took place. 

On the theoretical side a new space-time approach to BEG was developped 0, 
which is more appropiate to particle physics and which contains as a special case the 
QO formalism. In particular the two exponential feature of the correlation function is 
recovered. 

On the experimental side a new technique for the study of higher order correlations 
was developed, the method of correlation integrals which was applied ||3^ to a subset of 
the UAl data in order to test the above quoted QO fromalism. The fits were restricted 
to second and third order cumulants only. Again it was found that by extracting the 
parameters p and R from the second order data, the "predicted" third order correlation, 
this time by using a correct QO formula, differed significantly from the measured one. 

If confirmed, such a result could indicate that the QO formalism provides only a 
rough description of the data and that higher precision data demand also more realistic 
theoretical tools. Such tools are the QS space-time approach to BEG presented in section 
4.3. A further, but more remote possibility would be to look for deviations from the 
Gaussian form of the density matrix. However it seems premature to speculate along 
these lines given the fact that the procedure used to test the relation between the second 



and third order correlation functions may have to be qualified. Indeed in one did not 
perform a simultaneous fit of second and third order data to check the QO formalism. Such 
a simultaneous fit appears necessary before drawing conclusions, because as mentioned 
above (cf.(ii) and (iii)), the errors involved in "guessing" the form of the correlator, and 
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the fact that the variable Q does not characterize completely the two-particle correlation, 
limit the applicability of the theorem which reduces higher order correlations to first 
and second ones. As a matter of fact, it was found (cf. also section 6.3) in a 
comparison of the QS space-time approach with higher order correlation data, that the 
second order correlation data is quite insensitive to the values of the parameters which 
enter the correlator, while once higher order data are used in a simultaneous fit, a strong 
delimitation of the acceptable parameter values results. Thus there are several possible 
solutions if one restricts the fit to the second order correlation and the correct one among 
those can be found only by fitting simultaneously all correlations. If by accident one 
chooses in a lower correlation the wrong parameter set, then the higher correlations cannot 
be fitted anymore^. 

Before ending these phenomenological considerations in which the variable Q played 
a major part, a few remarks about its use may be in order. 

The invariant momentum difference Q BEC studies in particle physics use often a 
priviledged variable namely the squared momentum difference 

q2 = (fc^ _ = (k^ _ k2)2 - (El - E2)'. (2.39) 

It owes its special role to the fact that it is a relativistic invariant and was used already 
in the pioneering paper by GGLP ||^. It has also the advantage that it involves all four 
components of the momenta k simultaneously so that the intercept of the correlation 
function C2(k, k) coincides with C2{Q = 0). Thus by measuring C2 as a function of one 
single scalar quantity Q one gets automatically the intercept. This is not the case with 
other single scalar quantities used in BEC like yi — y2 or kj_^i — k±^2 which characterize 
only partially the intercept. 

On the other hand, Q suffers from certain serious diseases which make its use for 
practical interferometrical purposes questionable. 

^^The fact that in [34| the correlations were normahzed by mixing events rather than by comparing 
with the single inclusive cross sections in the same event, as prescribed by the definition of the correlation 
functions may also influence the applicability of the central limit theorem. 
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The first and most important deficiency of Q is the fact that it mixes time and space 
coordinates: the associated quantity R in the conventional parametrization of the correla- 
tion function C2 = 1 + A exp(— is neither a radius nor a lifetime, but a combination 
of these, which cannot be easily disentangled. Another deficiency of Q, which is com- 
mon to all single scalar quantities is the circumstance that it does not fully characterize 
the correlation function. Indeed the second order correlation function C2 is in general a 
function of six independent quantities which cannot be replaced by a single variable. 

An improvement on Q was proposed by Cramer with the introduction of coa- 
lescence variables which constitute a set of three boost invariant variables to replace, 
for a pair of identical particles, the single variable Q. They are related to Q by = 
2rn?{Cl + C|. + where C^, Ct, Cr denote longitudinal, tranverse and radial coales- 
cence respectively, and m the mass of the particle. They have the properties that (7^ = 
when ui = y2, Ct = when 0i = 02 and Cr = when either mi = m2 or ki = Here 
rrii is the transverse mass and (p the azimuthal angle in the transverse plane. It is shown 



m 



^1 that with these new variables a Lorentz invariant separation of the space-like and 
time-like characteristics of the source is possible, within the kinematical assumptions in- 
volved by the particular choice of the coalescence variables. This separation is however 
rather involved. In the coalescence variables are used for the introduction of Coulomb 
corrections into second and higher order correlation functions. 

Another way to compensate in part for the fact that one single variable does not 
characterize completely the two-particle system, but which maintains the use of Q is, as 
was explained above, to consider higher order correlations. 



3 Final state interactions of hadronic bosons 

One of the most important differences between the HBT effect in optics and the corre- 
sponding GGLP effect in particle physics is the fact that in the first case we deal with 
photons while in the second case with hadrons. 

While photons in a first approximation do not interact, hadrons do interact. This 
interaction has two effects: (i) it influences the correlation between identical hadrons and 
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(ii) it leads to correlations also between non-identical hadrons. 

This review deals only with correlations due to the identity of particles and in partic- 
ular with Bose-Einstein correlations. Therefore only effect (i) will be discussed 

Effect (i) is usually described in terms of final state interactions. In some theoretical 
studies (cf. e.g. and references quoted there) emission of particles at different times 
is also treated as an effective final state interaction. /^From the BEC point of view the 
final states interaction constitutes in general an unwanted background, which has to be 
substracted in order to obtain the "true" quantum statistical effect on which interferom- 
etry measurements are based. That this is not always a trivial task will be shown in the 
following. 

There are two types of final state interactions in hadronic interferometry: Electromag- 
netic, traded under the generic name of Coulomb interactions, and strong. Furthermore 
one distinguishes between one-body final state interactions and many body final state 
interactions. 

3.1 Electromagnetic final state interactions 

The plane wave two-body function used in the considerations above (cf. eq.( |2.5| )) applies 
of course only for non-interacting particles. 

As a first step towards a more general treatment consider charged particle interfer- 
ometry. As a matter of fact the vast majority of BEC studies, both of experimental and 
theoretical nature, refer to charged pions. For two-particle correlations, we will have to 
consider the interaction of each member of a pair with the charge of the source and the 
Coulomb interaction between the two particles constituting the pair. The first effect will 
affect primarily the single particle probabilities and is not expected to depend on the mo- 
mentum difference q. Attention has been paid so far mostly to the second effect, i.e the 
modification of the two-particle wave function due to the Coulomb interaction between 
the two particles. While initially, having in mind the Gamow formula, it was assumed 
^^The reader interested in correlations between non-identical particles is refered to |^ for the period 
up to 1990; for more recent literature cf. ||37|] , where correlations in low energy heavy ion reactions are 
reviewed. 
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that this effect is (for small q values) quite important, at present serious doubts about 
these estimates have arisen. The model dependence of corrections for this effect makes 
it almost imperative that experimental data should be presented also without Coulomb 
corrections, so that it should be left to the reader the possibility of introducing (or not 
introducing) corrections according to her/his own prejudice. 

3.1.1 Coulomb correction and "overcorrection" 

As usual one separates the center of mass motion from the relative motion. For the last 
one the scattering wave function reads: 



uir) e'"' + r-'g{e, <P)e"'' (3.40) 

where the relative position vector r has polar coordinates r, 6, 0. 

The form of the function g depends on the scattering potential. In the Coulomb case 
the corresponding Schrodinger equation can be solved exactly and the correction to the 
two-particle wave function and the correlation function can be calculated. So far three 
more and more sophisticated procedures were used for this purpose. In ||39| the value of 
the square modulus of the wave function u in the origin r = was proposed as a correction 
term G to the correlation function C2. Up to non- interesting factors this is the Gamow 
factor which reads 

i"(°'i^= exp(2:;')-i =g""- f^-^) 

Here 

rj = arriT^/q (3.42) 

and = 1^1 — ^21- ^-n is the pion mass. 



However as pointed out by Bowler and subsequently also by others, in BEC this 
approximation may be questionable. Indeed in a typical — e~ reaction e.g. the source 
which gives rise to BEC has a size of the order of 1/m which is a large number compared 
with a typical "Coulomb length" defined as the classical turning point where the kinetic 
energy balances the potential Coulomb energy: 

gV2m,ed = eVr-i (3.43) 
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where q = q/2 and rrired is the reduced mass of the pion pair. For a typical BEC momen- 
tum difference of g = lOOMeV^, one gets from (|3.43|) rt = 0.08/m. This suggests that by 



taking the value of the wave function at r = one overstimates the Coulomb correction. 

More recently Biyajima and collaborators (cf. also |^) have considered a further 
correction to the correction proposed by Bowler, which decreases even more the Coulomb 
effect and which is also of heuristic interest. In it is pointed out that the wave function 
( p.40| ) used by Bowler does not yet take into account the symmetry of the two-particle 
system. It has to be supplemented by an exchange term so that the rhs of ( |3.40| ) becomes 

^^igz ^ ^-igz^ ^ ^j^Q^ + J(7r - 0, + TT)]r-h''}' (3.44) 



The corrections due to this new effect are of the same order as those found in |50] and 
go in the same direction. 

Another approach to the Coulomb correction in BEC has been suggested by Baym and 



Braun-Munzinger |f43|| . Starting from the observation of |^ about the classical turning 
point these authors propose the use for heavy ion reactions of a classical Coulomb correc- 
tion factor arising in the assumption that the Coulomb effect of the pair is negligible for 
separations less than an initial radius tq. This model is tested by comparing its results 
with experimental data on tt+tt^, vr^p, n^p correlations in heavy ion collisions^ (Au-Au at 
AGS energies). The assumption behind this comparison is that the observed correlations 
are solely due to the Coulomb effect Indeed qualitatively this seems to be the case: 
thus the data for the vr+vr" and n'p correlations show a bunching effect characteristic 
for an attractive interaction while the data for the n^p correlation show an antibunching 
effect, characteristic for repulsion. After this test the authors compare their correction 
with the Gamow correction and find that the last one is much stronger. Herefrom they 
also conclude that the Gamow factor overestimates quite appreciably the Coulomb effect 



^^The measured tt+tt correlations were also used in two recent experimental papers |^ to 

estimate the Coulomb correction. Why such a procedure is questionable is explained below. 

^^This assumption has to be qualified among other things because final state strong interactions effects 

due to resonances can also influence these correlations. Furthermore there exists also a quantum statistical 

correlation for the tt+tt^ system (cf. sections 4.4, 4.7 and 4.8) which, however may be weak. 
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An even stronger conclusion is reached by Merlitz and Pelte |^ from the solution of 
the time dependent Schrodinger equation for two identical charged scalar bosons in terms 
of wave packets. These authors find that the "expected" Coulomb effect in the correlation 
function is obliterated by the dispersion of the localized states and is thus unobservable. 
This makes the interpretation of experimentally observed tt+tt^ correlations in terms of 
Coulomb effects even more doubtful. 

The theoretical studies of the Coulomb effect in BEC quoted so far are based on the 
solution of the Schrodinger equation and apply in fact only for the non-relativistic case. 
While one might argue that the relative motion of two mesons in BEC is for small q non- 
relativistic, this is not true for the single particle distributions (cf. below). Therefore in 
principle one should replace the solution of the Schrodinger equation in the Coulomb field 
used above by the corresponding solution of the Klein-Gordon equation. This apparently 



has not yet been done, with the exception of of a calculation by Barz^^ who investigated 
the influence of the Coulomb correction on the measured values of radii. He found an 
important change of these radii due to the Coulomb field only for momenta < 200 MeV. 

Finally, one must mention that the corrections of the wave function described above do not take into 
account the fact that the charge distribution of a meson is not point-like, but has a finite extension of 
the order of 1 fm. This means that in principle the Schrodinger (or the Klein-Gordon) equation has to 
be solved with a Coulomb potential modified by this finite size effect. Given the great sensitivity of BEC 
on small corrections in the wave function, this might be a worthwile enterprise for future research. As a 
matter of fact it is known from atomic physics (isotopic and isomeric shifts and hypcrfinc structure) that 
these finite size effects lead to observable consequences. 

Besides the wave function effect which influences the BEC due to directly produced 
particles or those originating from short lived resonances, one has to consider |^ the 
Coulomb overcorrection applied to pairs of which one particle is a daughter of a long lived 
state. This effect may bias the correction by up to 20%. 



^^Cf. however also ref.[Q where rescattering is added to the classical Coulomb effect and where 
somewhat different results are obtained. It is unclear whether the strong position-momentum correlations 
implicit in this rescattering model do not violate quantum mechanics. 
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Coulomb correction for higher order correlations The Coulomb corrections dis- 
cussed above were limited to single and two-body interactions. In present high energy 
heavy ion reactions we have already events with hundreds of particles and in the very 
near future at the relativistic heavy ion collider at Brookhaven (RHIC) this number will 
increase by an order of magnitude. Then many-body final state interactions may become 
important. Unfortunately the theory of many-body interactions even for such a "simple" 
potential as the Coulomb one is apparently still unmanageable. The long-range nature 
of the electromagetic interaction does not make this task simpler. Bowler |^ sketched 
a scenario for Coulomb screening based on the string model. In such a model particles 
are ordered in space-time, so that e.g. at least one 7r+ must be situated between the 
members of a tt^tt^ pair. While the net infiuence of this effect on tt+tt^ is expected to be 
small, for a tt^tt^ pair the situation is different, because instead of repulsion one obtains 
attraction. In the case of long range interactions such an effect may become important if 
the Til propagates together with the tt27t^ pair. This happens if Qu ~ Q13 ~ Q23- To 
take care of this effect Bowler suggests the replacement 

C(Q23) ^ C(Q23) < C-\Qu)C-\Qis) >fe3 (3.45) 

for the Coulomb tt^tt^ correction and 

C(Qi2) ^ C(Qi2) < C-\Q23)C{Qi3) >k3 (3.46) 

for the Coulomb tt+tt" correction. Here < ... >k-i symbolizes averaging with respect to 
the momentum of particle 3. While on the average 

< C-\Q23)C{Qr3) >k3 (3.47) 

is unity, at small Q the function C{Q) oscillates rapidly and therefore the factor 

< C-\Q,2)C-\Q,3) >k3 (3.48) 

is sensitive to the distribution of Qi2,Qi3 associated with the local source. According to 



^ this last factor for tt+tt pairs does not exceed 0.5% but for like sign pairs no estimate 



is provided. 



32 



Coulomb and resonance effect in single inclusive cross sections At a first look 
one might be tempted to believe that for single inclusive cross sections in heavy ion 
reactions the estimate of Coulomb effects is straightforward. Unfortunately this is not 
the case and so far there is no reliable theoretical estimate of this effect. This is so 
because the produced charged secondaries move not simply in the electromagnetic field of 
the colliding nuclei but interact at the same time with all the other secondaries. In view of 
this situation, recently an attempt has been made to put in evidence experimentally the 
Coulomb effect in the single inclusive cross section of pions in heavy ion reactions [|^. In 
this experiment an excess of negative pions over positive pions in Pb — Pb reactions at 158 



AGeV was observed which the authors of |^ atrributed to the Coulomb interaction of 
produced pions with the nuclear fireball. However this interpretation has been challenged 
in where, in a detailed hydrodynamical simulation it was shown that a similar excess 
in the 7r~/7r+ ratio is expected as a consequence of resonance (especially hyperon) decays. 
This qualification goes in the same direction as that mentioned above with respect to 
the exaggeration of the effects of Coulomb interaction in BEC. It also illustrates the 
complexity of the many body problem of heavy ion reactions even for weak interactions 
like the electromagnetic one, which in principle are well known. 



3.2 Strong final state interactions. 

This is a very complex problem because we are dealing with non-perturbative aspects of 
quantum chromodynamics. That there is no fully satisfactory solution to this problem 
can be seen from the very fact that we have at least three different approaches to it. 
As will become clear from the following, these different approaches must not be used 
simultaneously, as this would constitute double (or triple) counting. That is also why a 
fourth solution proposed here and which is of heuristic nature will appear in many cases 
more appealing and more efficient. 
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3.2.1 Final state interactions through resonances 

The majority of secondaries produced in high energy collisions are pions out of which a 
large fraction (between 40 and 80%) arise from resonances]^. Since the resonances have 
finite lifetimes and momenta, their decay products are created in general outside the 
production region of the "direct" pions (i.e., pions produced directly from the source) and 
that of the resonances. As a consequence, the two-particle correlation function of pions 
reflects not only the geometry of the (primary) source but also the momentum spectra 
and lifetimes of resonances [Q]. Kaons are much less affected by this circumstance p^ , 
however correlation experiments with kaons are much more difficult because of the low 
statistics. For a more detailed discussion of kaon BEC cf. section 5.1.3. 

The (known) resonances have been taken into account explicitely within the wave 
function formalism (cf. e.g. [|l3,|5^, [|18|), the string model (p2|, [Q) or within other 



variants of the Wigner function formalism (cf. e.g. [55|, |19|, |Q). The drawback 



of this explicit approach is that it is rather complicated, it is usually applicable only 
at small momenta differences q and it presupposes a detailed knowledge of resonance 
characteristics, including their weights, which, with few exceptions cannot be measured 
directly and have to be obtained from event generators^ or other models^. 

With these essential caveats in mind one finds that the distortion of the two-particle 
correlation function due to resonance decay leads to two obvious effects: (a) the effective 
radius of the source increases, i.e., the width of the correlation function decreases, and (b) 
due to the finite experimental resolution in the momentum difference the presence of very 
long-lived resonances leads to an apparent decrease of the intercept of the correlation 
function. Effect (b) is particularly important if one wants to draw conclusions from 
the intercept about a possible contribution of a coherent component in multiparticle 
production. 

In hydrodynamical studies of multiparticle production processes one considers reso- 
^^For experimental estimates of resonances cf. refs. 

^"^That event generators are not a reliable source of information for this purpose was demonstrated in 
the case of e+ — reactions in |5^, Q. 

^^In ref.pq] the weights were determined from thermodynamical considerations. 
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nances within the Wigner function approach (cf. section 5.1.3). 



3.2.2 Density matrix approach 



In ref. [^] one describes the effect of strong final state interactions by constructing a 
density matrix based on an effective Lagrangian of Landau- Ginzburg form as used in 
statistical physics and quantum optics (cf. also [§0|)- This is a more theoretical approach 
as it allows to study the effect of the strength of the interaction g on BEG, albeit in an 
effective Lagrangian description. 
One writes the density matrix 



P 



7-1 



vr 



(3.49) 



where F is the analogue of the Landau-Ginzburg free energy and the integrals are func- 
tional integrals over the field vr. This field is written as a superposition of coherent tt^ and 
chaotic fields iTch 

TT = TTciy) + T^ch{y) (3.50) 

The variable y refers in particular to rapidity. The total mean multiplicity < n > is 
related to the field tt by 



< n >=< 71^ >: 



(3.51) 



similar relations hold for the coherent and chaotic parts of vr. 

One assumes stationarity in y, i.e. the field correlator G{y, y') =< -K{y)-K{y') > depends 
only on the difference y — y' = Ay. One writes the Landau-Ginzburg form for F as 



Fin) 



dy 



an{y) + b 



d7i{y) 



dy 



(3.52) 



where a, b and g are constants. The strong interaction coupling is represented by g. The 
constants a, b can be expresed in term of < rich > and the "coherence length" which is 



defined through the correlator G (cf. eq.( |2.24| ). 

The main result of these rather involved calculations is the fact that while the inter- 
action does not play any significant role in the value of the intercept 6*2(0) it plays an 
important part in {C2{Ay 7^ 0). This situation is ilustrated in Figs. 3,4. Thus it is seen in 
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Fig. 3 that all C2 curves for various g coincide in the origin. The effect of the interaction 
in this approach is similar to that of (short lived) resonances, i.e. it leads to a decrease of 
the width of the correlation function. The sensitivity of C2 on g suggests that the shape 
of the correlation function can in principle be used for the experimental determination 
of g. One finds furthermore that there is no (^-dependence for purely chaotic or purely 
coherent sources. This observation suggests that for a strongly coherent or chaotic field 
the final state interaction does not manage to disturb the correlation. 

Note that the curves in Fig. 4 intersect at some Ay which depends on the chaotic- 
ity. This is characteristic for correlations treated by quantum statistics, when one has a 
superposition of coherent and chaotic fields, and is a manifestation of the appearance of 
two (different) functions in the correlation function . 

We recall (cf. section 2.2) that for the particular case of a Lorentzian spectrum one 
gets, in the absence of final state interactions, 

C^iy, y') = l + 2p(l - p)e-^y-y'^/^ + p2g-2||y-y'|/c ^3 53^ 

instead of the empirical relation 

C2{y, y') = 1 + Aexp(-|y - y'|/0- (3-54) 



As shown in Fig. 5 the parametrization ( |3.54| ) leads to parallel curves for various values 
of A, while the more correct parametrization of ref. |^ leads to intersecting curves. This 



is due to the fact that the relation for C2 derived in contains as a particular case 
(for g = 0) eq. (|3.53|) and retains the essential feature of eq. ( |3.53| ), which consists in the 
superposition of two exponentials. 

3.2.3 Phase shifts 

For charged pions the strong final state interactions can be described also by phase shifts. 
It is known that for an isospin 1 = 2 state (this the isopsin of a system of two identically 
charged pions) the corresponding strong interaction is repulsive. However it was suggested 
)l[| that the range of strong interactions is smaller than the size of the hadronic source 
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and therefore the correlation should be essentially unaffected by this effect^. Even for 
particle reactions like hadron-hadron or e^ — e~ the size of the source is of the order of 1/m 
while the range of interaction is only 0.2/m. On the other hand the effective size quoted 
above arises because of the joint contribution to BEG of direct pairs and resonances. So 
it is interesting to analyze these two contributions separately. 

Two identically charged pions are practically always produced together with a third 
pion of opposite charge . Then according to Bowler one has also to consider the / = 
attractive interaction between the oppositely charged pions and this compensates largely 
the 1 = 2 state interaction. Thus it appears that also in particle reactions only resonances 
play an important role in final state strong interactions. 

The considerations about final state interactions made in this subsection treat sepa- 
rately Coulomb and strong interactions. This is permitted as long as we deal with small 
effects or when the ranges of the two types of interactions do not overlap. For very small 
distances this is not anymore the case. Furthermore the Gamow and the phase shift 
corrections are based on the wave function formalism which ignores the possibility of cre- 
ating particles. However when entering the non-classical region the well known difficulties 
of the wave function formalism become visible (Klein paradox). To consider this effect, 
in ref.(|6^ the joint contribution of the strong interaction potential and the Coulomb 
potential are analyzed in a version of the Bethe-Salpeter equation for spinless particles. 
It is found that as expected also from the considerations presented above that strong 
interaction diminishes appreciably the Gamow correction. 

3.2.4 Effective currents 

As will be explained below the most satisfactory approach to BEG is at present the 

classical current approach, based on quantum field theory. There three types of source 

chracteristics appear: the chaoticity, the correlation lengths/times and the space-time 

dimensions. It is obvious that these quantities contain already information about the 
^^The separation between resonances and phase shifts is of course not rigorous because phase shifts 
reflect also the effect of resonances; however as long as phase shifts constitute a small effect, this should 
not matter. 
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nature of the interaction and therefore it is quite natural to consider them as effective 
parameters which describe all effects of strong final state interactions. 

This approach to strong final state interactions is probably the most recommendable 
at the present stage for BEC in reactions induced by hadrons or leptons because: 1. it 
is simpler; 2. it avoids double counting; 3. it avoids the use of poorly known resonance 
characteristics; 4. it avoids the use of the ill defined concept of final state interactions for 
strong interactions. For heavy ion reactions, when hydrodynamical methods are used the 
explicit consideration of resonances can be practiced up to a certain point without major 
difficulties and then the strong final state interactions can be taken into account through 
these resonances (cf. section 5.1.3). 

To conclude this discussion of final state interactions in BEC, it is interesting to note 
that in boson condensates the final state interactions might be different than in normal 
hadronic sources. In a condensate the Bose field becomes long range in configuration 
space. This can be understood as a consequence of the fact that in a condensate the 
effective mass of the field carrier vanishes. Indeed a calculation based on the chiral 
sigma model shows that the effective range of the pion field can increase several times due 
to this effect. 

4 Currents 

4.1 Classical versus quantum currents 

In section 2.2 we were concerned mainly with the properties of fields and did not ask the 
question where these fields come from. In the present chapter we shall put and try to 
answer this question. 

We start by recalling the definition of correlation functions within quantum field the- 
ory. 

Let aj(k) and ai(k) be the creation and annihilation operator of a particle of mo- 
mentum k, where the index i labels internal degrees of freedom such as spin, isospin, 
strangeness, etc. The n-particle inclusive distribution is 
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1 Jn-.il. ..in " / 

= (2vr)3«n2E,Tr(p^a^k0...aljkOa.„(kO...a.,(k0) (4.1) 

where 

is the invariant volume element in momentum space. With the notation 

Q--(ki,...,k.) = -^^-r- (4-3) 

0" auji...aujn 

the general n-particle correlation function is defined as 

The particles which the operators and a are associated with are the quanta of the 
field (which we will denote in general by 0); these particles as well as the density matrix 
p refer to the final state where measurements take place. On the other hand we usually 
know (or guess) the density matrix only in the initial state. Therefore we will have to 
transform the above expression so that eventually the density matrix in the final state pf 
is replaced by the density matrix in the initial state pi, while the fields will continue to 
refer to the final state. To emphasize this we wrote in eg. (|4.1| ) pf. 

Pf = SpiS^ (4.5) 

so that we have 

Gi{k) = {27T){2Ei)Tr{p,S^a\k)a{k)S} (4.6) 

G2{ku k^) = {2'Kf{2E^){2E2)Tr{piS^a\k2)a\k^)a{k^)a{k2)S} (4.7) 

Thus, if the initial conditions i.e. pi are given, in principle the knowledge of the S 
matrix suffices to calculate the physical quantities of interest. In one case the 5* matrix 
can even be derived without approximations. This happens when the currents are classical 
and we shall discuss this case in some detail in the present and the following section. 

Before doing this, we shall consider briefiy the more general case when the currents 
are not necessarily classical. 



39 



The 5* matrix is given by the relation 

S = Texpii f d^xU^'^x) \ (4. 



where the interaction lagrangian L^nt is a functional of the fields 0. T is the chronological 
time ordering operator; we shall use below also the antichronological time operator 7". 
Consider for simplicity a scalar field produced by a current J. Then 

= J(x)0(x). (4.9) 

Equations (|4.8| ) and (^4.9|) allow us now to calculate the correlations we are interested 
in in terms of the currents after eliminating the fields. One obtains thus 

Pi(fc) = rr{p,4(k)JH(k)} (4.10) 

P2(fei , fc2) = Tr {p,f [Jl,(ki) Ji,(k2)] T [J^(ki) J^(k2)]} (4.11) 

where the label H stands for the Heisenberg representation. Now the cross sections depend 
only on the currents and the density matrix in the initial state. The appearance of the 
time ordering operators T and T in eqs.( 4.ld|) ,( [4.11| ) is a reminder of the fact that the 



current J is here an operator. 
4.2 Classical currents 

Besides the fact that in this case an exact, analytical solution of the field equations is 
available and that the limits of this approximation are quite clear, the classical current has 
the important avantage that in it the space-time characteristics of the source are clearly 
exhibited and thus contact with approaches like the Wigner approach and hydrodynamics 
are made possible. 

The assumption that the currents are classical implies that J is a c number and 
then the order in eq.( [4.11| ) does not matter. This approximation can be used in particle 



physics when the currents are produced by heavy particles (e.g. nucleons) and/or when the 
momentum transfer q is small compared with the momentum K of the emitting particles. 
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0. In section 4.7 a new criterion for the applicability of the classical current assumption 
in terms of particle-antiparticle correlations will be presented. 

The classical current formalism was introduced to the field of Bose-Einstein correla- 
tions in B3, 153]. In this approach, particle sources are treated as external classical 



currents J{x), the fluctuations of which are described by a probability distribution P{J}- 
^From many points of view like e.g. the understanding the space-time properties of the 
sources or the isotopic spin dependence of BEC this approach is superior to any other 
approach. This has become clear only in the last years [Q, when a systematic inves- 



tigation of the independent physical quantities which enter the dynamics of correlation 
functions has been made (cf. below). 

The classical current formalism in momentum space is mathematically identical with 
the coherent state formalism used in quantum statistics and in particular in quantum 
optics (cf. section 2.2), the classical currents in fc-space J{k) being proportional to the 
eigenvalues of the coherent states \a >. This explains the importance of the coherent 
state formalism for applications in particle and nuclear physics. 

The density matrix is 

p = y PJ V{J} I J >< J| (4.12) 

where the symbol VJ denotes an integration over the space of functions J{x), and the 
statistical weight V{J} is normalized to unity, 

J VJV{J} = 1. (4.13) 

(The reader will recognize in eq.( [4.12| ) the P-representation introduced in section 2.2.) 



Expectation values of field operators can then be expressed as averages over the cor- 
responding functional of the currents, e.g., 

Tr (pat (ki) ... at(k„)a(k„+i) ... a(k„+„)) (4.14) 
fJi J(27r)3 2Ej ,=ra+i J{2^Y2E, 



^•^For multiparticle production processes this implies also a constraint on the multiplicity and/or the 
momenta k of the produced particles 
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In the following we shall discuss the special case when the fluctuations of the currents 
J{x) are described by a Gaussian distribution V{J}. The reasons for this choice were 
given in section 2.2. 

As in quantum optics we write the current J{x) as the sum of a chaotic and a coherent 
component, 

with 

Jcoherentix) = < J{x) > (4.16) 
Jchaotic{x) = J{x)- < J{x) > (4.17) 

By deflnition, < Jchaotic{x) >= 0. The case I{x) ^ corresponds to single particle 
coherence. 

The Gaussian current distribution is completely determined by specifying its flrst two 
moments: the flrst moment coincides, because of eg. ( [1. 17] ), with the coherent component, 

I{x) = < J{x) > (4.18) 

and the second moment is given by the 2-current correlator, 

D{x,x')=< J{x) J(x') > - < J{x) > < J{x') > 

= < Jchaotic{x) Jchaotic{x') > (4.19) 



4.3 Primordial correlator, correlation length and space-time dis- 
tribution of the source 

We come now to a more recent development of the current formalism which has 



shed new light on both fundamental and applicative aspects of BEG. 

There are two fundamental and in principle independent aspects of physics which 
come together in the phenomenon of BEG in particle and nuclear physics. One refers 
to the geometry of the source and goes back to the original Hanbury-Brown and Twiss 
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interefence experiment in astronomy. The "geometry" is characterized by the size of the 
source e.g. the longitudinal and transverse radius R\\ and R± respectively and the lifetime 
of the source Rq. 

The second aspect is related to the dynamics of the source and is expressed through 
correlation lengths. In the following we will use two correlation lengths L± and a 
correlation time Lq As a consequence of the finite space-time size of sources in particle 
physics one cannot in general separate the geometry from the dynamics in the second 
(and higher) order correlation function. This separation is possible (cf. below) only by 
using simultaneously also the single inclusive cross section. This remarkable fact |]53| is 
not yet realized in most of the theoretical and experimental papers on this subject]^. One 
of the main objects of the present review is to clarify this situation. 

Assuming Gaussian currents, one may take the point of view that the purpose of 
measuring n-particle distributions is to obtain information about the space-time form of 
the coherent component, I{x), and of the correlator of the chaotic components of the 
current, D{x,x'). In practice, because of limited statistics it is necessary to consider 
simple parametrizations of these quantities, and to use the information extracted from 
experimental data to determine the free parameters (such as radii, correlation lengths 
etc.). 

The approach considered in differs in some fundamental aspects from those applica- 
tions of the density matrix approach in particle physics which are performed in momentum 
(rapidity) space and are limited usually to one dimension (cf. however [[7^] where rapidity 
and transverse momentum were considered). The approach of ref. is a space-time ap- 
proach in which the parameters refer to the space-time characteristics of the source. This 
approach has important heuristic advantages compared with the momentum (rapidity) 

space approach as will be explained below. 

^''We refer here to short range correlations. Cf. section 6.1.3 for a distinction between short range and 
long range correlations. 

^^For systems in local equilibrium Makhlin and Sinyukov j6^] introduced a length scale (called in |^ 
"length of homogeneity") which characterizes the hydrodynamical expansion of the source and can be 
different from the size of the system. Further references on this topic can be found e.g. in [|68|, 
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Among other things, in the quantum statistical (QS) space-time approach the pa- 
rameters of the source as defined above can be considered as effective parameters which 
contain aheady the entire information which one is interested in and which one could 
obtain from experiment, and thus distinguishing between directly produced particles and 
resonance decay products could amount to double counting. The apparent proliferation 
of parameters brought about by the QS approach is compensated by this heuristic and 
practical simplification. Furthermore a new and essential feature of the approach of ref. 
as compared with previous applications of the current formalism [|65| , |39| which assume 
L = 0, is the finite correlation length (time) L. This fact has important theoretical and 
practical consequences. It leads among other things to an effective correlation between 
momenta and coordinates, so that e.g. the second order correlation functions does not 
depend only on the difference of momenta q = ki — k2 but also on the sum ki + k2. 
This non-stationarity property, which is observed in experiment, is usually associated 
with expanding sources and treated within the Wigner function formalism. However from 
the considerations presented above it follows that expansion is in general not a neces- 
sary condition for non-stationarity in q. It will be shown how expanding sources can be 
treated without the Wigner formalism, which restricts unnecessarily the applicability of 
the results to small q values (cf. section 4.8). 

The distinction between correlation lengths and radii is possible only in the current 
formalism; the Wigner formalism provides just a length of homogeneity. 

The results which follow from the space-time approach include: 

(i) The existence of at least 10 independent parameters that enter into the correlation 
function; 

(ii new insights into the problem of partial coherence; 

(iii) Isospin effects: Tr'^vr^ correlations are different from tt^tt^ correlations; there ex- 
ists a quantum statistical (anti) correlation between particles and antiparticles (tt+tt^ in 
this case). These effects are associated with the presence of squeezed states in the den- 
sity matrix, which in itself is a surprising and unexpected feature in conventional strong 
interaction phenomenology. 

It turns out that soft pions play an essential role in the experimental investigation of 
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BEC, both with respect to the effect of particle-antiparticle correlations predicted here 
as well as in the investigation of the coherence of the source. Depending on the relative 
magnitude of the parameters of the coherent and chaotic component, soft particles can 
either enhance or suppress the coherence effect. 

Consider first the case of an infinitely extended source. The correlation of currents at 
two space-time points x and y is described by a primordial correlator, 

<J{x) J{y)>o= C{x-y) (4.20) 

Note that C depends only on the difference x — y. The correlator C{x — y) contains some 
characteristic length (time) scales L, the so-called correlation lengths (times) In the 
current formalism used in L = 0. 

C{x — y) is a real and even function of its argument. In the rest frame of the source, 
it is usually parametrized by an exponential. 



C{x — y) = Cq exp 

or by a Gaussian, 

C{x — y) = Co exp 



Fo - Vol 



^0 



(4.21) 



(xo - yo^ (x-y 



i2 



2Ll 2L2 



(4.22) 



However, it should be clear that in principle any well behaved decreasing function of 
|x — y| is a priori acceptable, and in practice it is usually up to the experimenter to decide 
which particular form is more appropriate. The ansatze (|4.21| , ^^2^ ) need to be modified 
for the case of an expanding source (cf. sections 4.8 and 4.9) where each source element 
is characterized not only by a correlation length but also by a four velocity m^. As a 
matter of fact, the form of the function C is irrelevant as long as one is interested in the 
general statements of the theoretical quantum statistical (current) formalism. In practical 
applications, of course, in order to obtain concrete information about the source and the 
medium (i.e., about L) the form of C has to be specified. In principle, a full dynamical 
theory is expected to determine the functional form of the correlation function, however 
These correspond to the "coherence lengths" used in the quantum optical literature. 



45 



at present this "fundamentalist" approach is not yet apphcable 0. One uses instead a 
phenomenological approach hke that reflected in eqs. ([4.21| , f4.22|) . 

Effects of the geometry of the source are taken into account by introducing the space- 
time distributions of the chaotic and of the coherent component, fch{x) and fc{x), respec- 
tively. The expectation values of the currents, I{x) and D{x, x'), take nonzero values only 
in space-time regions where fc and fch are nonzero. Thus, one may write 

I{x)=Ux) (4.23) 

D{x, x')=U{x) C{x - x') U{x') (4.24) 
We turn now to an important new aspect of the current approach 

4.4 Production of an isospin multiplet 

Following we generalize the previous results to the case of an isospin multiplet and 
derive explicit expressions for the single inclusive distributions and correlation functions 
of particles that form an isotriplet (such as the 7r"'",7r~ and vr'^-mesons) . For the sake of 
definiteness, we will refer to pions in the discussion below, but it should be understood 
that the formalism is applicable to an arbitrary iso multiplet. 

Consider the production of charged and neutral pions, 7r+,7r^ and n^, by random 
currents. The initial interaction Lagrangian is written 

Cint=J+ix)7i-{x) + J_(x)7r+(x) + Jo(x)7r°(x). (4.25) 

The current distribution is completely characterized by its first two moments. They read 
for the case of an isotriplet 

Ii{x)=< Ji{x) > = +,-,0 (4.26) 

Dii'{x,x')=< Ji{x) Ji'{x') > — < Ji{x) > < Ji'{x') > 
^^Indeed one may hope that for strongly interacting systems lattice QCD may provide in future the 
corrrelation length L. 
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Invariance of the chaotic 2-current correlator under rotation in isospin space imphes 



Doq{x,x')=D^_{x,x') = D{x,x') 



(4.27) 



D++{x,x')=D__{x,x') = D+o{x,x') = Do_{x,x') = 0. 



The corresponding current distribution is 



P{Ji} = ^ exp[-A{J^}] 



where in coordinate representation 



and 



with 



^f ^ J J J VJ+VJJDJq exp[-^{Ji}] 



A{Ji}^ 1 1 d'^x d^y[{J+{x) - I+{x))M{x,y){J4y) - I^y)) 



+ -iUx) - Io(x))M(x, y)(Jo(y) - Io(y)) 



M{x,x') = D-\x,x') 



(4.28) 



(4.29) 



(4.30) 



(4.31) 



For the general case of a partially coherent source, the single inclusive distributions of 
pions of charge i (i = +, — , 0) can be expressed as the sum of a chaotic component and a 
coherent component, 



1 da' 1 da' 



with 



and 



a du a du 

1 da 



+ 



1 da' 



chaotic 



a du 



coherent 



a du) 
1 da' 



chaotic 



a du 



= D{k,k) 



\m\' 



coherent 

In general, the chaoticity parameter p will be momentum dependent. 



p{k) 



D{k,k) 



D{k,k) + \I{k)\'^ 



(4.32) 



(4.33) 



(4.34) 



(4.35) 
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To write down the correlation functions in a concise form, one introduces the normal- 
ized current correlators, 

DiKJ^ J D{K,-h) 

[D{kr,kr)-D{ks,ksp [D{kr,kr)-D{ks,hp 

where the indices r, s label the particles. Note that d is the same function as d but of the 
change of sign of one of its variables. One may express the correlation functions in terms 
of the magnitudes and the phases, 

Trs=T(^kr, kg) = /c^) I 

Trs=T(^kr, kg) = \d(^kr,ks)\ 

C,=(t)'\K,h) = Aigd{kr,ks) (4.37) 
4>';':=^'\K,h) = kigd{K,ks) 

and the phase of the coherent component, 

ct>l = <P%kr) = Arg I{kr) . (4.38) 
The same notation will be used for the chaoticity parameter, 

Pr = p{kr) (4.39) 

The two-particle correlation function is 



C++{k,,k2) = 1+ 2^pi(l-pi)-P2(l-P2)ri2 cos(0^^-0^ + 0^) + p^p,T^, (4.40) 

In higher order correlation functions up to and including order 5 are given. 
In the absence of single particle coherence the two-particle correlation functions for 
different pairs of tt^, vr^, tt" mesons read 

C2++(ki,k2) = l+Ml2r, 



|2 



Cr(ki,k2)=l + |rfi2| 
C2+°(ki,k2)=l, (4.41) 

C°°(ki,k2) = l+Ml2|'+M~i2|' 
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These results were surprising in that they disagreed with some of the preconcei- 
ved notions on Bose-Einstein correlations. For instance, it was commonly assumed that 
without taking into account final state interactions and in the absence of coherence, the 
maximum of the two-particle correlation of identical pions is 2 (for ki = k2). It was also 
assumed that there are no correlation effects among different kinds of pions because these 
particles are not identical. (This last assumption is even sometimes used in normalizing 
the experimental data on with respect to C2 .) 

The results ( |4.41| ) show that these assertions are not necessarily true. In particular, 
looking at the two pion correlations one can see that in addition to the familiar correla- 
tions of identical particles (the terms |c?i2p) there are particle-antiparticle - in this case, 
TT'^TT~ - correlations (the terms |(ii2p). The 7r° has both terms, as it is identical with 
its antiparticle. Essentially this last fact is the explanation for the appearance of the 
"surprising" effects. Obviously it is a specific quantum field effect. It will be shown in 
section 4.8, that for soft pions and for small lifetimes of the source the terms |(ii2p can 
in principle become comparable with the conventional terms |c?i2p. This implies that the 
distribution 6*2" (ki, k2) of two neutral pions can be as large as 3, and the maximum value 
of C2 is 2 (instead of 1). The corresponding limit of C^^'^ is 15 (instead of 6), and that 
of 6*3"^" is 6 (instead of 2). However, it should be noted that these are merely upper 
limits, which for massive particles are not reached except for sources of infinitesimally 
small lifetimes. For soft photons however the situation is different (cf. below). 

The "new" terms, proportional to < a^{ki)ae^{k2) > are due to the non-stationarity 
(in k space) of the source. While in quantum optics time-stationarity is the rule, in 
particle physics this is not the case because of the finite lifetime and finite radius of the 
sources. The existence of a non-vanishing expectation value of the products a{ki)a{k2) 
is what one would expect (cf. section 2.2) from two-particle coherence (squeezing), just 
as < a{k) >y^ follows from ordinary (one-particle) coherence (note that the latter has 
not been assumed here). The fact that squeezing which, as mentioned above, is quite 
an exceptional situation in optics, discovered only recently, is a natural consequence of 
the formalism for present particle physics, is possibly one of the most startling results 
obtained recently in BEG. A characteristic feature of squeezed states is the fact that for 
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static sources they lead anti-correlations. Finally, it should be pointed out again that the 
above results - in particular, the fact that in general differs from C2 - are consistent 
with isospin symmetry. 

In closing this section one should note that the existence of particle-antiparticle corre- 
lations is not restricted to pions but applies also to other systems, like e.g. neutral kaons. 
In principle thus there exist also KqKq quantum statistical correlations. However since 
Kq particles cannot be observed except in linear combinations with Kq in the form of 
Kg and Ki the QS particle-antiparticle correlation effect has to be disentangled from the 
KgKg or KiKi Bose-Einstein correlation [Kg and Ki are of course bosons and thus subject 
to BEC) which exist also in the "old fashioned" wave function formalism which ignores 
the intrinsic "new" KqKq correlation. As a matter of fact KgKs correlations have been 
observed experimentally, however no attempt has been made so far to extract from them 
the "surprising" effects (For more recent experiments cf. |7^) and for a theoretical 



analysis of the "old" effects and their possible application in CP violation phenomena cf. 

!)• 



^^Because of their larger mass as compared with that of pions, these effects may be even more quenched 

than in the case of pions, except for sources of very short hfetime (cf. section 4.8). 

^^The analysis in the preceding subsections referred to the production of an isotriplet assuming just 



symmetry between the isospin components (cf. eq.(4.25)) of the current. In principle, for strong inter 



actions the conservation of isospin / has also to be considered. While the chaotic part is not affected by 
this condition, the coherent component is influenced by conservation of isospin]?^. In particular this can 
lead to an additive positive term in the correlation function and thus to an increase of the bounds of BEC 
for pions. It remains to be seen whether this effect can be distinguished from the effect of long range 
correlations (cf. section 6.1.3). Moreover in hadronic reactions and in particular those involving nuclei 
such an effect would be strongly suppressed because of the following circumstance. The initial state has 
to be averaged over all components of isospin / which is a first "diluting" factor. Furthermore the effect 
is appreciable only for low total isospin. This total isospin has to be shared by the chaotic component Idi 
and the coherent component Ic'. I = Ic + Ich- The first one arises mostly from resonances with different 
isospin values, so that even if the total isospin takes its minimum value (/ = 0), Ich and therefore Ic can 
take larger values. 
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4.4.1 An illustrative model of uncorrelated point-like random sources 

To clarify the origin of different terms in the functions C|''(/ci, /C2), let us consider, a source 
consisting of N point-like random sources 

AT 

M^) = E jai^i) - ti) -^)^ a = +, -, (4.42) 

i=l 

and assume that the currents ja{xi) at different points Xi are mutually independent 
and have the same statistical properties, i.e. 



< Ja(^i) jh{xj) >= 5ij < fjb > . (4.43) 
We also assume in this section 



< ja >= 0, (4.44) 

i.e. ignore a possible coherent component < ja > to make the presentation more trans- 
parent. 

Now the one-particle distribution is 

<j:{k)Uk)>=N.<j:ja> (4.45) 

and the two-particle distribution takes the form 

<j:{ki) Jaih) j;{k2) Jb{k2) > = 



N 



i=l 
N 

+E < j*a{Xi)ja{Xi) > ■ < j;{Xj)jbiXj) > 

N 

N 



(4.46) 
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Let us consider separately the four different terms on the right-hand side of eq. ( [4.46| ) . 
The first term corresponds to two particles being emitted from a single point (Fig. 6a); 
it is proportional to the number of emitting points. 

The second term describes an independent emission of two particles from different 
points (Fig. 6b). 

The third term, being non-zero for a = b, describes an interference effect of direct and 
exchange diagrams, characteristic for identical particles, emitted from different points 
(Fig. 6c). This is the usual BE-correlation effect. 

The fourth term describes an interference of two-particle emissions from different 
points (Fig. 6d) ("two-particle sources"). It appears to be non-zero for real currents 
with a = b (7r*^7r^) and for complex currents with J* = Jf, (vr+vr"), that is for particle- 
antiparticle associative emission. In this simple model, it represents the "surprising" 
effects discussed above. 

This diagramatic illustration of the "surprising" BE-correlations is due to Bowler , 



who found that the "surprising" effects derived for the first time in ref. can be 

understood in terms of these more qualitative considerations. Bowler derived the "new" 
effects from the string model. 

4.5 Photon interferometry. Upper bounds of BEC. 

The advantage of photon BEC resides in the fact that photons are not influenced by final 
state interactions. Photons present also an interesting subject of theoretical research from 
the general BEC point of view, since they are spin-one bosons while pions and kaons used 
in hadronic BEC are scalar particles. We shall see below that this supplementary degree 
of freedom has specific implications in BEC. Last but not least photon correlations are 
for various reasons to be discussed below of particular interest for the search of quark 
matter. We present some of the results of |7^, which contain as a special case those of 



| 77| and where these topics are discussed. 

Consider a heavy ion reaction where photons are produced through bremsstrahlung 
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from protons in independent proton-neutron collisions^. The corresponding elementary 
dipole currents are 



1 P 

where p = pj — is the difference between the initial and the final momentum of the 
proton, eA is the vector of linear polarization and k the photon momentum; e and m are 
the charge and mass of the proton respectively. The total current is written 

N 

j\k) = J2e^'^^"j::ik). (4.48) 

n=l 

For simplicity we will discuss in the following only the case of pure chaotic currents 
< J^{k) >= and refer for coherence effects to the original literature |]76l, [^. In analogy 
to the considerations of the previous subsection the index n labels the independent nucleon 
collisions which take place at different space-time points x„. These points are assumed 
to be randomly distributed in the space-time volume of the source with a distribution 
function f{x) for each elementary collision. The current correlator is proportional to 
products of the form 

N 



< J''ik,)j'^{-k,) >= e\^{k^) < pl^pi >^ ei^ik,) 



(4.49) 



For central collisions due to the axial symmetry around the beam direction one has 
for the momenta the tensor decomposition 

< P>i >= l^J'' + SJ'P, (4.50) 
where / is the unit vector in the beam direction and (j„, 5„ are real positive constants. 



In IT3] an isotropic distribution of the momenta was assumed. This corresponds to the 



particular case (5„ = 0. The generalization to the form ( [4. 501) is due to |76|]. The summation 



over polarization indexes is performed by using the relations 



<{e\pi){e^.p,,)>=^{e\e^)6u' (4.51) 
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Photon emission from proton-proton collisions is suppressed because it is of quadrupole form. 



53 



and 



j2e\{k)ei{k)=6^^ -n^n\ (4.52) 

A=l 

where n = k/|k|. 

We write below the results for the second order correlation function 

C,ih,k,)= 'fl^^'l, (4.53) 

for two extreme cases: 

(1) Uncorrelated elementary currents (isotropy) (cr ^ 5) 



C2{k,, k^; a j^O, S = 0) = ! + ][! + (ni.na)^] \\f{k^ - ^2)^ + \f{k^ + fca)!'! , (4.54) 



4 

leading to an intercept 

C2{k,k) = ^- + ^\f{2k)\' (4.55) 

limited by the values (3/2,2). 
(2) Strong anisotropy (cr ^ 6) 

C^ik,, k^; a = 0, S^0) = 1 + \f{k, - k^)]' + \f{k, + k^)]'' (4.56) 

with an intercept 

C2{k,k)=2 + \~f{2k)\'' (4.57) 

limited this time by the values (2,3). Note that due to the form of the photon-current in- 
teraction ( 4.47|) in this (strong anisotropy) case the photons emerge practically completely 



polarized so that the summation over polarizations does not affect the correlation. 

These results are remarkable among other things because they illustrate the specific 
effects of photon spin on BEC. Thus while for (pseudo-) scalar pions the intercept is a 
constant (2 for charged pions and 3 for neutral ones) even for unpolarized photons the 
intercept is a function of k. 

For a graphical illustration and explanation of this fact cf. Fig. 7 from ref. [^). It 
is seen that to perform the summation over polarization implied by eq.( [4.51| ) only one 
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direction of the linear polarization can be chosen equal for both photons, while the other 
polarization direction differs by the angle 9 between the two momenta ki, k2. 

One thus finds that, while for a system of charged pions (i.e. a mixture of 50% positive 
and 50% negative) the maximum value of this intercept MaxC2{k, k) is 1.5, for photons 
MaxC2{k, k) exceeds this value and this excess reflects the space-time properties (repre- 
sented by f{k)) , the degree of (an)isotropy of the source represented by the quantities a 
and 6, and the supplementary degree of freedom represented by the photon spin. 

As a consequence of the fact that / is a decreasing function of its argument, in 
eqs.( [4.54D ,( [4.56D the terms with f{ki + k2) are in general smaller than the terms with 
f{ki — ^2), except for small momenta k. 

The fact that the differences between charged pions and photons are enhanced for soft 
photons reminds us of a similar effect found with neutral pions (cf. section 4.4). Neutral 
pions are in general more bunched than identically charged ones and this difference is more 
pronounced for soft pions. This similarity is not accidental, because photons as well as 
TT^ particles are neutral and this circumstance has quantum field theoretical implications 
which will be mentioned also below. 

We see thus that photon BEG can provide information both about the space-time 
form of the source represented by / and the dynamics which are represented by 5 Q Q. 

The results on photon correlations presented above refer to the case that the sources 



are "static" i.e. not expanding. Expanding sources were considered in [32| within a 
covariant formalism. 

Some of the results above, in particular eqs. (|4.54| , ^755| ), which had been initially 



derived by Neuhauser were challenged by Slotta and Heinz |^3|. Among other things, 
these authors claim that for photon correlations due to a chaotic source "the only change 
relative to 2-pion interferometry is a statistical factor | for the overall strength of the 



correlation which results from the experimental averaging over the photon spin" . In 
an intercept | is derived which is in contradiction with the results presented above and in 
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Cf. [k8|,|7MI where photon correlation experiments in low energy (100 MeV/nucleon) heavy ion 



reactions are reported. For a theoretical discussion of these experiments cf.|80 . 

■^^The relation between photon interferometry and the formation length of photons is discussed in pi[ 
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particular with eq.( [4.55D where besides the factor | there appears also the k dependent 

function ||/(2/c)p. Similar statements can be found in previous papers IQ, 

where more detailed applications concerning heavy ion reactions based on this assertion 



of are presented. 

Some of the papers quoted above were criticized immediately after their publication 
in ||8^ and the paper was written with the intention to "settle" this "con- 

troversy" . 

It should be pointed out here that the reason for the difference between the results 



of ||7^,0] on the one hand and those of ref. |^ on the other is mainly due to the fact 
that in ||83[ a formalism was used which is less general than that used in []7^ and 
and which is inadequate for the present problem. This implies among other things that 



unpolarized photons cannot be treated in the way proposed in p3| and that the results 
of and are correct, while the results of are not. 

In ||8^ the following formula for the second order correlation function is used: 

^^,(q,K)r^(-q,K) 



C(ki,k2) = l + 



(4.5^ 



Here g is the Fourier transform of a source function, q = ki — k2 and K = i(ki + k2). 
This formula is a particular case of a more general formula for the second order correlation 
function derived by Shuryak using a model of uncorrelated sources, when emission of 
particles from the same space-time point is negligible (cf. section 4.4.1). 

Since this equation is sometimes used in the recent literature without giving the reader 
the possibility of evaluating the approximations used in its derivation, we will sketch this 
derivation in the following. 



In ||6^ one starts with the current correlator 

< J*{xi)Jj{x2) >= 6ijJi{x,Ax), 
where Ji{x) is the current emitted by point x and 

{Xi+X2)/2, Ax = Xi-X2. 



X 



(4.59) 



(4.60) 



Eq. ( [4.59| ) assumes that the individual currents {i 7^ j) are uncorrelated. With the 
notation Ii{q,K) for the Fourier transform of Ii{x,Ax) and I{q,K) = J^i^iQ^K) the 
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inclusive single particle distribution reads 



W{k) 



J2mk) = iio,k) 



(4.61) 



and the two particle distribution is given by 

W{kik2) = 



ikixi+ik2X2 _j_ ^ikiX2+ik2X 



^)Ji{xi)Jj{x2)dXidX2 



(4.62) 



or finally 



W{k^, fc2)=/(0, fci)/(0, k2) + \i{q, K)\' 



+ E[< J:{ki)J:{k2)Mki)Mk2) > 



-i,{0,k,)i,{0,k2)-\Uq,K)\' 



(4.63) 



The first term on the rhs of eg. ( [4.63| ) is the product of one-particle distributions and 
the second term is the conventional interference term, corresponding to figure 6c. By 
going over to the Wigner source function (cf. also section 4.9) 



d'ye'^^y <j;{x + -y)Mx--y)> 



(4.64) 



these first two terms result in eq.( [4.58| ) of ||83|. However as is clear from eq.( |4.63D there 
exists also a third term, neglected in eq. ( |4.58| ) and which corresponds to the simultaneous 
emission of two particles from a single point (xj) as indicated in figure 6d. While for 
massive particles this term is in general suppressed, this is not true for massless particles 
and in particular for soft photons. In ||7^ and this additional term had not been 
neglected as it was done in and therefore it is not surprising that ref. could not 



recover the results of refs. [^] and . The neglect of the term corresponding to emission 
of two particles from the same space-time point is not permitted in the present case. As 
mentioned in section 4.4.1, in a model of uncorrelated point-like random sources like the 
present one, emission of particles from the same space-time point corresponds in a first 
approximation to particle-antiparticle correlations and this type of effect leads also to 
the difference between BEC for identical charged pions and the BEC for neutral pions. 
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This is so because neutral particles coincide with the corresponding antiparticles. (As 
a consequence of this circumstance e.g. while for charged pions the maximum of the 
intercept is 2, for neutral pions it is 3 (cf. section 4.4)). Photons being neutral particles, 
similar effects like those observed for vr^-s are expected and indeed found (cf. above). 

This misapplication of the current formalism invalidates completely the conclusions of 
ref.[^ and confirms and strengthens the criticism expressed in [|^, |^9|, of the papers 



8^-|8^. The fact that for unpolarized photons MaxC2{k, k) is 2 and not 1.5 as stated in 
can be understood by realizing that a system of unpolarized photons consists on the 



average of 50% photons with the same helicities and 50% photons with opposite helicities. 
The first ones contribute to the maximum intercept (of the unpolarized system) with a 
factor of 3 and the last ones with a factor of 1 (coresponding to unidentical particles). 

For the sake of clarification it must be mentioned that ref. |^ contains also other 
incorrect statements. Thus the claim in |Q that the approach by Neuhauser "does 
not correctly take into account the constraints from current conservation" is completely 
unfounded as can be seen from eq.( 4.51|) which is an obvious consequence of current con- 



servation (cf.e.g. eq.(7.61) in Q). Last but not least the statement that because the 
tensor structure in eq.(20) of ref. is parametrized in terms of ki and k2 separately 
"instead of only in terms of K, leading to spurious terms in the tensor structure which 
eventually result in their spurious momentum-dependent prefactor" has also to be qual- 
ified. As mentioned above, eq.( [4.58| ) to which this observation about the K dependence 



of [B3] refers is not general enough for the problem of photon interferometry. 



4.6 Coherence and lower bounds of Bose-Einstein correlations 

We mentioned in the previous subsections that the intercepts of the second and higher 
order correlation functions can deviate from the canonical values derived within the wave 
function formalism. This effect is important for at least three reasons: 

(i) it illustrates the limitations of the wave function approach 

(ii) it can in principle (provided other effects like final state interactions are taken into 
account) be used for the determination of the degree of coherence. 
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(iii) it can serve as a test of models of BEC, since the value of the intercept follows 
from very general quantum statistical considerations, in particular the Gaussian nature 
of the density matrix. 

In most BEC models the intercept is identical to the maximum of the correlation 
function and therefore it can be studied by limiting the discussion to chaotic sources as 
was done in the previous subsection where the upper bounds of correlation functions were 
investigated. On the other hand the minimum of the correlation functions is determined 
both by the form of the density matrix and the amount of coherence (cf.ref.(|^8[)) because 
coherence leads to a decrease of the correlation function. This will be illustrated below by 
discussing the lower bounds of this function. We will show among other things (cf. 0) 
that in the quite general case of a Gaussian density matrix, for a purely chaotic system 
the two-particle correlation function must always be greater than one. On the other hand, 
in the presence of a coherent component the correlation function may take values below 
unity. Some implications for experimental and theoretical results found in the literature 
will be discussed here as well as in section 5.1.6. 

We have seen in section 4.4 that for identically charged bosons (e.g., vr"*") the two- 
particle correlation function reads 



C++{k,,k2) = 1+ 2^pi(l-pi)-p2(l-p2)Ti2 COs(05^-05+0^) + PlP2T^2 (4-65) 

For neutral bosons like photons, or vr^'s, the terms d{kr, kg) also appear (cf. eq:( [4.36| ) 
in the BEC function: 



C°°(ki,k2) = l + 2V/Pl(l-Pl)-P2(l-P2)T12 COS(05^ - 0J + 0^2) + PlP2Tl^ 



+2y'pi(l-pi)-P2(l-P2) T12 cos(0^^ -<Pl- 0^) + P1P2 fl^ 

(4.66) 

Let us first consider the case of a purely chaotic source. Insertion of p{k) = 1 in eqs. 
( |4.65| ) and ( 4.66| ) immediately yields C2(ki,k2) > 1. In the case of partial coherence. 



the terms containing cosines come into play and consequently C2 may take values below 
unity. Eqs. ( ^ , ^1661) imply that C2~-(ki,k2) > 2/3 and C00(ki,k2) > 1/3. Because of 
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the cosine functions in ([4.65| , |4l66| ) one would expect C2 as a function of the momentum 
difference q to oscillate between values above and below 1. 

Such a behaviour of the Bose-Einstein correlation function has been observed in high energy e+ — 
colhsion experiments (cf., e.g., ref. ||9l[), but apparently not in hadronic reactions. This observation was 
interpreted as a consequence of final state interactions in ref. [75[ |. If final state interactions determine 
this effect, it is unclear why the effect is not seen in hadronic reactions. On the other hand, if coherence 
is responsible for it, this would be easier to understand. Indeed multiplicity distributions of secondaries 
in e+e^ reactions are much narrower (almost Poisson-like) than in pp reactions, which is consistent with 
the statement that hadronic reactions are more chaotic than e+ — reactions |9^p|. 

So far, two methods have been proposed for the detection of coherence in BEC: the 
intercept criterion (C2(k, k) < 2) and the two exponent structure of C2 Both 
these methods have their difficulties because of statistics problems or other effects. The 
observation of C2(ki,k2) < 1 could constitute a third criterion for coherence. 

In the two-particle correlation function has been calculated for photons emit- 

ted from a longitudinally expanding system of hot and dense hadronic matter created 
in ultrarelativistic nuclear collisions . For such a system, the particles are emitted from 
a large number of independent source elements (fluid elements), and consequently one 
would expect the multiparticle final state to be described by a Gaussian density matrix. 
However, although the system is assumed to be purely chaotic the correlation function 



calculated in |Q is found to take values significantly below unity. Clearly, this is in 
contradiction with the general result derived above from quantum statistics (C2 > 1 for 
a chaotic system). 

The reason for this violation of the general bounds derived for a purely chaotic source 
is in this concrete case the use of an inadequate approximation in the evaluation of the 
space-time integrals. However as pointed out in the expression for the two-particle 



inclusive distribution used in ref. (equation (3) of that paper), which in our notation 
^^This statemant is not necessarily in contradiction with the empirical observation that the A factor in 
— reactions appears in general to be larger than in p — p reactions, given the fact that A is not a 
true measure of coherence. 
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takes the form 

P2(ki,k2) = J (fxi j d'^X2 g{xi,ki) g{x2,k2) [1 + cos((A;i - A;2)(a;i - X2))], (4.67) 

is also unsatisfactory Q because for certain physical situations it can lead to values be- 
low unity for the two-particle correlation function even if the integrations are performed 
exactly. To see this, consider, e.g., the simple ansatz 

g{x,k) = const. exp[— «(x — /3k)^] 5(t — to) (4.68) 



where a and (3 are free parameters. The expression for P2(ki,k2) used in ref. |Q then 
yields 



C2(ki,k2) = 1 + exp 



cos[/?q 2] (4.69) 



2a _ 

Clearly, if (3 exceeds the above expression will oscillate and take values below unity. 
On the other hand, in the current formalism (cf. below) one obtains with the same ansatz 
for g 



C2(ki,k2) = 1 + exp 



> 1. (4.70) 



2a 

Thus, eq. ( |4.67| ) can lead to values C2 < 1 if there is a strong correlation between the 
momentum of a particle and the space-time coordinate of the source element from which it 
is emitted. Such correlations between x and k can occur in the case of an expanding source. 
The reason for this pathological behaviour is the fact that the simultaneous specification of 
coordinates and momentum as implied by eq.( |4.67| ) is constrained in quantum mechanics 



by the Heisenberg incertainty relation and any violation of this constraint leads necessarily 
to a violation of quantum mechanics. This violation manifests itself sometimes, as in the 
present case, through a violation of the conservation of probability. This phenomenon 
is also met when using the Wigner function, which for this reason cannot always be 
associated with a bona fide probability amplitude. We will discuss this problem also in 
section 5.1.6. 

The above considerations concerning bounds for the BEC functions refer to the case 

of a Gaussian density matrix. In general, a different form of the density matrix may yield 
■^^This formula appears apparently for the first time in and was criticized (for other reasons) already 
in |9^. It is nevertheless used in certain event generators for heavy ion reactions (cf. section 4.10). 
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correlation functions that are not constrained by the bounds derived here. For instance we 
have seen in section 2.2 that for a system of squeezed states C2 can take arbitrary positive 
values. Moreover, for particles produced in high energy hadronic or nuclear collisions, 
the fluctuations of quantities such as impact parameter or inelasticity may introduce 
additional correlations which may also affect the bounds of the BEC functions. 

4.7 Quantum currents 

The results derived in the previous section, in particular the isospin dependence of BEC, 
were obtained in the assumption that the currents were classical. The question arises 
up to what point these conclusions survive in a fully quantum treatment of the problem. 
It would also be important to get a more precise criterium for the phenomenological 
applicabilbity of the classical assumption, besides the no-recoil prescription. This question 
was discussed in where it was found that the "surprising" effects not only persist when 
the currents are quantum, but that they can serve as an experimental estimate of the size 
of the quantum corrections. We shall sketch briefly in the following the results of ref . ( [|96[ ) . 
As in the classical current case one starts with the interaction Lagrangian 

L,„i(x) = J(+)(x)7r(-)(a;) + J(_)(a;)7rW(a;) + Jo{x)iv\x) (4.71) 

The currents J^^\ J^~\ are operators which we assume again for simplicity not to 
depend on the vr*^^) and 7r° fields ([^7, tt] = 0). Taking into account the different isospin 
components in eqs.( [4.10" ) and ( 4.11| ) we find that as in the classical current case the single 



and double inclusive cross sections depend on these components, e.g.0 



G[-\k) = Tr[p,J^^\-k)J^^\k)} 



(4.72) 



G[-^\k, , fcs) = Tr {p,f [/+\-k^)/-\-k2)] r [/H\k,)ji;\k2)] } (4.73) 



^From now on we shall omit the label H and assume that all operators are written in 



35pQ]- j-gasons of notational simplicity we have replaced Jt(fc) by J{~k). 
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the Heisenberg representation. Assuming a Gaussian density matrix one gets 

G'C\k) = F''{k,k) (4.74) 

G[-\k) = G['^\k) = F'''{k,k) (4.75) 

Gi-\k^,k2) = G[-\k^)G[-\k2) + \F'\k^,k2)\'' (4.76) 

Gr(fci,fc2) = GS°^(fci)GS°)(fc2) + |F"(A;i, A;2)P + |$"(A;i, ^2)!' (4.77) 

Gi-^\k„k,) = G[-\k^)G^^\k2) + \<l>^\-k^,k2)\' (4.78) 

where the functions F and $ are defined for charged particles (upper index ch) and for 
neutral ones (upper index n) as follows: 



F'\k,, 


k2) 






-ki 


)jL'\k2)}:> 


= <jW(-^l)^^"n^2)> 


(4.79) 


<l>^''{ki, 


k2) 






-h 


)4-)(fc2)}» 


= «t{jW(-A;i)j(-)(A;2)}» 


(4.80) 


F'^ik,, 


k2) 








jS\k2)}:^- 


= <j(°H-^i)^^°n^2)> 


(4.81) 


$"(fci, 


k2) 


= <T, 




-kl) 


jS\k2)}:^- 


= <f {j(°)(-A;i)j(°)(A;2)}> 


(4.82) 



In contrast to the classical current approach which deals with only one type of two- 
current correlator we have here two different kinds of two-current correlators depending 
on their ordering prescriptions. Moreover and most remarkably, the difference between 

these two correlators is reflected by the difference between and H — correlations. It 

thus follows from |Q that the "surprising" effects found in i.e. the presence of 



particle-antiparticle Bose-Einstein type correlations and a new term in the Bose-Einstein 
correlation function for neutral particles are reobtained, but under a more general form 
which contains also the quantum corrections. These equations also prove that the above 
effects are not an artifact of the classical current formalism but have general validity. 

Morover and most remarkably, from the above equations follows that the difference 
between the effects of the classical and quantum currents resides in just these "new effects" 
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and in particular in the difference between 00 and correlations i. e. in the H — corre- 
lations. This result can serve as an estimate of the importance of quantum corrections 
to the classical current formalism of BEC. Since H — correlations are in general small, it 
follows that the classical current approach is a good approximation, except for very short 

lived sources, where the H — correlations become comparable to the correlations. It 

also follows that the experimental measurement of H — correlations is a highly rewarding 
task, since they are a rather unique tool for the investigation of two very interesting effects 
in BEC, namely squeezed states and quantum corrections. 

4.8 Space-time form of sources in the classical current formalism 

In two types of sources were considered, a "static" one which corresponds to a source 
in rest and an expanding one. We will present below some of the results, as they exem- 
plify certain important features of the space-time approach within the classical current 
formalism. 



A static source The space-time distributions of static sources, as well as the primordial 
correlator, are parametrized as Gaussians: 



/c/i(a;)=exp 



X, 



R2 



ch,. 



(4.83) 



and 



/c(x)=exp 



Xr] 



c,0 



xj 

r! 



Rl,± 



C{x-y) 



exp 



{xo-y^f (x||-?/||)2 (xx-y±)^ 



(4.84) 



(4.85) 



'2iLq ^^11 '2iLj^ 

Note that the term static here does not imply time independence but rather a specific 
time dependence defined by eqs. ([4.83|) , (|4.84|) corresponding to source elements being 
at rest. This is to be contrasted to the expanding source, discussed in the next section, 
which explicitly contains velocities of source elements. 

The main justification for this particular form of parametrization is mathematical 
convenience, because, as will be shown below, for this case the correlation functions in 
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momentum space can be calculated analytically and the physical implications can be read 
off immediately. 

In eqs.( |4.83| - f4.85| ), Rch,a and Rc,a{c( = 0, ±, ||) are the lifetimes, transverse radii and 
longitudinal radii of the chaotic source and of the coherent source, respectively, and 
La{<y = 0, _L, II) are the correlation time and the corresponding correlation lengths in 
transverse and in longitudinal direction. The relative contributions of the chaotic and the 
coherent component are determined by fixing the value of the (momentum dependent) 
chaoticity parameter p at some arbitrary scale (in this case, at k = 0): 



Po = P{k = 0) 



(4.86) 



The model contains 10 independent parameters: the radii and lifetimes of the chaotic and 
of the coherent source, the correlation lengths in space and time, and the chaoticity po- 
In it is assumed that L\\ = L± = L, i.e., that the medium is isotropic, which leaves us 
with 9 independent parameters. 
With the definitions 



d2 t2 



[a 



0, 



(4.87) 



one may write the single inclusive distribution in the form 



E 



a 



where 



and 



a (Pk 



Sch{k) 



Sc{k) 



E 



a 



a (Pk 



(Po Sch{k) + (1 ~Po) Sc{k)) 



exp 



exp 



fc=0 



k2 R2 



k2 R2 



(4.89) 



(4.90) 



2 2 2 

The scales which determine the mean energy-momentum of the coherently produced par- 
ticles are given by the inverse lifetime and radii, Rc\i of the coherent source. For the 
chaotically produced particles, these scales are given by the inverse of a combination of 
correlation lengths and dimensions of the chaotic source, R^i- Eq.( |4.87|) implies that 
RaL < Rch a- The radius of the chaotic source enters the single inclusive distribution only 
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in combination with the correlation length L. This feature which occurs also for higher 
order correlations leads to the important consequence that experimental measurements of 
BEC do not provide separately information about radii (lifetimes) of sources, nor about 
correlation- lengths (-times), but rather about the combination of these quantities as given 
by eq. ( [4.87| ). On the other hand, by measuring both the single and the double inclusive 
distribution one can determine radii and correlation lengths separately. 

It follows from eqs.( [4.88 , 4.8£ , 4.90 ) that in the presence of partial coherence in general 
(i.e., unless Rai = Rc,a) the single inclusive distribution is a superposition of two Gaus- 
sians of different widths. If the geometry of the coherent source is the same as that of the 
chaotic source, one has i?c,a = Rch.a > Ria, which would imply that coherently produced 
particles can be observed predominantly in the soft regime. However, if the coherent radii 
are small compared to the chaotic ones, this situation is reversed. 

As a next step, consider the correlation functions. The correlation function of two 
negatively charged pions is 



C2 {ki,k2) = 1+ 2Jpi{l-pi) ■p2{l-p2) T12 cos( 



p,P2 (4.91) 



For the Gaussian parametrizations all phases in the second order correlation function 
disappear 

0^2 = 0?2 = 0,' = 0,^ = (4.92) 

and 



Ti2=exp 



8 



8 



(ki^-k2^)2(i?,V-i?ij 



(4.93) 



Ti2=exp 



(^1 + E^nRLo - RIl) iKw + k^MRl,. - R\ 



+ k2±) {Rchi. - R±L 



(4.94) 



■^^This is not the case anymore for an expanding source e.g. (cf. below) or in general for higher order 
correlations 
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The two particle correlation function C2 is the sum of a purely chaotic term (oc TI2) 
and a mixed term (oc T12). The momentum dependence of the chaoticity parameter, 
p = p{k), implies a momentum dependence of the contribution of the mixed term relative 
to that of the purely chaotic term. To see how this affects the interplay between the two 
terms (i.e., the interplay between the two Gaussians), it is useful to explicitly insert the 
momentum dependence of the chaoticity parameter by writing 

Pr = p{kr) = ^ (r = l,2) (4.95) 



with 



and 



Ar = A{kr) = po + (1 -Po) Srr (r = 1,2) (4.96) 



(4.97) 



With this, 6*2 takes the form 



C--(.,, = 1 + -"^^ r,. + (4.98) 

The momentum-dependence of the relative contributions of the purely chaotic and of 
the mixed term is reflected in the factor 5*12. Depending on the sign of the combinations 
-^c,a ~ -^aL5 'S'12 may act either as a suppression factor or as an enhancement factor of 
the mixed term relative to the chaotic term. This is a consequence of the fact that, in 
contrast with the case of the correlation function C2 derived within the wave function 
formalism, where the C2 depends only on the difference of momenta ki — k2 now the 

correlation function depends also on fci + /c2- 

■^^Up to recently this desirable physical property, which is observed in most experimental data on BEC, 
was considered to be a consequence of the expansion of the source and used to be derived within the 
Wigner function formalism, which is also a particular case of the classical current formalism.) As shown 
in the example treated above (cf. jss)) it can be considered also a consequence of the (partial) coherence 
of a non-expanding source. 
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It is instructive to discuss the tilde terms that give rise to the particle-antiparticle 
correlations for the parametrization ( [4.83 - C85 ) of a static source. For the sake of trans- 



parency, consider only the purely chaotic case, po = 0. The correlation functions of like 
and unlike charged pions then take the form 

C2'{ki,k2) = l+T^2 (4.99) 

C+-{h,k2)=l+f^2 (4.100) 



^From (|4.93|) and ([4 .941) it can be seen that the "new" T12 terms that appear in the 



particle-antiparticle correlations are in general small compared to the "ordinary" T12 
terms that determine the particle-particle correlations. The term T12 gives rise to an 
ant icor relation effect due to the factor in eq. ( [4.94|) containing the sum ki + k2, if the first 



factor, containing Ei + E2, is not too small. The latter is possible, if the time duration 
of the pion emission process and/or the pion energies are small enough. We thus expect 
an enhanced contribution of the "new" terms for soft pions. The appearance of anti- 
correlations is, as mentioned above, a general property of squeezed states, which are 
present in the space-time formalism of 0. The tilde terms arise as a consequence of the 
non-stationarity of the source. In the limit of a stationary source, Rq — ^ 00, and T12 0. 
An upper limit is given by 

Ct = l + |Ti2r < l + exp[-{Rl-RlL)ml] (4.101) 

In the limit Lq » Rq on the other hand, Rq ~ Rql and reaches its maximum 

value {C2~)max = 2. We observe in the above equation that the contribution of the 
"tilde" terms increases with decreasing mass of the particles and reaches its maximum 
of 2 for massless particles at fixed and non- vanishing Rq — Rql- This is related to the 
observation made in the case of photon BEC where we saw that for unpolarized photons 
the maximum of C2 is also 2. Indeed the role of the charge degrees of freedom ((+,—) for 
pions is for unpolarized photons taken over by the spin. From this mass dependence or 
more general from the energy dependence of the anticorrelation follows that if the factor 
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Some authors have recently called them "back to back" correlations. 
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E1+E2 in eg. ([4 .941) could be decreased, an enhancement of the "new" terms would emerge. 
A possible mechanism for this could be the sudden transition mechanism considered in 
p!6| . Indeed as shown in this reference, for a chaotic source the correlator < a(ki)a(k2) > 



characteristic for the "new" terms turns out to be an increasing function of the parameter 
T = \ log{Ea/ Eb) where Ea, Eb are the energies of the particle in the vacuum and medium 
respectively. Thus by allowing for medium effects, which in a certain sense is equivalent 
to an effective change of mass 0, one can possibly enhance the anticorrelation effect Qin 
BEC. 



Expanding source High energy multiparticle dynamics suggests that the sources of 
produced particles are expanding. This property is reflected in particular in hydrody- 
namical models and also in string models. In terms of the current formalism this means 
that the correlators are velocity dependent. While many of the studies of Bose- Einstein 
correlations for expanding sources have followed, with slight variations, a Wigner func- 
tion type of approach, the use of the Wigner approach is in general too restrictive and 
is recommendable only in the case when a full-fledged hydrodynamical description of the 
system is performed. In the present section, following |^, we shall therefore start with a 
more general discussion of the expanding source which is based on the space-time current 
correlator and the space-time form of the coherent component and which is not affected by 
the semi-classical and small q approximations inherent in the Wigner function approach. 
We introduce the variables r, rj and x\\, with 

r = Jxl-xl 7] = iln^^^ (4.102) 

V " 2 xq — x\\ 

Here r is the proper time, x\\ the coordinate in the longitudinal direction (e.g. the collision 
axis in p-p reactions or the jet axis in — reactions) and rj the space-time rapidity. An 
ansatz which is invariant under boosts of the coordinate frame in longitudinal direction 
will be considered (Fig. 8). Physically this ansatz is motivated by the prejudice that the 
^^This particular possibility was suggested in Andreev [|8| suggested a time evolution scenario 

for the medium effect, which involves two modes k and —k. 

Anticorrelations in disoriented chiral condensates are considered in . 
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single inclusive distribution in rapidity is flat. The space-time distributions of the chaotic 
and of the coherent source and the correlator are then parametrized as 



/,,(x)~exp ( 1 exp ( ) (4.103) 

/.(x)~exp exp f-il (4.104) 



exp (-^^ 



C{ti - T2, rji - r]2, x_L,i - x^,2)=exp 



ilL^-^sinh^^^^-^^' 



[X^,l - X_L,2j 

2Ll 



(4.105) 



The model contains again 10 independent parameters: The proper time coordinates 
of the chaotic and the coherent source, ro,ch, to,c, their widths in proper time, Stch and 
5tc, the transverse radii, Rch and i?c, the correlation lengths Lr, L± and L^, and the 
chaoticity parameter po- 

In order to be able to obtain explicit expressions for the single inclusive distribution 
and the correlation functions, a further simplifying assumption is made, namely, that 
5Tch = 6tc = 0. Eqs.(|]l03|,|]T04D then take the form 



fch{x)r^S{T - TQ^ch) exp 



(4.106) 



fc{x)r^6{T - ro,c) exp - 



(4.107) 



Now the results no longer depend on the correlation length L^-, and one is left with 7 
independent parameters: tq^c/i, to.c, Rch, Rc, L±, and Pq. 

The Fourier integrations necessary to obtain D{ki, and I{k) can be performed by 
doing a saddle point expansion; this should provide a good approximation if 



»1, 



2^ 



> 1 



With the definitions: 



Rlh + L'i 



(4.108) 



(4.109) 
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7l2 



ro,chimi± - m2±) 



ro,chi'mi± + m2±) 

7l2 = 



(4.110) 



the single inclusive distribution can be written as the sum of a chaotic and a coherent 
term, 



E 



a 



a d?k 



{Po Sch{k) + (1 -Po) Sc{k)) E 



a 



(7 d^k 



k=0 



with 



Schik)-- 



Sc{k)-- 



m± 



exp 



exp 



2 



(4.111) 



(4.112) 



(4.113) 



where 'm± is the transverse mass of the pions emitted. The momentum dependence of the 
chaoticity parameter takes the form 



Pr = p{kr) = ^ 



1,2) 



with 



and 



A{kr) 



exp 



Po + (1 -Po) Srr (r = 1,2) 



(4.114) 



(4.115) 



(4.116) 



Unless Rc — Rl, the transverse momentum distribution is a superposition of two Gaus- 
sians of different widths. The rapidity distribution is uniform, dN/dy — const., as a 
result of boost-invariance. In opposition to what is assumed usually in simplified pseudo- 
hydro dynamical treatments, the transverse radius of the chaotic source, i?c/i, cannot be 
determined independently by measuring only the single inclusive distribution, as the quan- 
tity Rl which sets the scale for the mean transverse momentum of the chaotically produced 
particles is a combination of Rch and the correlation length L^. 
We recall that the second order correlation functions 

C2++(ki,k2) = 1+ 2^pi(l-pi)-p2(l-p2)ri2 COs(05^-05 + 0^) +PlP2Tl^ (4.117) 

are defined in terms of the magnitudes and phases of drs and dre-i Trs, Trs, 4>r'l and (pf^ of 
the chaotic source as well the phases of the coherent component, 0^. The expressions for 
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these quantities read for an expanding source: 



12- 



12- 



<l + 7i2 
<l+7i2 



-1/4 



exp 



-1/4 



exp 



1 + 7l2 

b 

1 + 712 



2 (yi - 2/2)^ - 



(yi - 2/2)^ 



and 



6712 



^?2 = 



1 + 7l2 
&7l2 



(2/1 - 2/2)^ - ro,c/^("^l^ 



"^2± 



7l2 

^J = - ^0,c"^j± 



2/2 J 



arctan7i2 



^o,ch("^i± + ?«2±) - -arctan7i2 



(4.118) 

(4.119) 

(4.120) 
(4.121) 



One thus finds again that the correlation functions do not depend separately on the 
geometrical radii R or on the correlation lengths L but rather on the combination Ri 
defined in ( [4.109| ). This expression reduces in the limit Rch ^ L to L and in the limit. 



Rch <^ L to R. The model considered in |^ is thus a particular case of the space-time 
approach [|| for L = 0. 

As in the static case the tilde terms give rise to the particle-antiparticle correlations. 
For a purely chaotic system the intercept of the tt'^ti^ correlation function is 



C2 {k, k) 



1 + 4 



2\ -1/2 



1 



1 + 4 



/ To,ch \ 



2\ -1/2 



(4.122) 



We conclude this section with the observation that in a correspondence between 
the correlation length L in the primordial correlator C{x — y) and the temperature T 
for a pion source that exhibits thermal equilibrium was established. In the limit of large 
volume V (X R^ and lifetime Rq of the system, it reads 



L 



(4.123) 



4.9 The Wigner function approach 

As mentioned previously, the experimental observation of the fact that the two particle 
correlation function depends not only on the difference of momenta q = ki — k2 but 
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also on the sum ki + k2 led to the introduction and the use ||100|| of a "source" function 



within the well known Wigner function formalism of quantum mechanics While it 
turned out later that this property of the correlation function can be derived within the 
current formalism without the approximations involved by the Wigner formalism, this 
formalism is still useful when applied within a hydro dynamical context. On the other 
hand it should be clear that there is no justification for using this formalism without a 
full hydrodynamical apparatus (cf. also below). This message has apparently not yet 
come through, as many phenomenological and experimental papers have continued to use 
this restrictive formalism. 

The Wigner function approach for BEG was proposed in a non-relativistic form in 
Ref. ||1UU[ and subsequently generalized in |101[ ], |0(cf. also |B^,|^). The Wigner function 



called also source function, g{x,k), may be regarded as the quantum analogue of the 
density of particles of momentum k at space-time point x in classical statistical physics. 
It is defined within the wave function formalism as 

g{^, k, t)= J d'x'r (x + ^x', ^ (x - ix' , t) e'"^^' 

= J d^k'ij* (k + ^k', t^^(k- ik', ?j e"^'"" (4.124) 
and is related to the coordinate and momentum densities by the relations 

n(x,t) = J d^kg{x,k,t) (4.125) 

and 

n(k,t) = J d^xg{x,k,t) (4.126) 

respectively. 
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An attempt to consider the correlation between coordinates and momentum was also performed 



earlier within the ordinary wave function formalism by Yano and Koonin |94 who proposed a formula 



for the second order correlation function of form (4.67). However this form turned out subsequently to 



have pathological features as it leads in some cases to a violation of the lower bounds of the correlation 
function (cf. section 5.1.6). The reason for this mis-behaviour was mentioned in section 4.6 and will also 
be discussed in the following. 
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Due to its quantum nature the function g{x, k) takes real but not necessarily positive 
values. Although eq. ( [4.124| ) is nothing but a definition which does not imply any approx- 
imation, its form suggests that it might be useful when simultaneous information about 
coordinates and momenta are desirable, provided of course that the limits imposed by 
uncertainty relations are not violated. As a matter of fact as will be shown below, the 
Wigner function is useful for BEC only if a more stringent condition is fulfilled, namely 
that the difference of momenta q of the pair is small, as compared with the individual 
momenta of the produced particles. It is thus clear that its applicability is more re- 
stricted than that of the classical current approach, where only the "no recoil" condition, 
i.e. small total momentum of produced particles, as compared with the momentum of 
incident particles, must be respected. This circumstance is often overlooked when com- 
paring theoretical predictions based on the Wigner approach with experimental data. In 
particular herefrom also follows that the application of the Wigner formalism to data has 
necessarily to take into account from the beginning resonances which dominate the small 
q region. It turns out that the use of the Wigner function for BEC is justified only in 
special cases as e.g. when a coherent hydrodynamical study is performed, i.e. when the 
observables are related to an equation of state and when simultaneously single and higher 
order inclusive dstributions are investigated. Unfortunately only very few papers, where 
the Wigner function formalism is used, are bona fide hydrodynamical studies. The ma- 
jority of "theoretical" papers in this context are "pseudo-hydrodynamical" (cf. sections 
5.1.5 and 5.1.6) in the sense that the form of the source function is expressed in terms 
of effective physical variables like temperature or velocity, which are not related by an 
equation of state. In this case the application of the Wigner approach is a "luxury" which 
is not justified. This is a fortiori true since, as will be shown in the following, the Wigner 
approach is mathematically not simpler that the classical current approach, of which it is 
a particular case. Thus the space-time model presented above (cf. section 4.8) is more 
general than the Wigner approach, albeit it is not more complicated and has not more 
independent parameters 

In second quantization g{x, k) is defined in terms of the correlator < a"''(ki)a(kj) > by 
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the relation 

< a\k,)a{kj) >= J d^x exp[-ix^{k,'' - k/)] ■ g[x, l{ki + kj)] (4.127) 



This is a natural generalization of (|4.126|) to which it reduces in the limit ki = kj. 



Accordingly, for the second order correlation function one writes 

(4.128) 

where K^^ = (fcf + k2)/2 and = k^ — fcj are the mean momentum and momentum 
difference of the pair 

The relation between this Wigner approach and the classical current approach is es- 
tablished by expressing the rhs of eq. (|4.127|) in terms of the currents. One has 



The derivation of the Wigner formalism from the classical current formalism has the 
important advantage that it avoids violations of quantum mechanical bounds as those 
mentioned previously. 

Note that in the rhs of eq.( [4.128D enters the off-mass shell average momentum \{ki + k2) 



which is not equal to the on-mass shell average K = E"^ — ml + where E is the 

total energy of the pair (1,2). This means among other things that in this approach it 

is not enough to postulate the source function g in order to determine the second (and 

higher order) correlation function C2, but further assumptions are necessary. Usually oe 

neglects the off-mass shellness i.e. one approximates E by the sum Ei + E2 where Ei 

are the on-shell energies of the particles (1,2), which means that one neglects quantum 

corrections Q which is permitted as long as ki — k2 = q is small 

''^For neutral particles, there are additional contributions to P2(ki,k2) which play a role for soft 
particles and which will be neglected here. 
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That these corrections can be important has also been shown in [ 102 



^''it is sometimes argued that the relevant q range in BEC is given by D^^ where Z? is a typical length 
scale of the source and therefore for heavy ion reactions this should be allowed. This is not quite correct, 
because the shape of the correlation function from which one determines the physical parameters of the 
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As mentioned already, the use of the Wigner formahsm is worthwhile within a true 
hydrodynamical approach when the relation with the equation of state is exploited. In 
this case the probability to produce a particle of momentum k from the space-time point 
X then depends on the fluid velocity, and the temperature, T(x), at this point, and 

one has 



EiEj < a^{ki)a(kj) > 



(2vr) 



exp 



exp[-zx^(fci^ - A:/)] (4.130) 



- 1 



Here, da^ is the volume element on the freeze-out hypersurface S where the final state 
particles are produced. We will discuss applications of this approach in section 5. 



4.9.1 Resonances in the Wigner Formalism 

For a purely chaotic source, the formalism to take into account the effects of resonance 
decays on the Bose-Einstein correlation function can be found, e.g. in [|1^, |5^. An 
extension of this approach is due to |^ and ||101|| which allows to consider also the 
effect of coherence and provides rather detailed and subsequently, apparently, confirmed 
predictions for heavy ion reactions. It is based on the Wigner function formalism. 

The correlation function of two identical particles of momenta ki and k2 can be written 

as 



C2(ki,k2) = 1 + 



Al2 A21 



(4.131) 



An A22 

where the matrix elements Aij are given in terms of source functions g{x, k) as follows 



A, = Je,Ej < at(k,)a(k,) > = / d'xg{x^, P) e'^ 



(4.132) 



A typical source function reads 



(4.133) 



res=p,ui,r),. 



source is not given just by the values of the correlation function near the origin, but depends also on its 
values at large q. 
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where the labels dir and res — > tt refer to direct pions and to pions which are produced 
through the decay of resonances (such as p, u, r], ... etc.), respectively. 

The contribution from a particular resonance decay is estimated in ||56|j,| ]101[ using 
kinematical and phase space considerations as well as the source function of that res- 
onance. The source distribution for the direct production of pions and resonances is 
calculated assuming local thermodynamcal and chemical equilibrium as is appropiate for 
a hydrodynamical treatment. 



Pfnx.p") ^ ^ / ,^f^'t-'^:!rr£ (4.134) 



(27r)3 



exp 



TfK) 

Here a denotes the particular resonance and da^ is the differential volume element and 
the integration is performed over the freeze-out hypersurface S. and Tf are the 

four velocity of the fluid element at point x and the freeze-out temperature, respectively. 
B and 5* are the baryon number and the strangeness of the particle species labeled a, 
respectively, and hb and ns are the corresponding chemical potentials. J is the spin of 
the particle. 



This approach is then extended ||101|| to include also a coherent component rsulting in 



a second order correlation function of the form 

C2(ki,k2) = l + 2p,ff (1 -Pe//) Red^^+p%f |di2p (4.135) 
where p^ff is an effective chaoticity related to the true chaoticity pdir via 

Peff = Pd^r{l-m + f'- (4-136) 

f^^^ is the fraction of particles arising from resonances. The form of this equation us 
the same as that derived previously for a partially coherent source within the current 
formalism and manifests the characteristic two component structure. 



45 



It is assumed that only directly produced particles have a coherent component. 
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The sensitivity of the correlation function on the chaoticity parameter pdir can be 
estimated e.g. from the intercept (cf. eq. ( ^4.135|) ) 

h = C2(k,k) = l + 2pe//-Pe// (4.137) 

Defining the fractions of pions produced directly (chaotically and coherently) and from 
resonances, 

fdir Pdir -^ii rdir i}- ~ Pdir) rres ^res=p,u,r),... qoN 

Jch = — A ) Jco = 2 ) / = 2 (4.id8j 



with 



A- 



Jch ~r Jco + J — ^ 

In fig. 9 the intercept of the correlation function is shown as a function of pdir and 
Z^*^**. In order to read off the fraction of direct chaotically produced particles, pdir, from 
the intercept of the correlation function, one has to extract the effective chaoticity Pejf 
according to eq. ( [4.1371 ) and then correct for the fraction of pions from resonance decays. 
Note that pe// < Pdir- In particular, if a large fraction of pions arise from resonance 
decays, Pe// 1 and it will take very precise measurements of the two-particle correlation 
function at small q to determine the true chaoticity, pdir- 

A further complication arises if a fraction of particles are the decay products of long- 
lived resonances This topic as well as the problem of misidentification are discussed 
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■^^In some papers [ 103 pions originating from long lived resonances are associated with a "halo" while 
those coming form short lived resonances or directly produced are related to a "core" . Then it is claimed 
among other things that the "core" parameters of the source (like radius and A factor) can be obtained 
from the data just by eliminating the small Q points and fitting only the remaining points. Even if such 
a separation would be clear cut (there are doubts about this because of the lu resonance) , it would be of 
course dependent on the resolution of the detector. 
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4.10 Dynamical models of multiparticle production and event 
generators 

Due to the lack of a full fledged theory of multiparticle production in strong interactions 
different models of multiparticle dynamics were proposed. Bose-Einstein correlations 
measurements have been used either to test a particular model or/and to determine 
some of its parameters. Among other things these models might be used to predict the 
dependence of the chaoticity on the type of reaction. In the following we will sketch the 
main theoretical ideas on which these models are based and mention briefly their relation 
to data. 

One of the first models of particle production from which definite predictions on BEG 
can be derived is the Schwinger model | 104 | for — e~ reactions. It visualizes the 
source as an one-dimensional string in a coherent state and thus predicts the absence of 
any bunching effect. A similar prediction follows from the bremsstrahlung model ||105|| . 
Recoilless bremsstrahlung can be described by a classical current which also corresponds 
to a coherent state. Given the fact that in all hadron production processes BEG have 
been seen, it follows that the above two models are ruled out by experiment. 

More complex predictions follow from a dual topological model due to Giovannini and 
Veneziano [|106|| which associates the processes — ^ hadrons to a unitarity cut in 
one plane, reactions due to Pomeron exchange to a cut in two planes, and annihilation 
reactions p — p to a cut in three planes. This model predicts then among other things 
that for TT^TT^ BEG the intercepts C2{k,k) of the second order correlation functions for 
the above reactions should satisfy the following relation: 

[Cf^-{k,k)-2]/[Cr{k,k)-2]/[Crik,k)~2] = (4.139) 

(A similar, but quantitatively different relationship is predicted for vr+vr" correlations.) 

Despite the fact that since the publication of this paper in 1977 many experimen- 
tal BEG studies of these reactions were performed, the above predictions could not yet 
be tested quantitatively in a convincing manner. This is due among other things to 
experimental difficulties (cf. below, in particular factors (a),(b)) and illustrates the un- 
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satisfactory status of experimental BEC investigations. A qualitative remark can however 
be made: The expectation that the annihilation reaction leads to more bunching than 
other reactions is apparently confirmed (cf. e.g. ref.[|T2] and section 2.1.1). As to the 



difference between — e~ reactions and hadronic reactions the experimental situation is 
rather confused (cf. also below). 

A somewhat related dynamical model based on Reggeon theory, was proposed already 
in [107]. A straightforward extension of of this formalism to heavy ion reactions does not 



work as it predicts that the longitudinal radius is of "hadronic" size [108 



A different approach to BEC based on the classical current formalism is proposed in 



109|[ . The currents are associated with the chains of the dual parton model, and contrary 
to what is assumed in other applications of the classical current formalism, all the phases 
of these elementary currents are fixed, so that the source is essentially coherent. This is a 
special case of the classical current approach presented in sections 4.2, 4.3 and 4.8 where 
allowance is made both for a chaotic and coherent component. The model is intended to 
work for p — p reactions where the authors state that resonances do not play an important 
role. It explains, according to the authors, the dependence of the A parameter in the 
empirical formula for the second order correlation function 

C2 = 1 + A exp(-i?2g2) (4.140) 

on the multiplicity and energy. Unfortunately the claim that in p — p reactions pions are 
only directly produced is unfounded. Furthermore there are other factors which influence 



the multiplicity dependence of A (cf. section 6.2) which are not considered in ||109|[ and 
which are of more general nature. 

An orthogonal point of view for the interpetation of the same A factor (also for directly 
produced pions, only,) is due to [[110|[ . In this approach the source is made up of totally 



chaotic elementary emitting cells which are occupied by identical particles subject to Bose- 
Einstein statistics. Different cells are independent so that correlations between particles 
in different cells lead to A = 0, while correlations between particles in the same cell are 
characterized by A = 1. From the interplay of these two types of correlations, one obtains 
with an appropiate weighting, large A values in e"*" — e~ reactions and small A values in 
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p — p reactions, as in ||109|| , but within a completely different approach. 

We conclude the discussion of these two approaches by the following remarks. Besides 
the reservations about the role of directly produced pions in BEC expressed above and 
which presents the two approaches in a rather academic light, it is unclear whether the 



A factor in e+ — e~ reactions is larger than in p — p reactions as assumed in | 110 |. This 
issue awaits a critical analysis of the specific experimental set-ups. The fact that quite 
different approaches lead to similar conclusions about the A factor confirms that the 
parametrization of the second order correlation function in the form ( ^.140| ) is, as pointed 
out already in section 2.2, inadequate. 

We discuss now other two, closely related, approaches, which make more detailed 



predictions about the form of the correlation function in e"*" — e~ reactions: , |]75[ on 
the one hand, and ||23[| , ||1 1 1|| , ||112|| on the other. Both approaches are based on a variant 
of the string model (for a more extended review of this topic cf. e.g. [|113|| ). Such a 
string represents a coloured field formed between a quark q and an antiquark q, which 
tend to separate. Because of confinement the breakup of the string can be materialized 
only through creation of new qq pairs, which are the mesons produced in the reaction. 
The difference between the Schwinger model of confinement on the one hand and the 
models of Bowler and Andersson-Hofmann-Ringner on the other is that in the former the 
field couples directly and locally to a meson, while in the latter ones the quarks, which 
constitute the meson, are created at different points. This feature destroys the coherence 
inherent in the Schwinger model and makes possible the Bose-Einstein bunching effect. 
For massless quarks the second order correlation function can be approximated by the 
relation |^ 

^2 = 1+ < cos(kAA)/ cosh(6AA/2) > (4.141) 

where AA denotes the difference between the space-time areas of coloured fields spanned 
by the two particles, k, is the string tension and b a parameter characterizing the decay 
probability of the string. For massive quarks the formulae become more involved and 



were approximated analytically in or calculated numerically in |^,[ni], | 112 |. In 



this model the correlation function depends both on the difference of momenta ki — k2 
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as well as on their sum fci + ^2, reflecting the correlation between the momentum of the 
particle and the coordinate of its production point. This is a consequence of the fact that 
string models use a Wigner function type approach. From the above equations follows 
that there are two length scales in the problem, one associated with n and another with 
h. Phenomenologically these correspond to q\\ and g_L. Both these lengths are correlation 
lengths rather than geometrical radii. (As a matter of fact there is no geometrical radius 
in the string model.) Their magnitudes are quite different. 



In both string approaches |7^, p3| one obtains a difference between BEC for identically 
charged and neutral pions as found in |7T|]. However while in there is room for 
coherence, this is apparently not the case for ||112|| , which predict a totally 

chaotic source. Furthermore in p2| an energy dependence of the BEC is predicted (the 
correlation function is expected to shrink with increasing energy), while in p5[, |lll| |,[ |112|| 
the correlation function does not depend on energy. 



A rather discordant note in this string concerto ||22[, [^] is represented by the paper 



by Scholten and Wu [ 114 |. These authors, using a different hadronization mechanism 
conclude that dynamical correlations, at least in e"*" — e~ reactions dominate over BEC 
correlations so that BEC cannot be used to infer information about the size and lifetime 
of the source 0. This point of view seems too extreme, as it is contradicted by some 
simple empirical observations: in — e~ reactions, as well as in all other reactions, cor- 
relations between identical particle are observed which are much stronger than those of 
non-identical ones, the correlation functions are (in general) monotoneously decreasing 
functions of the momentum difference q and in nuclear reactions the "radii" obtained 
from identical particole correlations increase with the mass number of the participating 
nuclei. All these observations are in agreement with what one would expect from BEC, 
which suggests that dynamical correlations cannot distort too much this picture. How- 



"'^To make the model more realistic in ref. [|2^ resonances were included according to the variant of 
the Limd model (JETSET) in use at that time (1986) and agreement with e+ — data was found. 
Subsequently however it was pointed out in [^8| that some resonance weights used in |p3| were incorrect, 
so that the agreement mentioned above is probably accidental. 

'^^What concerns — e~ reactions similar skepticism was expressed by Haywood M. 
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ever a thorough comparison of BEC in different reactions, using the same experimental 
techniques, appears highly desirable. 



Event generators The model was implemented by Sjostrand ||115|| into JETSET 
under the name LUBOEI by modifying a posteriori the momenta of produced pions so 
that identical pairs of pions are bunched according to . This manipulation "by hand" 
violates energy-momentum conservation which was imposed at the beginning in JETSET. 
To compenate for this, the momenta are rescaled so that energy-momentum conservation 
is restored. However this rescaling introduces spurious long-range correlations, which 
bias the BEC. Nevertheless, in general this program leads to a reasonable description of 
the bunching effect in the second order correlation function 0. More refined features of 
BEC which reflect the quantum mechanical essence of the effect, cannot be obtained of 
course. One reason for this, of rather technical nature, is due to the fact that the ad- 
hoc modification of two-particle correlations does not yet include many-body correlations, 
reflected in the symmetrization (or antisymmetrization) of the entire wave-function. 

Another reason of fundamental character is the fact that event generators like any 
Monte Carlo algorithm deal in general with probabilities and therefore cannot account 
for quantum effects, which are based on phases of amplitudes The event generator 
JETSET was further developped by Lonnblad and Sjostrand , ( ||118|| , ||119|| ) and used 
to estimate the influence of BEC on the determination of the mass of W in e+ — 
reactions, a subject of high current interest for the standard model and in particular 
for the search of the Higgs particle. This effect was also studied using different event 



generators in | 120 | and [ 121 |. The argument of Lonnblad and Sjostrand is the following. 
Consider the reaction e"*" — e~ — > W^W~ <ii<i2(l?,<iA when both W-s decay into hadrons. 



Then according to ||122|| the typical space-time separation of the decay vertices of the 



^^Cf. however ref. |116| 

^°For heavy ion reactions the "Quantum Molecular Dynamics" (QMD) model |117 attempts to surpass 

this deficiency by using wave functions rather than probabilities as input. However this model also neglects 

(anti)symmetrization effects and cannot be used for interferometry studies. 
^^It is interesting to mention that for one string Andersson and Hofmann (|2^)proposed a formulation 

of the BEC effect in terms of amplitudes. However this procedure cannot be used to generate events. 
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and the W~ is less than 0.1 fm (at LEP 2 energies) and thus much smaller than a typical 
hadronic radius (~ 0.5 fm. There will thus be a Bose Einstein interference between a pion 
from the and a pion (with the same charge) from the W~ and one cannot establish 
unambigously the "parenthood" of these pions. This prevents then in this model a precise 



determination of the invariant mass of the W-s. In | 118 | algorithms for the inclusion of 
this effect into the determination of the mass of the W are proposed and for certain 
scenarios mass corrections of the order of 100 MeV at 170 GeV cm. energy are obtained. 
However, as emphasized in ||118|| other scenarios with less or no effect of BEC on the mass 
determination of the W are possible. Thus in [ |120| | and | |121| | effects of the order of only 
20 MeV are found. For more details we refer the reader to the original literature. 

This aspect of BEC is interesting in itself as it illustrates the possible applications of 
this effect in electroweak interactions, a domain which is beyond the usual application 
domain of BEC, i.e. that of strong interactions. 

The Lund model was applied also to heavy ion reactions and then extended to in- 
clude BEC (e.g. the SPACER version |^ of the Lund model). The topic of event 
generators for heavy ion reactions is of current interest because of the ongoing search for 



quark gluon plasma. Padula, Gyulassy and Gavin |^ suggested to use for this purpose 
the Wigner function formalism in order to take into account explicitly the correlation 
between momenta and coordinates, as implied by the inside-outside cascade approach. 
This is evidently another way of expressing the non-stationarity of the correlation func- 
tion mentioned above. (An explicit introducton of momentum-coordinate correlations in 
particle physics event generators like JETSET/LUBOEI is not necessary, because the 
non-stationarity is delivered "free house" by the string model used in the LUND genera- 
tor.) 

On the other hand the Wigner formalism may present also another advantage as em- 



phasized more recently by Bialas and Krzywicki ||123|| . This has to do with the important 
difficulty mentioned also above and which is inherent in all event generators, namely the 
probabilistic nature of Monte Carlo methods. The Wigner function has in certain limits 
the meaning of a wave function and thus provides quantum amplitudes. The proposal of 
Bialas and Krzywicki consists then in starting from the single particle distribution no(k) 
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constructed from non-symmetrized particle wave functions as produced by conventional 
event generators and writing the Wigner function 



^(k;x) = (]o(kV(k;x) (4.142) 

where w(k;x) is the conditional probability that given that the particles with momenta 
ki, , ^2, ■■■kn are present in the final state, they are produced at the points xi, ...Xn- Then 
the art of the model builder consists in guessing the probability w(k;x). This may be 
easier than guessing from the beginning the exact Wigner function. For example a simple 
ansatz would be to assume that the likelyhood to produce a particle from a given space 
point is statistically independent of what happens to other particles. This means that 
^(k; x) can be factorized in terms of the individual particles. Implementations of this 
scheme were discussed in |[L24| , ||125|| 

When using the Wigner formalism or any model (like those used in event genera- 
tors) which specifies momenta and coordinates simultaneously, one must of course watch 
that the correlations between momenta and coordinates do not become too strong. This 
apparently has not always been done.Q 

That such a procedure is dangerous since it can lead to unphysical antibunching effects, 
i.e. to violation of unitarity was already mentioned in (cf. section 4.6). This point 
has been reiterated recently e.g. in [ |1271 ],| |128| ,| |129 |. 



Concluding this section one should emphasize that event generators are just an ex- 
perimental tool, sometimes useful in the design of detectors or for getting rudimentary 
information about experimentally unaccessible phase. Often they are however ab-used, 
e.g. to search for "new" phenomena: if agreement between data and event generators 
is found, one states that no "new " physics was found. Such a procedure is injustified, 
because agreement with a model or an event generator is often accidental. Furthermore, 
for the reasons mentioned above event generators cannot be used to obtain the "true" 
correlation function, i.e. they are no substitute for a bona fide HBT experiment (cf. also 
130|| for a critical analysis of transport models from the point of view of interferometry) . 



^^E.g. it is unclear to us whether the "pure" multiple scattering approach of |126| satisfies the above 
constraint. 
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4.11 Experimental problems 

The confrontation of model predictions with experimental BEC data has been hampered 
by two major facts: (i) most models are idealizations, i.e. they use assumptions which 
are too strong. Examples of such assumptions are: neglect of final state interactions, 
boost invariance, particular analytical forms of the correlation functions; (ii) for various 
reasons in almost no BEC experiment so far a "true" or "complete" correlation function 
was measured, i.e. a correlation function as is defined by 

Here P2 and Pi are the double and single inclusive cross sections respectively. Note in 
particular that the single inclusive cross sections in the denominator are defined in terms of 
the same density matrix as the double inclusive in the numerator. This means that these 
cross sections and the corresponding ratios have to be measured for each event separately 
and only after that one is allowed to average over events. What is measured usually 
instead is a function which differs from eq.( |4.143| ) in two respects: the normalization of 
C2 is not done in terms the product of single inclusive cross sections Pi but in terms 
of a "background" double inclusive cross section, which is obtained either by considering 
pairs of (identically charged) particles which come from different events, or by considering 
oppositely charged particles, or by simulating P2 with an event generator which does not 
contain BEC. This means that the structure of the event is falsified and therefore this 
kind of normalization biases the results. Furthermore, one does not (yet) measure the full 
correlation function C2 in terms of its six independent variables, but rather projections of 
it in terms of single variables like the momentum difference q, rapidity difference yi — 1/2 
etc. Last but not least, the intercept of the correlation function, which contains important 
information about the amount of coherence, cannot be really be measured at present 
because (a) one does not yet control sufficiently well the final state interactions which 
contribute to the intercept; (b) its experimental determination implies an extrapolation 
to g = 0. Such an extrapolation can be performed only if the analytical form of the 
correlation function at g > is known, which is not the case. 

For these reasons at present it is difficult to test quantitatively a given model, except 
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when its predictions are very clear-cut. This circumstance hmits certainly the usefulness 
of BEC as a tool in determining the exact dynamics of a reaction. 

5 Applications to ultrarelativistic nucleus-nucleus col- 
lisions 

5.1 BEC, hydrodynamics and the search for quark-gluon plasma 

The use of BEC in the search for quark-gluon plasma is in most cases based on hydrody- 
namics. This is so because the space-time evolution of the system is given by the equations 
of hydrodynamics the solutions of which are different depending whether a QGP is formed 
or not. In this way hydrodynamics also provides information about the equation of state 
(EOS). QGP being a (new) phase it is described by a specific EOS which is different 
from that of ordinary hadronic matter. The proof that this phase has been seen must 
include information about its EOS and thus the combination of hydrodynamics with BEC 
constitutes the only consistent way through which the formation of QGP can be tested. 

QCD predicts that the phase transition from hadronic matter to QGP takes place only 
when a critical energy density is exceeded. To measure this density we need to know the 
initial volume of the system. While via photon intcrfcrometry (cf. section 4.5 and 4.6) 
one can in principle measure the dimensions (and thus the energy density) of the initial 
state, hadron interferometry yields information only about the final freeze-out stage when 
hadrons are created. To obtain information about the initial state with hadronic probes, 
hydrodynamical models have to be used in order to extrapolate backwards from the final 
freeze-out stage where hadrons are created, to the interesting initial stage. The lifetime 
of the system as given by BEC is also an important piece of information for QGP search. 
Indeed, in order to decide whether we have seen the new phase, we have to measure the 
lifetime of the system. Only lifetimes exceeding significantly typical hadronic lifetimes 
(10~^^sec.) could prove the establishment of QGP (cf. below) 
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5.1.1 General remarks about the hydro dynamical approach 

Besides the main advantages of hydrodynamics related to the information about initial 
conditions, freeze-out and equation of state, for the study of BEC in particular hydro- 
dynamics is very useful because it provides also the single inclusive distributions which 
are intimately connected with higher order distributions as well as the weights and the 
space-time and momentum distributions of resonances, which influence strongly the cor- 
relations. The phenomenological applications of the hydrodynamical approach to data 
are however hampered by two circumstances. 

(i) While the ultimate goal of BEC is the extraction of the minimum set of parameters 
which include radii and coherence lengths both for the chaotic and coherent components 
of the source, in practice, mainly because of limited statistics (but also because of an 
inadequate analysis of the data) one had to limit oneself to the determination of a re- 
duced number of parameters, which we will call in the following "effective radii" -Re// and 
"effective chaoticity" pef / . In reality Rej / is a combination of correlation lengths^ L and 
geometrical lengths R as introduced in sections 4.3 and 4.8. Only in the particular case 
that one length scale is much smaller than the other, one can assume that one measures 
a "pure" radius or a "pure" correlation length. For simplicity in the following we shall 
assume that this is the case and in particular we will assume that R ^ L so that, Reff 
reduces to L. This limit might perhaps correspond to what is seen in experiment, if one 
considers the expansion of the system in the hadronic phase. (In the high temperature 
limit L ^ T^^ (cf. section 4.8). 

(ii) The presentation of the data is yet biased by theoretical prejudices. Instead of a 
consistent hydrodynamical analysis, much simplified models are used (cf. section 5.1.5 
where these models are presented under the generic name of pseudo-hydrodynamics) for 
this presentation and therefore to obtain the real physical quantities, one would have to 
solve a complicated mathematical "inverse" problem, i.e. one would have to reconstruct 
the raw data from those presented in the experimental papers and then apply the correct 
theoretical analysis to these. This has not been done so far and even if the statistics would 

^•^In the following the concept of length refers to space-time. 



88 



be sufficient for this purpose, the outcome is questionable because of the difficulties implied 
by the numerics. (That is why it would be desirable that experimentalists and theorists 
perform a joint analysis of the data or at least that the data should be presented in "raw" 
form.) 

The nearest approximation to the solution of the "inverse" problem found in the 
literature, is that of | |131|| based on the application of the HYLANDER code by the 



Marburg group: It consists in fitting the results of the hydrodynamical calculation to the 
Gaussian form used by experimentalists, 

C2(ki, ks) = 1 + A exp[-igji?j - |gL^L - l^de^Le] (5-1) 

and comparing with the inverse width of the correlation function as presented in the 
experimental papers Here = pi^ — P2^, = ^{pi^ + P2^) and q\\ and K^i denote 
the components of q and K in beam direction, and q± and K± the components transverse 
to that direction; qout is the projection of the transverse momentum difference, q_L on 
the transverse momentum of the pair, 2K^, and qside the component perpendicular to 
K^. (For a source with cylindrical symmetry, the two-particle correlation function can 
be expressed in terms of the 5 quantities K\\, K±, gy, qside and qout-) R\\, Rside, Rout are 



effective parameters, associated via eq. (|5.1|) to the corresponding q components. 

Eq.( |5.1| ) is equivalent to an expansion of the correlation function C(q, K) for small q. 
The use of eq.( ^.l|) for the representation of correlations data implies then that one does 
not measure the geometrical radius of the system but the length of homogeneity, which 
means that energy density determinations based on BEC are an overestimate To take 
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In 1 132 1 it was recommended that experimentalists should use the more complete formula 



C2(ki,k2) = 1 + XeM-hm ~ HutRlut ~ ^9Le^Le " 2g„„,g,[ i?^^^,,, ] (5.2) 
where 2(7oMt9||^o„(|| is called "cross" term. For a more detailed discussion of its meaning and dependence 



on the coordinate system cf. ref. |133|. Of course, in view of the defhciencies of the entire phenomenolog- 
ical procedure outlined above and discussed in greater length in section 5.1.5, these details are of limited 
importance. In particular they do not affect the conclusions discussed here. 
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In |134], |135] a distinction is made between the "local" length of homogeneity Lh{x,k) and the 



"hydrodynamical" length L;i(a;) which is the ensemble average of the former. 
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into account the fact that the correlation function depends in general not only on the 
momentum difference q = ki — k2 but also on the sum K = ^{ki + /C2) , the parameters 
R and A are assumed to be functions of K and rapidity + ^2)- 

Hadron BEC refer to the freeze-out stage. This stage is usually described by the 



Cooper- Frye formula |136 



dk (27r)3 J. exp C-^'^;-^' ) - 1 ' ^ ' 

which describes the distribution of particles with degeneracy factor Qi and 4-momentum 
emitted from a hypersurface element da^ with 4- velocity m'^ ^ After the cascading of 
the resonances we obtain the final observable spectra. 

5.1.2 Transverse and longitudinal expansion 

The equations of hydrodynamics are non-linear and therefore good for surprises. An 
illustration of this situation is represented by the realization, to be described in more detail 
below, that the naive intuition about the role of transverse expansion in the determination 
of the transverse and longitudinal radius may be completely misleading. Only a systematic 
analysis based on 3+1 dimensional hydrodynamics could clarify this issue. 

In the present section we will discuss Bose-Einstein correlations of pions and kaons pro- 
duced in nuclear collisions at SPS energies in the framework of relativistic hydrodynamics. 
Concrete applications were done for the symmetric reactions S + S and Pb + Pb at 200 
AGeV. Many of the theoretical results were predictions at the time they were performed. 



^^In most applications particles produced with momenta pointing into the interior of the emitting 
isotherm {jj^da^ < 0) were assumed to be absorbed and therefore their contribution to the total particle 



number was neglected. In ref |137| this effect was indeed estimated to be negligible and recent attempts 
to reconsider it could not change this conclusion. Another effect is the interaction of the freeze-out system 
with the rest of the fluid. This effect can be estimated by comparing the evolution of the fluid with and 
without the frozen-out part. This is done by equating the frozen-out part with that corresponding in the 
equation of hydrodynamics to the case = 0. The fluid parameters are modified by this procedure at a 



level not exceeding 10% [138]. The influence of the freeze-out mechanism on the determination of radii 



via BEC has been discussed recently in several papers; cf.e.g. 139 and references quoted there 
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These predictions were subsequently confirmed in experiment. In [|140|| it was found tliat 
tlie transverse radius extracted from data on Bose-Einstein correlations (BEC) for 0+Au 
at 200 AGev reached in the central rapidity region a value of about 8/m. It was then 
natural to conjecture that this could be an indication of transverse flow ||141| ] , | |142| | , p43| . 
In the meantime the experimental observation in itself has been qualified | |144 ] and it now 
appears that the transverse radius obtained from the BEC data does not exceed a value 
of 4-5 fm (see however ref. ||143|| ). Motivated by this situation in ref. ||131|| an investigation 
Q of the role of 3-dimensional hydrodynamical expansion on the space-time extension of 
the source was performed and compared with a 1+1 d calculation. 

Contrary to what one might have expected it was found that transverse flow does not 
increase the transverse radius. On the other hand, a strong dependence of the longitudinal 
radius on the transverse expansion was estblished. 

Fig. 10 shows two typical examples of the Bose-Einstein correlations as functions of 
gy and q± for the 1- and the 3-dimensional hydrodynamical solution. The dependence 
of C2{q\\) on transverse expansion agrees qualitatively with what one would expect. For 
a purely longitudinal expansion, the effective longitudinal radius of the source is larger 
than in the case of 3-dimensional expansion, which is reflected in a decrease of the width 
of the correlation function (see also in Fig. 11 below). 

On the other hand, the results for C2{q±) were at a first glance rather surprising. 
Naively one might have expected that the transverse flow would lead to an increase of the 
transverse radius, i.e., to a narrower correlation function C2(g_L). However, in Fig. 10 the 
curves that describe the 1-dimensional and the 3-dimensional results are almost identical. 
If anything, one would conclude that the effective transverse radius is smaller in the 
presence of transverse expansion. 

This effect can however be explained if one takes a closer look at the details of the 
hydrodynamic expansion process as investigated in ref. ||131|| . Due to the strong correlation 
between the space-time point where a particle is emitted, and its energy-momentum, 



^^In ref. [145 the dependence of the effective transverse radius on the transverse velocity field was 
investigated for a fixed freeze-out hypersurface. The important effects of transverse expansion on the 
shape and position of the hypersurface are not considered there. 
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the effective radii obtained from Bose-Einstein correlation data present a characteristic 
dependence on the average momentum of the pair, K^. Fig. 11 shows the dependence of the 
effective radii R\\ , Rside and Rout on the rapidity uk and the mean transverse momentum 
of the pair i^^, both for the 1-dimensional and for the 3-dimensional calculation. The 
longitudinal radius R\\ becomes considerably smaller (by a factor of 2-3) if transverse 
expansion is taken into account. For the 1-dimensional case, an approximate analytic 
expression has been derived for the yx- and the i^'j.-dependence of the longitudinal 'radii'" 
in Refs.lll) f^,^^^ (cf. also |5l): 



fl„ = r^-^, (5.4) 

V m_L cosh(?/i^) 

where mj_ = + m_|_2)/2 is the average transverse mass of the two particles, Tj is the 
freeze-out temperature and Tq = {du\\/ dx)^^ is the inverse gradient of the longitudinal 
component of the 4-velocity in the center (at x = 0) 0. In [|131|| one finds that this 
approximate expression describes R\\ {K^, yx) for S + S reactions quite well, both for the 
1-dimensional and the 3-dimensional case (see Fig. 11). However for Pb + Pb reactions the 
same formula fails to account for the data |p.48| |. This is not surprising, because eq. (|5.4|) 
is based on the assumption of boost invariance, i.e. no stopping. This assumption is not 
justified at SPS energies where there is considerable stopping. The inelasticity increases 
with atomic number and this may explain the breakdown of the above formula. This 
exemplifies the limitations of the boost-invariance assumption, an assumption which must 
not be taken for granted but in special circumstances. 



In ||149|| it was proposed to use the information obtained from fitting the single inclusive 
distribution to constrain the parameters that enter into the hydrodynamic description, 
and then to calculate the transverse radius directly. Indeed let S denote the hypersurface 



^^In this reference eq.(5.4) is used to disprove the applicability of hydrodynamics to p — p reactions 
in the ISR energy range (-y/s = 53GeV). Such a conclusion seems dangerous given the approximations 
involved both in the derivation of this formula as well as in the interpretation of the BEC measurements 
in the above reactions. 

^^The expression (5^) for i?|| denotes in fact the length of homogeneity Lh mentioned above. It refers 
to the region within which the variation of Wigner function is small. By definition < R where R is 
the geometrical radius. 
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in Minkowski space on which hadrons are produced. Then one can define e.g. a transverse 
radius 



= (5.5) 

exp[{pM„^-)^)/T]-l 

where u^, T and /i denote the four- velocity, temperature, and chemical potential on the 
hypersurface S, respectively as in eq. (|5.3|) . 

It is interesting to note that this method for the determination of transverse radii 
based on the single inclusive cross sections provides a geometrical radius while the use of 
the second order correlation function provides a coherence length (length of homogeneity). 
Comparing the effective transverse radius i?^ extracted from the Bose-Einstein correlation function to 



the mean transverse radius as calculated directly in [149| according to eq. (5.3) , one finds in |131| that 



the two results agree to an accuracy of about 10%. This conclusion is confirmed and strengthened in a 



more recent study by Schlei |131| for kaon correlations. 

Of course, this approach can be used only if a solution of the equations of hydrody- 
namics is available; with pseudo-hydodynamical methods this is not possible. 

5.1.3 Role of resonances and coherence in the hydrodynamical approach to 
BEC 

This problem was investigated using an exact 3+1 d numerical solution of hydrodynamics 

am. 



m 



The source distribution g{x,k) was determined from a 3-dimensional solution of the 
relativistic hydrodynamic equations. Fig. 12 illustrates the effect of successively adding 
the contributions from p, u, A and t] decays to the the BEC correlation functions of 
directly produced (thermal) tt^ (dotted lines), in longitudinal and in transverse direction. 
The width of the correlation progressively decreases as the decays of resonances with 
longer lifetimes are taken into account, and the correlation looses its Gaussian shape. 
The long-lived rj leads to a decrease of the intercept. 

Pion versus kaon interferometry Ideally, a comparison of pion and kaon interferom- 
etry should lead to conclusions concerning possible differences in the space-time regions 
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where these particle decouple from the hot and dense matter. It was proposed that kaons 



may decouple (freeze out) at earlier times and higher temperatures than pions 151 



Indeed, preliminary results had indicated that the effective longitudinal and transverse 
source radii extracted from vrvr correlations were significantly larger than those obtained 



from KK correlations ||152|| . However, as we have seen in fig. 12, the BEC of pions are 
strongly distorted by the contributions from resonance decay. It was pointed out in ref. 



52 1 in a study based on the Lund string model that such distortions are not present for 



the BEC of kaons, and that consequently for the effective transverse radii one expects 
R_i_{K^) < i?_|_(7r^), even in the absence of any difference in the freeze-out geometry 
of directly produced pions and kaons. These conclusions were confirmed in [p.5CI|| within 
the hydrodynamical approach and one found furthermore that this effect is even more 
pronounced if one considers longitudinal rather than transverse radii. Furthermore the 



interplay between coherence and resonance production which was not considered in 



was studied in ||150|| . There are also some striking differences between |^2| and ||150|| in 
the resonance production cross sections used. 

In fig. 13, BEC functions of vr" (solid lines) and of K~ (dashed lines) are compared, 
at = and kj_ = 1 GeV/c, respectively. The dotted lines correspond to the BEC 
function of thermally produced vr^. It can be seen that the distortion due to the decay 
contributions from long-lived resonances disappear only at large k±. Fig. 14 shows the 
effective radii R\\, Rside and Rout as functions of rapidity and transverse momentum of 
the pair, both for tt^tt^ (solid lines) and for K''K~ pairs (dashed lines). For comparison, 
the curves for thermally produced pions (dotted lines) are also included. The effective 
longitudinal radii extracted from tt^tt^ correlations are considerably larger than those 
obtained from K~K~ correlations. In the central region the two values for R\\ differ by a 
factor of ~ 2. For the transverse radii, the factor is ~ 1.3. A comparison between results 
for K~ and thermal vr" shows that part of this effect can be accounted for by kinematics 
(the pion mass being smaller than the kaon mass; see also eq. ( |5.4| )). Nevertheless, the 
large difference between the widths of pion and kaon correlation functions is mainly due 
to the fact that pion correlations are strongly affected by resonance decays, which is not 
the case for the kaon correlations. In the hydrodynamic scenario of ref. about 50% 
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of the pions in the central rapidity region are the decay products of resonances ||149| 
while less than 10% of the kaons are created in resonance decays {K* Kit dominates 
contributing with about 



In ref. ||150|| the problem of coherence within the hydrodynamical approach to BEC was 
also investigated. Fig. 15 shows the 7r~7r~ correlation functions in the presence of partial 
coherence. In order to extract effective radii from Bose-Einstein correlation functions in 
the presence of partial coherence, eq.( |5.1|) must be replaced by the more general form 



C2(ki,k2) =1+ A - plff exp -- ^(gi?) 



+ v^- 2peff (1 - Peff) exp - - Y^i^l^f 



5.1.4 Comparison with experimental data 



(5.6) 



Some of the predictions made in for S* + 5* reactions could be checked experimentally 

in refs. [|153||||154|| , in particular the rapidity and transverse momentum dependence of 

radii and remarkable agreement was found. In | 155 | the hydrodynamical calculations 

were extended to Pb + Pb reactions and compared with S + S reactions and, where data 

were available, with experiment. The calculation of Bose-Einstein correlations (BEC) 

was performed using the formalism outlined in refs. ||156| , |150| , |131j| including the decay 

of resonances. The hadron source was assumed to be fully chaotic. 

Figs. 16 and 17 show the calculations for the effective radii R\\, Rside and Rout as 

functions of rapidity yx and transverse momentum K± of the pion pair compared to 

the corresponding NA35 and preliminary NA49 data ||157| , [L58|| , respectively. All these 

calculations, which in the case of S* + 5* had been true predictions, agree surprisingly well 

with the dataQ. This suggests that our understanding of BEC in heavy ion reactions has 

made progress and confirms the usefulness of the Wigner approach when coupled with 
^°The EOS used in the hydrodynamical studies quoted above included a phase transition from QGP 
to hadronic matter. How critical this assumption is for the agreement with data is yet unclear and 
deserves a more detailed investigation. On the other hand the very use of hydrodynamics is based on 
local equilibrium, and this equilibrium is favoured by the large member of degrees of freedom due to a 
QGP. 
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full fledged hydrodynamics. An important issue in comparing data with theory is the 

detector acceptance of a given experiment. This is also discussed in detail in [|155|| . 

Another application of hydrodynamics to the QGP search in heavy ion reactions is due to Rischke 



and Gyulassy |159| who investigate the ratio r = Rout/ Rside- Based on considerations due to Pratt 
|100|, this quantitiy had been proposed by Bertsch and coUaborators ]16(| as a signal of QGP. Under 
certain circumstances one could expect that for a long lived QGP phase r should exceed unity, while for 
a hadronic system, due to final state interactions, the out and side sizes should be comparable. 



The authors of |15!;] performed a quantitative hydrodynamical study of r in order to check whether 
this signal survives a more realistic investigation, albeit they did not take into account resonances. For 
directly produced pions it is found that r indeed reflects the lifetime of an intermediate (QGP) phase. 



However we have seen from |56|] , [ 155 ] that for pion BEG, when resonances are considered, hydrodynamics 
with an EOS containing a long lived QGP phase, leads (in agreement with experiment) to values of r of 
order unity. 



To avoid this complication, in |161| it was proposed to consider kaons at large k±. However even this 
proposal may have to be qualified, besides the fact that it will be very difficult to do kaon BEG at large 
fcj^ . Firstly one has to recall that the entire formalism on which the r signal is based and in particular the 



parametrization (5.1) is questionable. Secondly it remains to be proven that this signal survives if one 
imposes simultanoeusly the essential constraint due to the single inclusive distribution^. Furthermore 
it is unclear up to what values of fc^ the Wigner formalism, (which is a particular case of the classical 
current formalism) on which the theory is based, is applicable. For these reasons the determination of 
the lifetime of the system via "pure" hydrodynamical considerations is certainly an alternative which 
deserves to be considered seriously, despite its own difhculties. 

5.1.5 Bose-Einstein correlations and pseudohydrodynamics 

As mentioned already, the initial motivation for proposing the Wigner function formalism 
for BEC was to explain why the experimentally observed second order correlation func- 
tions 6*2 were depending not only on the momentum difference q = ki — k2 but also on their 
sum K = ki + k2. However in the mean time it was shown |3^,0| that this feature follows 



from the proper application of the space-time approach in the current formalism even 
without assuming expansion. Furthermore the Wigner formalism is useful only at small q 



^^I am indebted to B. R. Schlei for this remark. 
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and cannot be applied in the case of strong correlations between positions and momenta 
while the current formalism is not limited by these constraints. As explained above the 
use of the Wigner formalism can be defended if combined with bona-fide hydrodynamics 
and an equation of state. 

This notwithstanding, besides a few real, albeit numerical, hydrodynamical calcula- 
tions, most phenomenological papers on BEC in heavy ion reactions (cf . e.g. , [ |162| | , |163| , 



||102|| , ||164|| - 1|168| , [pq] , [|69| , ||169|| - ||T78| have used the Wigner formalism without a proper 
hydrodynamical treatment, i.e. without solving the equations of hydrodynamics; hy- 
drodynamical concepts like velocity and temperature were used just to parametrize the 
Wigner source function. While such a procedure may be acceptable as a theoretical ex- 
ercise, it is certainly no substitute for a professional analysis of heavy ion reactions. This 
is a fortiori true when real data have to be interpreted^. 

As exemplified in previous sections such a procedure is unsatisfactory, among other 
things because it can lead to wrong results. 

The use of this "pseudo-hydrodynamical" approach is even more surprising if one 
realizes the fact that the Wigner formalism not only is not simpler than the more general 
current [^formalism but it is also less economical. The number of independent parameters 
^^A recent experimental paper where such a procedure is used is a good iUustration of the limits 
of pseudo-hydrodynamical models. Despite the fact that the statistics are so rich that "the statistical 
errors on the correlation functions are negligible" , the outcome of the analysis is merely the resolution 
of the ambiguity between temperature and transverse expansion velocity of the source. It is clear that 
such an ambiguity is specific to pseudo-hydrodynamics and is from the beginning absent in a correct 
hydrdynamical traetment. Moreover even this result is questionable given some doubtful assumptions 
which undcrly this analysis. To quote just two: (i) The assumption of boost invariance made in p5[ 
decouples the longitudinal expansion from the transverse one. This not only affects the conclusions drawn 
in this analysis but prevents the (simultaneous) interpretation of the experimental rapidity distribution, 
(ii) The neglect of long lived resonances which strongly influences the A factor and thus also the extracted 
radii. Of course, despite the claimed richness of the data, no attempt to relate the observations to an 

equation of state can be made within this approach. 

^■^We remind that the classical current formalism, in opposition to the Wigner function formalism, does 

not apply only to small q and to semiclassical situations. Furthermore it allows also for coherence and 

new phenomena like particle-antiparticle correlations. 
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necessary to characterize the BEC within the Wigner formahsm of ref. [[166|] is 10, i.e. it 



is as large as that in the current formahsm. However the 10 parameters of ||166|| describe 



a very particular source as compared with that of the current formalism: besides 
the fact that the correlation function source is assumed to be Gaussian, it is completely 



chaotic and it can provide only the length of homogeneity Lh |37|. For the search of 
quark gluon plasma, however, the geometrical radius R is relevant, because the energy 
density is defined in terms of R and the use of instead of R leads to an overestimate 
of the energy density 0. Furthermore the physical significance of the parameters of the 
Wigner source is unclear if the Gaussian assumption does not holdQ. Not only is there 
no a priori reason for a Gaussian form, but on the contrary, both in nuclear and particle 
physics as well as in quantum optics, there exists experimental evidence that in many 
cases an exponential function in |g| is at small q a better approximation for C2 than a 
Gaussian. Furthermore, in the presence of coherence, no single simple analytical function, 
and in particular no single Gaussian is expected to describe C2. This is a straightforward 
consequence of quantum statistics. 

Given the fact that good experimental BEC data are expensive both in terms of 
accelerator running time and man-power, the use of inappropiate theoretical tools, when 
more appropiate ones are available, is a waste which has to be avoided. For the reasons 
quoted above we will not discuss in more detail the numerous and sometimes unnecessarily 
long papers which use pseudo-hydrodynamical methods. 

Concluding remarks. In a consistent treatment of single and double inclusive cross 
sections for identical pions via a realistic hydrodynamical model, resonances play a major 



^"^For each value of K a Gaussian ellipsoid is described by three spatial extensions, one temporal exten- 
sion, three components of the velocity in the local rest frame and the three Euler angles of orientation. 



^^To consider such an approach as "model independent" [165|, |179|| is misleading. 

^^One may argue that the "length of homogeneity" Lh defined in terms of the Wigner function is 

a particular case of the correlation length L defined in terms of the correlator (cf. section 4.3). While Lh 

is always limited from above by R, L can be either smaller or larger than R. 



^^Even if a Gaussian form would hold for directly produced pions, resonances would spoil it [ p6[ , |155| 
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role leading to an increase of effective radii of sources. Effective longitudinal radii are 
more sensitive to the presence of resonances than transverse ones. /^From the hydrody- 
namical treatment we learn that the hadron source (the real fireball) is represented by 
a very complex freeze-out hypersurface (cf. ref. ||155|| ). The longitudinal and transverse 



extensions of the fireball change dynamically as a function of time, rather than show up in 
static effective radii. Thus, the interpretation of BEC measurements is also complicated. 
This is a theoretical task and if only for this reason experimentalists would be well advised 
to include in their teams also qualified theorists. 

For heuristic applications, when no quantitative comparison with data is intended, be- 
sides the current formalism, analytical approximations of the equations of hydrodynamics 
can be useful, because they allow a better qualitative understanding of hydrodynamical 
expansion. However when a quantitative interpretaion of experiments is intended and in 
particular a connection with the equation of state is looked for, the only recommendable 
method is full-fledged hydrodynamics. 

5.1.6 Photon correlations, pseudo-hydrodynamics and pseudo-Wigner for- 
malism 

Photon correlations have been investigated within the context of quark-gluon plasma 
search, since they present certain methodological advantages as compared with hadrons. 
While experimentally genuine photon BEC in high energy heavy ion reactions have not 
yet been unambigously identified, because of the strong 7r° background and the small 
cross sections for photon production, there are several theoretical studies devoted to this 
topic. The advantage of photon BEC resides in the fact that, while correlations between 
hadrons are influenced by final state interactions, photon correlations are "clean" from 
this point of view. For high energy physics photons present another important advantage 
due to the fact that they can provide direct information about the early stages of the 
interaction where quarks and gluons dominate and hadrons have not yet been created. In 
particular photon BEC contain information about the lifetime of the quark-gluon plasma 
[ 180 |, P^ . Among other things it was argued e.g. in ||180|| and confirmed in 
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89| using a more correct formalism, (cf. below), that the correlation function C2 in the 
transverse direction can serve as a signal for QGP as it is sensitive to the existence of a 
mixed phase. 

Unfortunately, some of these studies (refs.p4|-||87|, [ p.8CI|| ) besides suffering from the 
more general disease of pseudo-hydrodynamics, use an input formula for the second order 
BEG, which is essentially incorrect even within the Wigner-formalism (cf. section 4.6). 



Besides this, some approximations made in |Q e.g. are inadequate. This question was 
analyzed in more detail in ||8^ where it was found that only some of the results of ref. 
( |^4|| ) survive a more critical analysis. 

The formula for the two-particle inclusive probability used in (|84]-[B7],|180[) reads 



P2(ki,k2) = j d^xi j d^X2 g (xi, ki) g{x2,k2) [1 + cos{{ki - k2){xi - X2))], (5.7) 
while ref.fiSPIl uses the more correct formula 



P2(kl,k2) 



d'^xi j d!^X2 g \xi, 



g 1^x2, ^ [1 + cos((A;i-/i;2)(a;i-X2))], 

(5.8; 



(From eq. (|5.8| ) one gets for the second order correlation function eq.( 4.5S )). It is found 
in ||8^ that two rather surprising properties of the two-photon correlation function pre- 
sented in are artifacts of inappropriate approximations in the evaluation of space-time 
it was claimed that the BEC function in the longitudinal direction (a) 



integrals. In 

oscillates and (b) takes values below unity. As property (b) is inconsistent with general 
statistical bounds, it was important to clarify the origin of this discrepancy. On the 
other hand, it was confirmed in that the correlation function in the transverse direc- 
tion does exibit oscillatory behaviour in the out component of the momentum difference. 
Furthermore, in this reference a change of the BEC function in Ay from Gaussian to 
a two-component shape with decreasing transverse photon momentum was found which 
may serve as evidence for the presence of a mixed phase and, hence, as a QGP signature. 

However even after correcting the wrong input BEG formula it is questionanble whether 
the other approximations made in refs.([Q-|^^,|180[,[^) may not invalidate the above 
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result. Besides the use of a simplified hydro djTiamical solution one has to recall that 
(i) the Wigner formalism like the more general classical current formalism, is limited to 
small momenta k of produced particles (no recoil approximation) (ii) Besides this general 
limitation to small k the Wigner formalism is specifically limited to small differences of 
momenta q. (iii) Although Eq. ( [4.58|) does not suffer from the violation of unitarity dis- 
ease mentioned in section 4.6 of chapter 3, it is based on an approximation which makes it 
sometimes inapplicable for photons (cf. section 4.5). /^From the above discussion one may 
conclude that the experimental problems of photon BEC are matched by yet unsolved 
theoretical problems. 



5.2 Pion condensates 

One of the most interesting phenomena related to Bose-Einstein correlations is the effect 
of Bose condensates. The remarkable thing about this effect is that it is not specific for 
particle or nuclear physics, but occurs in various other chapters of physics, in particular in 
condensed matter physics like superconductivity and superfluidity. Moreover, recently not 
only the condensation of a gas of atoms has been experimentally achieved ||185|| but the 



quantum statistical coherence of these systems has been experimentally proven through 
Bose-Einstein correlations [^]. The proposal to use BEC for the detection of condensates 
was made a long time ago [^. The more recent developments in heavy ion reactions 
have made this subject of current interest. Already in experiments at the SPS (e.g., 
Pb + Pb at Ebeam = 160 AGeV) secondary particles are formed at high number densities 
in rapidity space ||181|| and in future experiments at RHIC and the LHC one expects to 
obtain even higher multiplicities on the order of a few thousand particles per unit rapidity. 
If local thermal (but not chemical) equilibrium is established and the number densities are 
sufficiently large, the pions may accumulate in their ground state and a Bose condensate 
may be formed^. A specific scenario for the formation of a Bose condensate, namely, the 



decay of short-lived resonances, was discussed in Ref. ||183|| where conditions necessary 



^Recently a different type of pion condensate, the disordered cliiral condensate, has received much of 



attention in the literature. It has been argued |182 that such a condensate would lead to the creation of 
squeezed states. 
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for the formation of a Bose condensate in a heavy ion colhsion were investigated. In Ref. 



183|| it was found that if a pionic Bose condensate is formed at any stage of the colhsion, 
it can be expected to survive until pions decouple from the dense matter, and thus it can 
affect the spectra and correlations of final state pions. 



In ref. ||184|| one investigated the influence of such a condensate on the single inclusive 
cross section and on the second order correlation function of identically charged pions 
(Bose-Einstein correlations BEC) in hadronic reactions for expanding sources. A hydro- 
dynamical approach was used based on the HYLANDER routine. The Bose-condensate 
affects the single inclusive momentum distributions EdN/d^k, the momentum-dependent 
chaoticities p and the Bose-Einstein correlation functions C2 only over a limited mo- 
mentum range. This is due to the fact that in a condensate there exists a maximum 
velocity (which implies also a maximum momentum difference qmax) and leads to a very 
characteristic structure in single and double inclusive spectra. 

In Fig. 18 the results of the numerical evaluations of the Bose-Einstein correlation 
functions C2 are shown for a spherically and for a longitudinally expanding source. The 
presence of a Bose-condensate of only 1% results in a decrease of the intercept by about 
15%. Furthermore due to a limited value of Qmax a part of the tail of the two-particle 
correlation functions is not affected by the pionic Bose-condensate and a peak appears. 
To what extent such peaks can be observed in experimental data depends among other 
things on the size of the source, the details of the freeze-out, width of the momentum 
distribution in the bosonic ground state, and detector acceptance. 

Plasma droplets? If the phase transition from hadronic matter to QGP and 
viceversa is of first order then one could expect the formation of a mixed phase, in which 
QGP and hadronic matter coexist. Such a mixed phase manifests itself in the hydrody- 
namical evolution of the system ||186|| and it influences among other things the transverse 



momentum distribution of photons as we have seen in section 5.1.6. It was suggested 



by Seibert ||187|| that the mixed phase could also lead to a granular structure which 
might be seen in the fluctuations of the velocity distributions of secondaries produced 
in the hadronization stage. Pratt, Siemens and Vischer ||188|| (cf. also [ |189|| proposed 



102 



subsequently that in Bose-Einstein correlations, too, one might see a signature of this 
granularity. 



6 Correlations and Multiplicity Distributions 

6.1 From correlations to multiplicity distributions 

An important physical observable in multiparticle production is the multiplicity distribu- 
tionQ P{n), i.e., the probability to produce in a given event n particles. The link between 
the multiplicity distribution P{n) and BEC is represented by the density matrix p , since 
it is the same p which appears in the definition of P{n) 

P{n) =< n\p\n > (6.1) 

and the definition of correlation functions (cf.e.g. eq. ([4.4|) ). Usually one expresses p in 
terms of the V{a) representation (cf. section 2.2 ) and by using for V{a) simple analytical 
expressions one is able to derive the most characterisitc forms of P{n) in an analytical 
form, like the Poisson or the negative binomial representation (cf. e.g. [^). 

Sometimes there exist however physically interesting cases where no analytical expres- 
sion for the multiplicity distribution P{n) exists, but instead the moments of P are given. 
^From the phenomenological point of view the approach to MD via moments presents 
sometimes important advantages because it allows the construction of an effective density 
matrix from the knwoledge of a few physical quantities like the correlation lengths and 
mean multiplicity, which in turn can be obtained from experiment (cf.e.g. ||190|| ). In the 



following we will address this aspect of the problem, the more so that in this way the link 
between correlations and MD becomes clearer. 

We start by recalling some definitions. Besides the normal moments of the MD given 

by 

<n'^ > = P{n)n'i (6.2) 

n 



^The multiplicity distribution will be denoted sometimes in the following also by MD 
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one uses frequently the factorial moments 



= {n{n-l)- ...-{n-q + l)) (6.3) 

n 

These can be expressed in terms of the inclusive correlation functions pq through the 
relation 

^q = [ duji... [ dujg p,(ki,...,kg) = ( ) (6.4) 

Jn Jn [n — q)l 

Eqs. ( |6.3| ) and ( p.4|) illustrate the fundamental fact that the inclusive cross sections pq and 
thus the correlation functions determine the moments of the MD, which are nothing else 
but the integrals of Pg. Although this relation between moments of MD and correlation 
functions is a straightforward aspect of multiparticle dynamics, the connection between 
MD and correlations was often overlooked. This is in part due to the fact that mea- 
surements of correlations are, for reasons of statistics and other technical considerations, 
frequently performed in different (narrower) regions of phase space than measurements of 
MD However the importance of the use of the relationship between MD and correlations 
can hardly be overemphasized, just because of the different experimental methods used in 
the investigation of these two observables. In the absence of a theory of multipaticle pro- 
duction, the form of the correlators and the amount of chaoticity are unknown and have 
therefore to be parametrized and then determined experimentally in correlation experi- 
ments. Both the parametrization and the measurements are affected by errors. Similar 
considerations apply for MD, but because of the different experimental conditions un- 
der which correlation measurements and MD measurements are made, the corresponding 
errors are different. (Cf. also the discussion of the importance of higher order correla- 
tion, section 2.2.1). Therefore from a phenomenological and practical point of view, MD 
™This is e.g the case when MD are measured with no proper identification of the particles, while BEC 



refer of course to identical particles. Thus the UA5 experiment [191|, which discovered the violation 
of KNO scaling in MD-s, measured only charged particles, without distinguishing between positive and 
negative charges. However the rapidity region accessible in this experiment was much broader than the 



corresponding region in the UA-1 experiment |192] where correlation measurements were performed and 



where a distinction between positive and negative charges could be made. 
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and correlations are rather complimentary and have to be interpreted together. In the 
following we will exemplify the usefulness of this point of view. 



6.1.1 Rapidity dependence of MD in the stationary case 

The dependence of moments of multiplicity distributions P(n) on the width of the bins 
in momentum or rapidity space has been in the center of multiparticle production studies 
for the last 15 years. It got much attention after: 1) the experimental observation [|191| 
by the UA5 collaboration that the normalized moments of P{n) in the rapidity plateau 
region increase with the width of the rapidity window Ay; 2) the proposal by Bialas and 



Peschanski |p.93|| that this behaviour, which at a first look was power-like may refiect 



"intermittency" , i.e. the absence of a fixed scale in the problem, which could imply that 
self-similar phenomena play a role in multiparticle production. 

However soon after this proposal was published it was pointed out in ||194|| that the 



quantum statistical approach, presented in section 2.2 and which implies a fixed scale 
0, predicts a similar functional relationship between the moments of MD and Ay. In 



Fig. 19 from | 194 | some examples of this behaviour are plotted and compared with 
experimental data. For small —6y the (semilogarithmic) plot can be approximated by 
a power function as indicated by the data. Recalling that the QS formalism applies to 
identical particles, it follows that BEC could be at the origin of the so called intermittency 
effect. This point of view was corraborated subsequently also in ||195|| and ||196|| and was 



confirmed experimentally by the observation that the "intermittency" effect is strongly 
enhanced when identical particles are considered Q and /or when studied in more than 
one dimension For a more recent very clear confirmation of this point of view cf. the 
studies by Tannenbaum [197|. Further developments related to "intermittency" will be 



""^This scale is the correlation length ^ of eqs. ( |2.24| ) or ( |2.25| ). 

^^Thc UA5 data refer to a mixture of equal numbers of positive and negative particles; this dilutes the 
BEC effect. 

^■^This last observation also supports the idea that BEC is the determining factor in intermittency be- 
cause the integration over transverse momentum implied by a a one dimensional y investigation diminishes 
the BEC. 
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discussed in section 6.3. 



6.1.2 Rapidity dependence of MD in the non-stationary case 

The assumption of translational invariance in rapidity permitted to apply the quantum 
optical formalism, in which time has the analogous property, to MD and led to a simple 
interpretation of the observed broadening of the MD with the decrease of the width of the 
rapidity window in high energy reactions. However stationarity in rapidity is expected to 
hold only in the central region (and only at high energies). Indeed experimental data on 
proton(antiproton)-proton collisions in the energy range 52 < ^/s < 5A0GeV show that if 
one considers shifted rapidity bins along the rapidity axis the MD in these bins depend on 
the position of the bin: in the central region the MD is broader and can be decribed by a 
negative binomial distribution while in the fragmentation region it is narrower and can be 



described by a Poisson MD ||198|| . The mean multiplicity in the central region increases 
faster (approximately like s^^^) with the energy than that in the fragmentation region. 
This was interpreted in [|198|| as possible evidence for the existence of two sources, one of 



chaotic nature localized in the central region and another coherent in the fragmentation 
region. This interpretation is in line with the folklore that gluons which interact stronger 
(than quarks) form a central blob which may be equilibrated, while the fragmentation 
region is populated by througoing quarks associated with the leading particles. 

A few years later the NA35-collaboration |pO(]| ] measured BEC in ^^O — Au reactions at 
200 GeV/nucleon in a relatively broad y region and found evidence for a larger and more 
chaotic source in the central rapidity region, and a smaller and more coherent source in the 
fragmentation region. It was then natural to correlate (cf. ref.[^) the two observations, 
i.e. that refering to MD and that refering to BEC. Taken together the credibility of 
the conjecture made in ref . ||198|| is strongly enhanced, just because we face here different 



physical observables and different experiments, each with its own specific corrections and 
biases. Moreover it was pointed out in |2^ that experiment [pOOj did not necessarily imply 



that the two sources were independent, but could also be interpeted as due to a single 



*For a microscopic interpretation of this effect in terms of a partonic stochastic model, cf. ref. [ 199 1. 
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partially coherent source. 

Indeed consider in a simplified approach as used e.g. in quantum optics a superposition 
TT of two fields one coherent denoted by vr^ and another chaotic denoted by tTch so that 
TT = 7Tc{k^±^ + TTch{k^±^). whrc Q± = k^l^ — k^f' . Assuming boost invariance the correlator 
depends only on k±^ and we have for the second order correlation function 

where the transverse "radius" i?^ plays the role of the correlation length and p = 

<kch>P> 

kcP + <|7rch>P>- 

Assume now that the chaoticity p is rapidity dependent so that in one rapidity region, 
denoted by (A), p{A) ^ 1. In that region then the third term in the equation (|6.5|) 
dominates, i.e. 

l+p'(A)e-«i^i(^) (6.6) 



In the parametrization used in ||20CI|| , according to which we have two independent sources, 
this suggests that the effective radius of source A is R^. Conversely, for the more coherent 
region denoted by B, p{B) is small and C2 reads 

C2 ^ 1 + 2p{B){l - p(5))e~<3i^i(^)/2 (6.7) 

with an effective radius R±{B) / \/2 < R±(A). 

This corresponds qualitatively to the observations made in ref. [pOCI|| . Unfortunately 
these observations have not yet been confirmed by another, independent experiment so 
that the reader should view these considerations with prudence In any case they 
prove the usefulness of a global analysis which incorporates both BEC and MD. On the 
other hand, a dedicated simultaneous investigation of the rapidity dependence of these 
two observables appears very desirable. 



''Spor various caveats concerning the analysis of BEC and MD data cf. refs.( [^00| ,|191|,|198|,|2g 
20l|). 
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6.1.3 Energy dependence of MD and its implications for BEC; long range 
fluctuations in BEC and MD 



This subject has been discussed recently in [p02|| . Besides the rapidity dependence, the 



dependence of MD on the center of mass energy of the collision, ^/s, constitutes an 
important topic in the study of high energy multiparticle production processes. This 
energy dependence is usually discussed in terms of the violation ||191|| of KNO scaling ||203| . 



KNO scaling implies that the normalized moments {ri"^) / {n)"^ are constant as a function 
of s (for high energies, i.e. large (n), these moments coincide with the normalized factorial 
moments). For charged particles it turned out that while KNO scaling is approximately 
satisfied over the range of ISR energies (20 GeV < y/s < 60 GeV), it is violated if one 
goes to SPS-CoUider energies (200 GeV < y/s < 900 GeV), i.e. one finds a considerable 
increase of multiplicity fiuctuations with increasing energy. In the following we will show 
how BEC can be used to to understand the origin of this s-dependence. To do this, one 
needs to distinguish between long range (dynamical) correlations (LRC) and short range 
correlations (SRC)[^. For identical bosons, one important type of SRC are BEC which 
refiect quantum statistical interference. In addition, there exist dynamical SRC like final 
state interactions, which however are quite difficult to be separated from BEC. Up to 1994 
one usually assumed that LRC do not play an important part in BEC measurements (see 
however ||195|| , [|196|| ) in the sense that for Q > 1 GeV the two-particle correlation functions 
do not significantly exceed unity. 

However in ref. ||202|| evidence, based on observational data, was presented showing that 



this is not the case and that new and important information about LRC is contained in 
the BEC data obtained by the UAl-Minimum-Bias Collaboration ||192|. In principle, the 



observed increase of multiplicity fiuctuations with could be due to a change of the SRC 
as seen in BEC, i.e. of the chaoticity and radii/lifetimes. This possibility was discussed in 
||207|| but could not be tested because of the lack of identical particle data for multiplicity 
distributions at Collider energies. At that time only the UA5 data ||191|| for multiplicity 



''^The importance of making this distinction was pointed out, among other things, in | p04| , |205[ . In 
06| rapidity correlations were measured for events at fixed multiphcity in order to get rid of the effect 
of LRC. 
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distributions of charged particles were available. Furthermore, up to this point, the effect 
of LRC had only been studied in terms of two-particle correlations as a function of rapidity 
difference, i.e. in one dimension. With the advent of the newly analyzed UAl data | |192| 
for identical particles in three dimensions (essentially in = —{ki — k2Y) this situation 



changed. The analysis of |P02|] led to the conclusion of the existence in BEC data of long 
range fluctuations in the momentum space density of secondaries and to the realization 
that the increase with energy of multiplicity fluctuations is to a great extent due to an 
increase of the asymptotic values of the m-particle correlation functions C^**^™"^*', i.e. their 
values in the limit of large momentum differences where BEC do not play a role Q. 



In ref. ||192|| , the UAl collaboration presents the two-particle correlation of negatively 
charged secondaries as a function of the invariant momentum difference squared = 
(ki — k2)^ — {El — E2Y. The data (cf. Fig. 20) have two unusual features: (1) at large 
the correlation function saturates above unity, and (II) at small it takes on values above 
2. The higher order correlation functions also exhibit property (I) ||192|| ^ By comparing 
the asymptotic values of the correlation functions at large momentum differences C^'^™'^^' 
[m = 2, ...,5), with the normalized factorial moments, (pm = {n{n — l)-...{n — m+l))/ (n)"^ 
in the momentum space region|?/| < 3, A;^ > 0.15 GeV one finds that the contribution of 
the BE interference peak to the moments is negligible for such large rapidity windows. 
Herefrom one concludes in ||202|| that (I) indicates the presence of LRC in the momentum 
space density of secondary particles and that it is quite plausible (cf. below) that (II) has 



to a great extent the same explanation. We sketch here the arguments of ref. [ p02|| . 

In general, LRC may be related to fluctuations in impact parameter or inelasticity, or 
fluctuations in the number of sources. In what follows, let us label these fluctuations by 
a parameter a. The m-particle Bose-Einstein correlation function at a fixed value of a is 
given by 

r-i n 1 I ^ Pm(ki, km a) m o\ 

C^{ki,...,k^\a) = — -- (6.8) 



^^The effect of LRC on BEC was discussed in ref. [|208| , where several models specific for nucleus- nucleus 

collisions were considered, but at that time no evidence for this effect could be found. 

''^For higher order correlations the equivalent of property (II) is Cm > m)-- (cf. also below). The values 

of Cm(k, k) for TO > 2 are apparently not yet available. 
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where pi^ki, ...,ke\a) are the ^-particle inclusive distributions. 

The fluctuations in a are described by a probability distribution h{a) with 



da h{a) = 1. (6.9) 
If the experiment does not select events at fixed a, the measured inclusive distributions 

are 

Pm(ki, ...,k^) = (p„(ki, ...,k„|a)) , (6.10) 
where the symbol (...) denotes an average over the fluctuating parameter a, i.e. 

(X(a)) ^ fda h{a) X{a) (6.11) 



The m-particle correlation function at the intercept reads 

C^(k,...,k) = ml (6.12) 

where the symbols <> refer to averaging with respect to h{a). At large momentum 
differences one has 



C^(ki,...,k^) ^ f-f = Cr™^*- for |k,-k 



■J I 



oo 



{i 7^ j, hi = 1, ■■■m), (6.13) 

i.e., the m-particle correlation functions can have intercepts above m! and saturate at 
values above unity for large momentum differences. 

The most obvious candidate for the function h{a) is the inelasticity distribution which 
describes the event to event fluctuations of the inelasticity K. With the identification 
a = {n{K)), where {n{K)) is the mean multiplicity at inelasticity K, one obtains from 
the above considerations a first "experimental" information about this important physical 
quantity at collider energies. Previous experimental information about this distribution 
was derived in ||209|| from the data of ref. ||210|| at ^/s 17 GeV. 



The conclusions obtained in ||202|| about LRC are based among other things on the 
different normalizations used in different experiments. Some tests related to these con- 
clusions are proposed: 
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• Analysis of BEC at lower energies (NA22 range) as well as at ^/s — 1800 GeV 
with the same normahzation as that used by the UAl Collaboration. The values 
of C2 at > 1 GeV"^ and possibly also at very small obtained in this way 
should exceed those obtained with the fixed multiplicity normalization in the same 
experiments. The following inequalities for C2{Q^ > 1 {GeV)"^), should be observed 
if this normalization is used: 

C2(NA22) < C2(UA1) < C2(Tevatron). 

• Analysis of BEC at UAl energies with the same normalization as that used so far by 
the NA22 and Tevatron groups (fixed multiphcity) . The enhancement of C2 at large 
Q (and possibly also at small Q) observed so far, should disappear to a great extent. 

6.2 Multiplicity dependence of Bose-Einstein correlations 

The operators for the field (intensity) and number of particles do not commute. This 
means that measurements of "ideal" BEC can be performed only when no restriction on 
the multiplicity n, which fluctuates from event to event, is made. In practice, however, 
very often such restrictions are imposed, either because of technical reasons or because 
of theoretical prejudices. To the last category belong considerations imposed by the 
search for QGP in high energy heavy ion reactions. Thus one expects that by selecting 
events with n > rimin, where rimin is in general an energy dependent quantity, one gets 
information about the interesting "central" collisions. 

Another reason why multiplicity constraints are of practical importance for QGP ex- 
periments is the need to compare various QGP signals in a given event and at the same 
time determine for that event the radius, lifetime and chaoticity of the source, among 
other things in order to be able to estimate the energy density achieved in that event. 
This means that for QGP search it is interesting to perform interferometry measurements 
for single events, which of course have a given multiplicity. 

For these reasons the investigation of the multiplicity dependence of BEC constitutes 
an important enterprise which we shall address in the following section. 
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6.2.1 The quantum statistical formalism 

Correlation functions defined in quantum statistics and used in quantum optics refer to 
ensemble averages of intensities / of fields vr, where 

/(kx,y) = |7r(k^,|/)|2. (6.14) 
The total multiplicity n of identical particles over a given phase space region is given 

by 

n = JdklJ dy\7r{k^,y)\^ (6.15) 

Both the field 7r(k^, y) and the intensity I(k±, y) are stochastic variables. Averaging over 
an appropriate ensemble, we get the mean total multiplicity 

<n>= J dklj dy <\n{k^,y)\'^ > . (6.16) 



In pil|| the dependence of BEC within the QS formalism on the total multiplicity n 
and on n^m was investigated and it was found that the size of this effect is (especially at 
low < n >) surprisingly large and must not be ignored, as had been done before. Both the 
Umin constraint and the < n > constraint lead to a decrease of the correlation function C2 
at fixed yg = yi~ y2, i-e. to an antibunching effect. The last effect can be approximated, 
except for very large n and small yg by a simple analytical formula 

where denotes the correlation function at fixed n and /2 is the reduced factorial 
moment. 

These results show that BEC parameters like radii, lifetimes, and chaoticity do depend 

on the particular experimental conditions under which the measurements are performed. 

It is worth mentioning that a multipHcity dependence of BEC was observed experimentahy for p — p 



and a — a reactions at Ecm — 53 GeV and 31 GeV respectively already in [212|. There it was found 



that the transverse radius increases with the multiplicity of charged particles ficharged- This effect was 



interpreted by Barshay [213| to be a consequence of the impact parameter dependence. The same effect 
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Herefrom the name "intensity interferometry" . 
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was seen in heavy ion reactions |20C|| and got a similar interpretation in [176| q[ The interpretation in 
terms of impact parameter dependence could be checked directly in heavy ion reactions since here the 
impact parameter can be determined on an event by event basis. 



Another mechanism for the increase of radius with multiplicity was proposed by Ryskin |214|. He 
pointed out that in high multiplicity events, which are associated by many authors to large transverse 
momenta of partons and thus to a regime where perturbative QCD applies, one expects that the size of 
the hadronization region should increase with the multiplicity like y^. 

A related topic is the dependence of BEC on the rapidity density d = An/ Ay, which 



has also been observed experimentally in p — p reactions at the CERN SPS collider [pi5 
and at the Fermilab tevatron | |216| |. Using a parametrization 

C2 = 1 + Xexpi-R^Q^) (6.18) 

where Q is the invariant momentun transfer, it was found that R increased with d while A 
decreased with d. This last observation is compatible with the results from |[^11| 



a more 



quantitative comparison would be possible only if, among other things, the data were 
parametrized in a way more consistent with quantum statistics. Thus the four momentum 
difference Q is not an ideal variable for BEC (cf. section 2.1.5) and the coherence effect 
has to be taken into account as outlined in section 2.2 and not by the simple empirical A 
factor. 

6.2.2 The wave function formalism; "pasers"? 

The dependence on total multiplicity of BEC was investigated also within the wave func- 
tion formalism. In section 2.1 where the GGLP theory was presented it was pointed out 
that the wave function formalism may be useful for exclusive processes or for event gen- 
erators. Indeed, in a first approximation, the wave function ipn of a system of n identical 
bosons e.g. can be obtained from the product of n single particle wave functions tpi by 
symmetrization. Then the calculation of 0^2^ is in principle straightforward and follows 
the lines of GGLP. 
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The approach of |176| combines simplified (1-d) hydrodynamics with a multiple scattering model, 



which also exploits the impact parameter dependence. 
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However when the multiphcities n become large (say n > 20) the exphcit symmetriza- 
tion of the wave function formahsm becomes difficult. This lead Zajc | |217| | to use numerical 
Monte Carlo techniques for estimating n particle symmetrized probabibilities, which he 
then applied to calculate two-particle BEC. He was thus able to study the question of 
the dependence of BEC parameters on the multiplicity. For an application of this ap- 
proach to Bevalac heavy ion reactions, cf. ||218|] . Using as input a second order BEC 
function paramerized in the form ( |6.18| ) Zajc found, (and we have seen above that this 
was confirmed in |pil|| ,) that the "incoherence" parameter A decreased with increasing n 

n 

However Zajc did not consider that this effect means that events with higher pion 
multiplicities are denser and more coherent. On the contrary he warned against such an 
interpretation and concluded that his results have to be used in order to eliminate the 
bias introduced by this effect into experimental observations 0. 

This warning apparently did not deter the authors of [^19|| and [|22(]|| to do just that. 



Ref. |pl9| | went even so far as to deduce the possible existence of pionic lasers from con- 
siderations of this type 0. 

Ref. ||219| | starts by proposing an algorithm for symmetrizing the wave functions 
which presents the advantages that it reduces very much the computing time when using 
numerical techniques, which is applicable also for Wigner type source functions and not 
only plane wave functions, and which for Gaussian sources provides even analytical results. 
Subsequently in ref. |p21|] wave packets were symmetrized and in special cases the matrix 
density at fixed and arbitrary n was derived in analytical form. This algorithm was then 
applied to calculate the infiuence of symmetrization on BEC and multiplicity distributions. 

®^In | |217t the clumping in phase space due to Bose symmetry was also illustrated; 

*^The same interpretation of the multiplicity dependence of BEC was given in |211|. In this reference 

the nature of the "fake" coherence induced by fixing the multiplicity is even clearer, as one studies there 

explicitly coherence in a consistent quantum statistical formalism. 

^■^The philosophy of rcf. [219 is reflected e.g. in statements like "multiparticle interference can enhance 

the emission of pions" . This language is apparently motivated by the practice of event generators where 

BEC (or in general quantum statistical effects) are "added" to the usual schemes at the end of the 

calculation. That this does not reflect the true sequence of the physical processes considered is obvious. 
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As in ||217|| it is found in ||219|| that the symmetrization produces an effective decrease of 
the radius of the source, a broadening of the multiphcity distribution P{n) and an increase 
of the mean multiphcity as compared to the non-symmetrized case. What is new in [pl9|| 



is (besides the algorithm) mainly the meaning the author attributes to these results. 

In a concrete example Pratt considers a non-relativistic source distribution S in the 
absence of symmetrization effects: 

^(^' ^) = (2.i?ir)3/2 - 2^) ^'-''^ 

where 

ko/T = ky2A^ (6.20) 

Here T is an effective temperature, R an effective radius, m the pion mass and A a 
constant with dimensions of momentum. 

Let ?7o aiid rj be the number densities before and after symmetrization, respectively. 
In terms of S{k,x) we have 

Vo = J S{k,x)d^kd^x (6.21) 

and a corresponding expression for 7] with 5* replaced by the source function after sym- 
metrization. Then one finds [ [21 9] ] that r] increases with tjq and above a certain crtitical 
density rj^"^, rj diverges. This is interpreted by Pratt as "pasing" . 

The reader may be rightly puzzled by the fact that while t] has a clear physical signif- 
icance the number density tiq and a fortiori its critical value have no physical significance, 
because in nature there does not exist a system of bosons the wave function of which is 



not symmetrized. Thus contrary to what is alluded to in ref. | 219 |, this paper does not 
does address really the question how a condensate is reached. Indeed the physical factors 
which induce condensation are, for systems in (local) thermal and chemical equilibrium 
0, pressure and temperature and the symmetrization is contained automatically in the 
distribution function 

f = (6.22) 



*For lasers the determining dynamical factor is among other things the inversion. 
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in the term -1 in the denominator; E is the energy and /i the chemical potential. 

To realize what is going on it is useful to observe that the increase of tjq can be achieved 
by decreasing R and/or T. Thus tiq can be substituted by one or both of these two physical 
quantities. Then the blow-up of the number density r] can be thought of as occuring due 
to a decrease of T and/or R. However this is nothing but the well known Bose- Einstein 
condensation phenomenon and does not represent anything new. While from a purely 
mathematical point of view the condensation effect can be achieved also by starting with 
a non- symmetrized wave function and symmetrizing it afterwards "by hand" , the causal 
i.e. physical relationship is different: one starts with a bosonic i.e symmetrized system 
and obtains condensation by decreasing the temperature or by increasing the density of 
this bosonic system. 

Another confusing interpretation in | |219 | relates to the observation made also in [ ^17]| 
that the symmetrization produces a broadening of the multiplicity distribution (MD). 
In particular starting with a Poisson MD for the non-symmetrized wf one ends up after 
symmetrization with a negative binomial. While Zajc correctly considers this as a simple 
consequence of Bose statistics, ref. [pl9|| goes further and associates this with the so called 
pasing effect. That such an interpretation is incorrect is obvious from the fact that for 
true lasers the opposite effect takes place. Before "condensing" i.e. below threshold 
their MD is in general broad and of negative binomial form corresponding to a chaotic 
(thermal) distribution while above threshold the laser condensate is produced and as 
such corresponds to a coherent state and therefore is characterized by a Poisson MD. 
Last but not least the fact that this broadening increases with n is not due, as one might 
be mislead to believe from [pi9|| , p20|] , to the approach to "lasing criticality" , but simply 
to the fact that the larger n the larger the number of independent emitters is and the 
better the central limit theorem applies. This theorem states (cf. section 2.2) that the 
field produced by a large number of independent sources is chaotic. 

Finally a terminological remark appears necessary here. We believe that names like 
"paser" or pionic laser used in the papers quoted above are unjustified and misleading. 
The only characteristic which the systems considerd in these papers possibly share with 
lasers is the condensate property i.e. the bunching of particles in a given (momentum) 
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state. However lasers are much more than just condensates; one of the main properties 
of lasers which distinguishes them from other condensates is the directionality, a problem 
which is not even mentioned in the "paser" literature. Q 

To conclude this excursion into what in our opinion is a surrealistic interpretation of 
multiparticle correlations, it is hard to find new physics in the "paser" papers, although 
they probably represent progress in the mathematical problem of symmetrization of wave 
functions or wave packets[^. 



6.3 The invariant Q variable in the space-time approach; higher 
order correlations; "intermittency" in BEC? 

The issue of apparent power-like rapidity dependence of moments of MD was discussed 
in section 6.1.1 where it was pointed out that this dependence could in principle be 
understood within the QS formalism without invoking the idea of intermittency. However 
this was not the end of the story because: 1) it was observed ||192|| that the power like 
behaviour extends also to BEC data (in the invariant variable Q). This was surprising 
because up to that moment BEC data could usually be fitted by a Gaussian or exponential 
function, albeit these data did not extend to such small Q values as those measured in 
192|; 2) Bialas |p27|| (cf. also |228|, |229|) proposed that the source itself has no fixed 



size, but is fluctuating from event to event with a power distribution of sizes. Since the 
measurements made in [p.92|| were in the invariant variable Q rather than Ay, one had to 
understand whether the QS approach which implies fixed scales does not lead to a similar 
behaviour in Q. 



Refs. ||230|| , ||231|| addressed this question and proved that, indeed, by starting from 
a space-time correlator with a fixed correlation length and a source distribution with 
a fixed radius, one gets after integrations over the unobserved variables a correlation 
function which is power-like in a limited Q range. 
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For a model of directional coherence, not necessarily related to pion condensates, cf.[222|; experimen- 



tal hints of this effect have possibly been seen in [223 

*^For another investigation of the effect of symmetrization on the single inclusive cross section cf. [224 



and [225|. In [225[ second order correlation functions are also considered. (Cf. also ref.[226[ for this topic.) 
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In the following we shall sketch how this happens. The two-particle Bose-Einstein 
correlation function is defined as 

C2(k„k,)= (6.23) 

where pi(k) and p2(ki,k2) are the one- and two-particle particle inclusive spectra, res- 
pectively. 

The two-particle correlation function projected on is 

with the integrals 

Xn(g')= / duj, I duj2 S + (A;f - k!^f] pi(ki)pi(k2) 

UQ')= J duJ, I duJ2 5 + (fcf - ^2^)2] (p2(ki,k2) - pi(ki)pi(k2)) (6.25) 

where du = d^k/ {2txY2E is the invariant phase space volume element. Two types of 
sources, a static one and an expanding one, were considered. The parameters employed 
are physically meaningful quantities in the sense that they give the lifetime, radii and 
correlation lengths of the source. This is not the case for the ad-hoc parametrization of 
C2(g^), e.g., with a Gaussian 



C2(g2) = 1 + Age-^^^' (6.26) 

To illustrate the behaviour of the correlation function as a function of g^, one applied 
the formalism to describe UAl data in the phase space region \y\ < 3.0 and k_^_ > 150MeV. 
In | 230 ], P^ C2(g^) was calculated for the static and for the expanding source by Monte 



Carlo integration, (for the static source, approximate analytical results could be obtained 



only for \y\ > 1.5, where they agree with the numerical results [p31j| ). Fig. 21 shows the 

results of fits to the UAl dataQ for the static source (dashed line) and for the expanding 
^^The C2((5^) data of | [192[ were normalized so that at large Q^, C2{Q^) ~ 1- An explanation for the 
experimental observation of |192| that C2{Q^) exceeds by a multiplicative factor of ~ 1.3 both the upper 



and lower "conventional" limits of 2 and 1, respectively, is discussed in ref. |202]. 
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source (solid line), which were obtained under the assumption of a purely chaotic source, 
= 1 (the data are consistent with an amount of < 10% coherence^, but the sensitivity 
to po was not sufficient to further constrain the degree of coherence within these limits). 
For comparison, the result of a power law fit (dotted line) as suggested in ref. ||192|| 

is also plotted. 

One finds that the data can be well described both with the static and the expanding 
source model with reasonable values for the radii, lifetime and correlation length. 

The results of above show that the power-like behaviour of C2 (the same holds for 
higher order correlations) can be reproduced by assuming a conventional space-time source 



with fixed parameters, i.e., without invoking "intermittency" |_|. This conclusion of [ 23C ] 
is strengthened by an explicit consideration of resonances in p31|| . 

The advantages of the QS formalism as compared with the wave function formalism 
emerge clearly also in the problem of higher order correlations. In the QS formalism higher 
order correlations are treated on the same footing as lower order ones and emerge just as 
consequences of the form of the density matrix. Therefore questions found sometimes in 
the literature like "what is the infiuence of higher order correlations on the lower one" 
do not even arise in QS and in fact do not make sense. We emphasized in section 2.2.1 
the importance of higher order correlations for the phenomenological determination of the 
form of the correlation function and this applies in particular when a single variable like 
Q is used instead of the six independent degrees of freedom inherent in the correlator. 
For this reason the space-time integration at fixed Q has been used recently [^| also 
in the study of higher order correlations and applied to the NA22 data [ p33|| . It was 
found among other things that an expanding source with fixed parameters as defined in 



^The older UA-1 data [bTI were limited to larger values. Furthermore the quantum statistical 



interpretation |p2| of these data is different from that in [230|, which is based on space-time concepts. 



This explains the different values of chaoticity obtained in [ ^30[ , [ 231 1 on the one hand and in refs. , |Q 
on te other hand 



'A similar point of view is expressed in |232| for the particular case of bremsstrahlung photons. 
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chapter 4 can account for the data up to and including the fourth order, confirming in 
a first approximation the Gaussian form of the density matrix. The fact that previous 
attempts in this direction hke | p3^| and [^] met with difficulties may be due to the fact 
that in these two experimental studies the QO formalism in momentum space was used 
and possibly also because no simultaneous fit of all orders of correlations was performed, 
as was the case in |^ . 



7 Critical discussion and outlook 

For historical reasons related to the fact that the GGLP effect was observed for the first 
time in annihilation at rest when (almost) exclusive reactions were studied, the theory of 
BEG was initially based on the wave function formalism. This formalism is not appropiate 
for inclusive reactions in high energy physics, among other things because it yields a cor- 
relation function which depends only on the momentum difference q and not also on the 
sum of momenta, it does not take into account isospin, and cannot treat adequately coher- 
ence. This last property is essential as it leads to one of the most important applications 
of BEG, i.e. to condensates. Furthermore this property strongly affects another essential 
application of BEG, namely the determination of sizes, lifetimes, correlation lengths and 
correlation times of sources. The adhoc parametrization of the correlation function under 



the form ( |2.8| ) where A is supposed to take into account (in) coherence is unsatisfactory. 
An improvement based on an analogy with quantum optics leads to an additional term 
(with the same number of parameters). A more complete and more correct treatment of 
BEG is provided by the space-time formalism of classical currents. 

Almost all studies of BEG assume a Gaussian density matrix. Possible deviations from 
this form could and should be be looked for by studying higher order correlations. 

The classical current approach is based on an exact solution of the equations of quan- 
tum field theory and it can be considered at present the most advanced and complete 
description of BEG. It introduces a primordial correlator of currents which is character- 
ized by finite correlation lengths and correlation times. The geometry of the source is an 
independent property of the system and it is here that the traditional radii and lifetimes 
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enter. The form of the geometry and of the correlator is not given by the theory and it is 
up to experiment to determine these. 

For a Gaussian density matrix the space-time approach within the classical current 
formalism leads to a minimum of ten independent parameters; these include geometrical 
and dynamical scales as well as the chaoticity. Phenomenologically these scales can be 
separated only by considering simultaneously single and double inclusive cross sections. 
Experimentally this separation as well as the determination of all parameters has not yet 
been performed and constitutes an important task for the future. This should be done not 
only for particle reactions but also for heavy ion reactions, as an alternative to the pseudo- 
hydrodynamical approach which cannot provide this separation. It would also be desirable 
to extend the parametrization for an expanding source by dropping the assumption of 
boost invariance. 

Besides the conventional +- 1- or pion correlations there exist also H — correlations, 

which become important for sources of small lifetimes. These "surprising" field theoretical 
effects represent squeezed states, which unlike what happens in optics, appear in particle 
physics "for free" . These effects can be used to investigate the difference between classical 
and quantum currents. Their detection constitute one of the most important challenges 
for future experiments. 

BEG are influenced by final state interactions. Goulomb final state interactions do 
not play a major part except at very small q. These small values have apparently not 
yet been reached in experiment and it is questionable whether they will be reached in 
the foreseable future. In heavy ion reactions this is due to the large values of radii which 
multiply q and even for typical scales of 1 fm the present resolution of detectors is not 
sufficient to make this effect very important. Nevertheless since Goulomb corrections have 
been studied so far only within the wave function formalism and usually by applying the 
Schrodinger equation, it would be desirable to extend this study by considering the Klein- 
Gordon equation which is more appropiate for mesons. Even more interesting would be 
to study Goulomb corrections within the classical current formalism. 

Resonances play an essential part in final state interactions and progress has been 
achieved in their understanding, in particular in heavy ion reactions, where their influence 
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has been investigated using solutions of the equations of hydrodynamics within the Wigner 
function formahsm 

BEC have been investigated in e"*" — e~, hadron-hadron and heavy ion reactions, how- 
ever a systematic comparison of resuhs using the same paramterization and normahzation 
awaits still to be doneQ. 

Correlations are intimately related to multiplicity distributions which can serve as 
complementary tools in the determination of the parameters of sources. Therefore a 
systematic investigation of BEC and multiplicity distributions in the same phase space 
region is desirable. How useful this can be has been shown by proving in this way the 
influence of long range correlations on BEC. 

BEC can be useful in heavy ion reactions and in particular for the search of quark- 
gluon plasma if one of the two conditions are satisfied: 

(A) the investigation is based on full-fledged hydrodynamics, impying the solution of 
the equations of hydrodynamics with explicit consideration of the equation of state. 

(B) the investigation is based on the classical current space-time approach. 

Case (A) has the advantage that the dependence of the equation of state on the phase 
transition may be reflected also in single inclusive cross sections and in the BEC. (So far, 
however, the sensitivity of BEC on the equation of state could not yet been proven with 



present data ||234|| .) It has however the disadvantage that its applicability is restricted 
because it is based on the Wigner function, which is a particular case of the classical 
current formalism, for small q and not too strong correlations between momenta k and 
coordinates x. 

Case (B) has the disadvantage that there is no contact with the equation of state. 
However it has the advantage it is not restricted to small q and weak correlations between 
k and x. 

Unfortunately many of the "theoretical" papers on BEC in heavy ion reactions do not 
satisfy either condition (A) or condition (B). This is the case with most of the pseudo- 



^°At present even for the same type of reaction different normalizations are frequently used and this 
can lead to apparently different results, as exemplified in the case of NA22 and UAl data. Tests for a 
better understanding and elimination of these discrepancies have been proposed (cf. section 6.1.3) 
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hydrodynamical papers which use a parametrization of the source function based on 
quahtative hydrodynamical considerations without the use of an equation of state and 
without solving the equations of hydrodynamics. This pseudo-hydrodynamical approach 
has the disadvantages of (A) and (B) but none of their advantages]^. 

/^From the above considerations it is seen that there are many usolved problems in 
the investigation of BEC, some of them of theoretical, but most of them of experimental 
nature. 

While an analysis of the experimental BEC deserves a special review, some of the 
obvious reasons for this unsatisfactory experimental situation are: 

1) Most of the BEC experiments performed so far use inadequate detectors, because 
they are not dedicated experiments but rather byproducts of experiments planned for 
other purposes. What is needed amng other things is track by track detection and im- 
proved identification of particles. 

2) Insufficient statistics. An improvement of statistics especially at small q by at least 
one order of magnitude is necessary to address some of the problems enumerated above. 

3) Incorrect or incomplete parametrizations of the correlation functions. Very often 
and in part because of 2) not all six independent variables of C2 are measured, but 
projection of these. Very popular among these projections is the relativistically invariant 
variable Q. This is not a good variable for BEC studies, because among other things 
it mixes the space and time variables in an uncontrollable way. Furthermore, in most 
parametrizations, coherence is not (or inadequately) considered. 

4) Inadequate normalizations. Practically all BEC experiments use a normalization 
procedure of the correlation function which does not correspond to its definiton. This 
definition refers to the single inclusive cross sections obtained in the same event as the 
double inclusive cross section. The use of an "uncorrelated background", instead, biases 
the results. 

^^Pseudo-hydrodynamics as compared with hydrodynamics has also the supphmentary disadvantage 
that it does not aUow a separation between geometrical radii and correlation lengths. Moreover it does 
not have even the excuse of simplicity, since the number of free parameters in the psudo-hydrodynamical 
approach is as large as that in the classical current space-time approach. 
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The solution of the problems mentioned above will make of boson interferometry what 
it is supposed to be: a reliable method for the determination of sizes, lifetimes, correlation 
lengths, and coherence of sources in subatomic physics. 

A more pedagogical presentation of the theory of Bose-Einstein correlations, which 
discusses also its quantum optical context, including a comparison between the HBT 
and the GGLP effects and between photon and hadron intensity interferometry can be 
found in the book by the author [§. Some of the most representative theoretical and 
experimental papers on EEC and which are frequently quoted within the present review 
have been reprinted in a single volume in ref. | |236| . 

Acknowledgements 1 am indebted to D. Strottman for a careful reading of the 
manuscript and for many helpful comments. 
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8 Figure captions 

Fig. 1 The GGLP experiment schematically. 

Fig. 2 Second order correlation function C2{k,k) = g'^'^^ as a function of the squeezing 
parameter r for pure squeezed states (from, ref . [p^ ) . 

Fig. 3 Second order correlation function as given by the quantum statistical formalism 
of ref. 1^ for various values of the coupling constant g for k = 0.5 (from [0). The 



parameter k is related to the chaoticity p via the relation k =< ric > / < rich >= 
1/p — 1 where p =< rich I n > and n = ric + rich is the total multiplicity. Y is the 
maximum rapidity. 

Fig. 4 The same for various values of k at g = (from [p9[| ). 

Fig. 5 Second order correlation function for various values of the "incoherence" param- 
eter A as given by the phenomenological equation ( |3.54| ). 



Fig. 6 Diagrams contributing to different terms of the two-particle correlator in the 
point-like random source model (from ref. 0). 

Fig. 7 The linear polarization vectors €iambda{k) for two photons with momenta ki and 
k2 (from ref.[§2[). 

Fig 8 Geometry of the boost-invariant source. 

Fig. 9 Intercept of two-particle correlation function in the presence of coherence and 
resonances (from ref. 



Fig. 10 Bose-Einstein correlation functions in longitudinal and transverse direction, for 
the 3-dimensional (solid lines) and the 1-dimensional calculations (dashed lines) 
(from ref. pi). 



Fig. 11 Dependence of the longitudinal and transverse radii extracted from Bose-Einstein 
correlation functions on the rapidity yx and average momentum K± of the pair. 
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As before, solid lines correspond to the 3- dimensional and dashed lines to the 1- 
dimensional results. The open circles indicate values of R\\ obtained from eq. (|5.4|) 
(fromref. pij ). 

Fig. 12 Bose-Einstein correlation functions of negatively charged pions, in longitudinal 
and transverse direction. The separate contributions from resonances are succes- 
sively added to the correlation function of direct (thermal) tt^ (dotted line). The 
solid line describes the correlation function of all vr^ (from ref. ||150|| ). 

Fig. 13 Correlation functions of all tt~ (solid lines), thermal tt~ (dotted lines) and all 
(dashed lines), for k± = and for kj_ = 1 GeV/c (from ref. | |150| ]). 

Fig. 14 Dependence of the longitudinal and transverse radii extracted from Bose-Einstein 
correlation functions on the rapidity and average momentum k± of the pair, for 
all 7T~ (solid lines), thermal it" (dotted lines) and all K~ (dashed lines). The full 
circles were obtained by substituting the value < tf{z = 0,r±) >= 2 fm/c for the 
average lifetime of the system (calculated directly from hydrodynamics by averaging 
over the hypersurface) into eq. (|5.4|), with Tj = 0.139 GeV (from ref. 



Fig. 15 TT TT Bose-Einstein correlation functions in the presence of partial coherence 



(from ref. [|15CI||) 



Fig. 16 Effective radii extracted from Bose-Einstein correlation functions as a function 
of the rapidity yt of the pair and the transverse average momentum K± of the pair 
for all pions ( from ref. ||155|| ). 



Fig. 17 Effective radii extracted from Bose-Einstein correlation functions as a function 
of the transverse average momentum of the pair for all pions compared with 
data (from ref. ||155|| ). 



Fig. 18 Two-particle BE correlation functions for a spherically and a longitudinally ex- 
panding source. The different line styles correspond to different condensate densities 



Tico compared to the thermal number densities rith (from ref. [ [184|| ). 
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Fig. 19 Normalized factorial moments $g of order q in finite (pseudo) rapidity windows 
of width 5y around ycMS = 0, plotted against 5rj (from ref. [|194[ ). 

Fig. 20 Second order correlation function for negative particles at ^/s = QSOGeV (from 
ref. p2| ). 



Fig. 21 C2{Q^) for the static source model (dashed line) and for the expanding source 
model (solid line) compared to the UAl data ||192|| . The dotted hne shows a power 
law fit. 
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